Commit 73a293dd authored by Jait Dixit's avatar Jait Dixit
Browse files

WIP: make transforms unitary

parent 93d877c3
...@@ -277,6 +277,7 @@ class Field(object): ...@@ -277,6 +277,7 @@ class Field(object):
def power_synthesize(self): def power_synthesize(self):
# check that all spaces in self.domain are real or instances of power_space # check that all spaces in self.domain are real or instances of power_space
# check if field is real- or complex-valued # check if field is real- or complex-valued
pass
# ---Properties--- # ---Properties---
......
...@@ -26,8 +26,7 @@ class FFTOperator(LinearOperator): ...@@ -26,8 +26,7 @@ class FFTOperator(LinearOperator):
if target is None: if target is None:
target = utilities.get_default_codomain(self.domain[0]) target = utilities.get_default_codomain(self.domain[0])
self._target = self._parse_domain( self._target = self._parse_domain(target)
utilities.get_default_codomain(self.domain[0]))
self._forward_transformation = TransformationFactory.create( self._forward_transformation = TransformationFactory.create(
self.domain[0], self.target[0] self.domain[0], self.target[0]
......
...@@ -135,4 +135,10 @@ class RGRGTransformation(Transformation): ...@@ -135,4 +135,10 @@ class RGRGTransformation(Transformation):
# correct for inverse fft # correct for inverse fft
Tval = self._transform.codomain.weight(Tval, power=-1, axes=axes) Tval = self._transform.codomain.weight(Tval, power=-1, axes=axes)
if self._transform.codomain.harmonic:
Tval *= 1.0 / (np.sqrt(self._transform.domain.dim))
else:
Tval *= np.sqrt(self._transform.domain.dim)
return Tval return Tval
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment