diff --git a/odl/trafos/backends/pyfftw_bindings.py b/odl/trafos/backends/pyfftw_bindings.py index 0150c754df8..75348a22563 100644 --- a/odl/trafos/backends/pyfftw_bindings.py +++ b/odl/trafos/backends/pyfftw_bindings.py @@ -203,7 +203,10 @@ def pyfftw_call(array_in, array_out, direction='forward', axes=None, else: fftw_plan = fftw_plan_in - fftw_plan(array_in, array_out, normalise_idft=normalise_idft) + if not normalise_idft and direction=='forward': + fftw_plan(array_in, array_out, normalise_idft=True) + else: + fftw_plan(array_in, array_out, normalise_idft=normalise_idft) if wexport: try: diff --git a/odl/trafos/fourier.py b/odl/trafos/fourier.py index 03af04f02ff..0ce5a6a982b 100644 --- a/odl/trafos/fourier.py +++ b/odl/trafos/fourier.py @@ -675,7 +675,8 @@ def _call_pyfftw(self, x, out, **kwargs): halfcomplex=self.halfcomplex, planning_effort=effort, fftw_plan=self._fftw_plan, normalise_idft=True) - # Need to normalize for 'forward', no way to force pyfftw + # Need to normalize for 'forward', pyfftw before version 0.13 + # does not offer a way to do this. if self.sign == '-': out /= np.prod(np.take(self.domain.shape, self.axes))