Tiny bug on restart for IRN_TV_CGLS

CERN · Dec 6, 2022 · b89fa59 · b89fa59
1 parent 4bbcf34
commit b89fa59
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 2 deletions.
diff --git a/MATLAB/Algorithms/IRN_TV_CGLS.m b/MATLAB/Algorithms/IRN_TV_CGLS.m
@@ -65,7 +65,7 @@
 end
 qualMeasOut=zeros(length(QualMeasOpts),niter);
 resL2 = zeros(1,niter);
-
+remember=0;
 iter=0;
 % for iii = 1:niter_outer
 while iter<niter

diff --git a/Python/demos/d09_Algorithms04.py b/Python/demos/d09_Algorithms04.py
@@ -189,9 +189,19 @@
     delta=np.array([-0.005]),
 )
 
+
+
+#  IRN-CGLS-TV: CGLS with TV regularization
+# ==========================================================================
+# ==========================================================================
+#
+# Should converge fairly fast to a TV regularized solution
+imgCGLSTV=algs.irn_tv_cgls(noise_projections,geo,angles,10,lmbda=5,niter_outer=2)
+
+
 #%% plot results
 
 # plot images
 tigre.plotimg(
-    np.concatenate([imgAWASDPOCS, imgOSASDPOCS, imgASDPOCS, imgOSSART], axis=1), dim="z", step=2
+    np.concatenate([imgAWASDPOCS, imgOSASDPOCS, imgASDPOCS, imgOSSART, imgCGLSTV], axis=1), dim="z", step=2
 )