Commit 9076afb2 authored by Jakob Knollmueller's avatar Jakob Knollmueller
Browse files

fix critical filter

parent 8efbc327
Pipeline #22268 passed with stage
in 4 minutes and 47 seconds
......@@ -104,16 +104,12 @@ if __name__ == "__main__":
flat_power = ift.Field.full(p_space, 1e-8)
m0 = ift.power_synthesize(flat_power, real_signal=True)
def ps0(k):
return (1./(1.+k)**2)
t0 = ift.Field(p_space,
for i in range(500):
# FIXME: Do not use ps0 for the inference! This is cheating!
# My critical filter version breaks, if I put power_spectrum=t0.
S0 = ift.create_power_operator(h_space, power_spectrum=ps0)
S0 = ift.create_power_operator(h_space, power_spectrum=ift.exp(t0))
# Initialize non-linear Wiener Filter energy
ICI = ift.GradientNormController(verbose=False, name="ICI",
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