From ff2eaec73371fd8b4dd83fa596a25c8ecda3e1e9 Mon Sep 17 00:00:00 2001
From: Martin Reinecke <martin@mpa-garching.mpg.de>
Date: Thu, 28 Dec 2017 18:52:58 +0100
Subject: [PATCH] cosmetics
---
demos/wiener_filter_easy.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/demos/wiener_filter_easy.py b/demos/wiener_filter_easy.py
index a0a30eb24..9ed7be0bf 100644
--- a/demos/wiener_filter_easy.py
+++ b/demos/wiener_filter_easy.py
@@ -57,7 +57,8 @@ if __name__ == "__main__":
IC = ift.GradientNormController(verbose=True, iteration_limit=500,
tol_abs_gradnorm=0.1)
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!
D = ift.InversionEnabler(D, inverter)
m = D(j)
--
GitLab