Fix scaling bug in CG (#45)

* Fixes scaling bug in CG from #38 Co-authored-by: Brendt Wohlberg <bwohlberg@users.noreply.github.com>
lanl · Oct 18, 2021 · c844535 · c844535
1 parent b8cba12
commit c844535
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/scico/admm.py b/scico/admm.py
@@ -185,7 +185,7 @@ def internal_init(self, admm):
         )
         if admm.f is not None:
             # hessian = A.T @ W @ A; W may be identity
-            lhs_op = lhs_op + 2.0 * admm.f.scale * admm.f.hessian
+            lhs_op = lhs_op + admm.f.hessian
 
         lhs_op.jit()
         self.lhs_op = lhs_op