Commit 1f83d30d authored by Philipp Arras's avatar Philipp Arras

More fixups

parent 6a8f8cc5
Pipeline #30845 passed with stages
in 1 minute and 23 seconds
......@@ -40,7 +40,8 @@ IC = ift.GradientNormController(tol_abs_gradnorm=1e-3, iteration_limit=N_iter)
m, samps = ift.library.generate_krylov_samples(D_inv, S, j, N_samps, IC)
m_x = sky(m)
inverter = ift.ConjugateGradient(IC)
curv = ift.library.WienerFilterCurvature(S=S, N=N, R=R_p, inverter=inverter)
curv = ift.library.WienerFilterCurvature(S=S, N=N, R=R_p, inverter=inverter,
sampling_inverter=inverter)
samps_old = [curv.draw_sample(from_inverse=True) for i in range(N_samps)]
plt.plot(d.to_global_data(), '+', label="data", alpha=.5)
......
......@@ -76,7 +76,7 @@ if __name__ == "__main__":
ctrl = ift.GradientNormController(name="inverter", tol_abs_gradnorm=0.1)
inverter = ift.ConjugateGradient(controller=ctrl)
wiener_curvature = ift.library.WienerFilterCurvature(
S=S, N=N, R=R, inverter=inverter)
S=S, N=N, R=R, inverter=inverter, sampling_inverter=inverter)
m_k = wiener_curvature.inverse_times(j)
m = ht(m_k)
......
......@@ -50,7 +50,7 @@ if __name__ == "__main__":
ctrl = ift.GradientNormController(name="inverter", tol_abs_gradnorm=1e-2)
inverter = ift.ConjugateGradient(controller=ctrl)
wiener_curvature = ift.library.WienerFilterCurvature(
S=S, N=N, R=R, inverter=inverter)
S=S, N=N, R=R, inverter=inverter, sampling_inverter=inverter)
m_k = wiener_curvature.inverse_times(j)
m = ht(m_k)
......
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