Mercurial > repos > ImageTools / changeset

--- a/src/Denoise.jl	Thu Sep 19 21:29:59 2024 -0500
+++ b/src/Denoise.jl	Thu Sep 19 21:30:26 2024 -0500
@@ -83,6 +83,9 @@
     # Run the algorithm
     ####################
 
+    q = norm₂²(b-x)/2
+    cum = 0.0
+
     v = iterate(params) do verbose :: Function
         ω = params.accel ? 1/√(1+2*γ*τ) : 1
         
@@ -97,14 +100,18 @@
         if params.accel
             τ, σ = τ*ω, σ/ω
         end
-                
+
+        qnew = norm₂²(b-x)/2
+        cum = qnew - q + dot(∇y, x) - dot(∇x, y)
+        q = qnew
+        
         ################################
         # Give function value if needed
         ################################
         v = verbose() do            
-            ∇₂!(Δy, x)
-            value = norm₂²(b-x)/2 + params.α*γnorm₂₁(Δy, params.ρ)
-            value, x
+            #∇₂!(Δy, x)
+            #value = norm₂²(b-x)/2 + params.α*γnorm₂₁(Δy, params.ρ)
+            cum, x
         end
 
         v