update set_residuals

MaxenceGollier · MaxenceGollier · commit 9e4e9ffbd25b · 2026-02-20T13:04:44.000-05:00
diff --git a/src/LMTR_alg.jl b/src/LMTR_alg.jl
@@ -506,6 +506,6 @@ function SolverCore.solve!(
   end
 
   set_solution!(stats, xk)
-  set_residuals!(stats, zero(T), sqrt_ξ1_νInv)
+  set_residuals!(stats, zero(T), min(sqrt_ξ1_νInv, norm_s_cauchydν))
   return stats
 end
diff --git a/src/LM_alg.jl b/src/LM_alg.jl
@@ -497,6 +497,6 @@ function SolverCore.solve!(
   end
 
   set_solution!(stats, xk)
-  set_residuals!(stats, zero(T), sqrt_ξ1_νInv)
+  set_residuals!(stats, zero(T), min(sqrt_ξ1_νInv, norm_s_cauchydν))
   return stats
 end
diff --git a/src/R2DH.jl b/src/R2DH.jl
@@ -483,6 +483,6 @@ function SolverCore.solve!(
   end
 
   set_solution!(stats, xk)
-  set_residuals!(stats, zero(eltype(xk)), sqrt_ξ_νInv)
+  set_residuals!(stats, zero(eltype(xk)), min(sqrt_ξ_νInv, norm_sdν))
   return stats
 end
diff --git a/src/R2N.jl b/src/R2N.jl
@@ -536,7 +536,7 @@ function SolverCore.solve!(
   end
 
   set_solution!(stats, xk)
-  set_residuals!(stats, zero(eltype(xk)), sqrt_ξ1_νInv)
+  set_residuals!(stats, zero(eltype(xk)), min(sqrt_ξ1_νInv, norm_s_cauchydν))
   return stats
 end
 
diff --git a/src/R2_alg.jl b/src/R2_alg.jl
@@ -544,6 +544,6 @@ function SolverCore.solve!(
   end
 
   set_solution!(stats, xk)
-  set_residuals!(stats, zero(eltype(xk)), sqrt_ξ_νInv)
+  set_residuals!(stats, zero(eltype(xk)), min(sqrt_ξ_νInv, norm_sdν))
   return stats
 end
diff --git a/src/TRDH_alg.jl b/src/TRDH_alg.jl
@@ -544,6 +544,6 @@ function SolverCore.solve!(
   end
 
   set_solution!(stats, xk)
-  set_residuals!(stats, zero(eltype(xk)), sqrt_ξ_νInv)
+  set_residuals!(stats, zero(eltype(xk)), min(sqrt_ξ_νInv, norm_sdν))
   return stats
 end
diff --git a/src/TR_alg.jl b/src/TR_alg.jl
@@ -548,5 +548,5 @@ function SolverCore.solve!(
   end
 
   set_solution!(stats, xk)
-  set_residuals!(stats, zero(eltype(xk)), sqrt_ξ1_νInv)
+  set_residuals!(stats, zero(eltype(xk)), min(sqrt_ξ1_νInv, norm_s_cauchydν))
 end
