Skip to content
Snippets Groups Projects
Commit 1fc0eeed authored by Martin Reinecke's avatar Martin Reinecke
Browse files

Merge branch 'nifty2go_fix_CF' into 'nifty2go'

fix critical filter

See merge request ift/NIFTy!205
parents 8efbc327 9076afb2
No related branches found
No related tags found
1 merge request!205fix critical filter
Pipeline #
...@@ -104,16 +104,12 @@ if __name__ == "__main__": ...@@ -104,16 +104,12 @@ if __name__ == "__main__":
flat_power = ift.Field.full(p_space, 1e-8) flat_power = ift.Field.full(p_space, 1e-8)
m0 = ift.power_synthesize(flat_power, real_signal=True) m0 = ift.power_synthesize(flat_power, real_signal=True)
def ps0(k):
return (1./(1.+k)**2)
t0 = ift.Field(p_space, t0 = ift.Field(p_space,
val=ift.dobj.from_global_data(np.log(1./(1+p_space.k_lengths)**2))) val=ift.dobj.from_global_data(-7.))
for i in range(500): 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=ift.exp(t0))
S0 = ift.create_power_operator(h_space, power_spectrum=ps0)
# Initialize non-linear Wiener Filter energy # Initialize non-linear Wiener Filter energy
ICI = ift.GradientNormController(verbose=False, name="ICI", ICI = ift.GradientNormController(verbose=False, name="ICI",
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment