# HG changeset patch
# User Tuomo Valkonen <tuomov@iki.fi>
# Date 1726799426 18000
# Node ID fe38c518e67e9fad3364e1c850eac50a443b2f29
# Parent  a20429780845aee6844084dc4e8900477bf1b1c6
Print ∑_{i=0}^N} 𝒢(x^{k+1}; x^k) as value

diff -r a20429780845 -r fe38c518e67e src/Denoise.jl
--- 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