Commit de783f8d authored by Philipp Arras's avatar Philipp Arras

Use HartleyOperator

parent fa765677
......@@ -54,7 +54,7 @@ class StarbladeEnergy(ift.Energy):
h_space = position.domain[0].get_default_codomain()
else:
h_space = correlation.domain[0]
FFT = ift.HarmonicTransformOperator(h_space, position.domain[0])
FFT = ift.HartleyOperator(h_space, position.domain[0])
self._ptanh = ift.PositiveTanh()
a = self._ptanh(position)
......
......@@ -39,7 +39,7 @@ def build_starblade(data, alpha=1.5, q=1e-40, cg_iterations=500):
s_space = ift.RGSpace(data.shape, distances=len(data.shape) * [1])
data = ift.Field.from_global_data(s_space, data)
h_space = s_space.get_default_codomain()
FFT = ift.HarmonicTransformOperator(h_space, s_space)
FFT = ift.HartleyOperator(h_space, s_space)
binbounds = ift.PowerSpace.useful_binbounds(h_space, logarithmic=False)
correlation = ift.power_analyze(FFT.inverse_times(ift.log(data)),
binbounds=binbounds)
......
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