Mercurial > repos > pointsource_algs / file diff

--- a/src/pdps.rs	Fri May 08 16:47:58 2026 -0500
+++ b/src/pdps.rs	Fri May 15 14:40:02 2026 -0500
@@ -221,13 +221,13 @@
         let ω = pdpsconfig.acceleration.accelerate(&mut τ, &mut σ, γ);
 
         // Do dual update
-        // y = y_prev + τb
-        y.axpy(τ, b, 1.0);
-        // y = y_prev - τ(A[(1+ω)μ^{k+1}]-b)
-        opA.gemv(&mut y, -τ * (1.0 + ω), &μ, 1.0);
-        // y = y_prev - τ(A[(1+ω)μ^{k+1} - ω μ^k]-b)
-        opA.gemv(&mut y, τ * ω, &μ_base, 1.0);
-        y = phistar.prox(τ, y);
+        // y = y_prev - σb
+        y.axpy(-σ, b, 1.0);
+        // y = y_prev + σ(A[(1+ω)μ^{k+1}]-b)
+        opA.gemv(&mut y, σ * (1.0 + ω), &μ, 1.0);
+        // y = y_prev + σ(A[(1+ω)μ^{k+1} - ω μ^k]-b)
+        opA.gemv(&mut y, -σ * ω, &μ_base, 1.0);
+        y = phistar.prox(σ, y);
 
         // Give statistics if requested
         let iter = state.iteration();