solve variable scope

Author: MaxenceGollier

src/TRDH_alg.jl

@@ -405,7 +405,8 @@ function SolverCore.solve!(
   norm_sdν = norm_s / ν
   atol_step += rtol_step * norm_sdν # make stopping test absolute and relative
 
-  solved = (ξ1 < 0 && sqrt_ξ_νInv ≤ neg_tol) || (ξ1 ≥ 0 && sqrt_ξ_νInv ≤ atol_decr) || (norm_sdν ≤ atol_step)
+  solved = reduce_TR ? (ξ1 < 0 && sqrt_ξ_νInv ≤ neg_tol) || (ξ1 ≥ 0 && sqrt_ξ_νInv ≤ atol_decr) || (norm_sdν ≤ atol_step) :
+    (ξ < 0 && sqrt_ξ_νInv ≤ neg_tol) || (ξ ≥ 0 && sqrt_ξ_νInv ≤ atol_decr) || (norm_sdν ≤ atol_step)
   set_status!(
     stats,
     get_status(