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

cosmetics

parent c3bcd006
No related branches found
No related tags found
1 merge request!207new operator convenience functionality
Pipeline #
Hide whitespace changes
Inline Side-by-side
demos/wiener_filter_easy.py
+ 2
− 1
View file @ ff2eaec7
...@@ -57,7 +57,8 @@ if __name__ == "__main__": ...@@ -57,7 +57,8 @@ if __name__ == "__main__":
IC = ift.GradientNormController(verbose=True, iteration_limit=500, IC = ift.GradientNormController(verbose=True, iteration_limit=500,
tol_abs_gradnorm=0.1) tol_abs_gradnorm=0.1)
inverter = ift.ConjugateGradient(controller=IC) inverter = ift.ConjugateGradient(controller=IC)
D = (R.adjoint*N.inverse*R + fft.adjoint*Sh.inverse*fft).inverse S_inv = fft.adjoint*Sh.inverse*fft
D = (R.adjoint*N.inverse*R + S_inv).inverse
# MR FIXME: we can/should provide a preconditioner here as well! # MR FIXME: we can/should provide a preconditioner here as well!
D = ift.InversionEnabler(D, inverter) D = ift.InversionEnabler(D, inverter)
m = D(j) m = D(j)
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment