Merge pull request #125 from MultiSimOLab/curvature_tolerance

Configures Newton-Raphson curvature tolerance

Author: jmartfrut

README.md

@@ -154,7 +154,7 @@ src="https://github.com/MultiSimOLab/HyperFEM/blob/main/docs/imgs/sims_.png?raw=
 
 In order to give credit to the HyperFEM contributors, we ask that you please reference the paper:
 
-C. Perez‐Garcia, R. Ortigosa, J. Martínez‐Frutos, and D. Garcia‐Gonzalez, **Topology and material optimization in ultra-soft magnetoactive structures: making advantage of residual anisotropies.** Adv. Mater. (2025): e18489. https://https://doi.org/10.1002/adma.202518489
+C. Perez‐Garcia, R. Ortigosa, J. Martínez‐Frutos, and D. Garcia‐Gonzalez, **Topology and material optimization in ultra-soft magnetoactive structures: making advantage of residual anisotropies.** Adv. Mater. (2025): e18489. https://doi.org/10.1002/adma.202518489
 
 
 along with the required citations for [Gridap](https://github.com/gridap/Gridap.jl).

src/Solvers/NonlinearSolvers.jl

@@ -11,12 +11,13 @@ struct Newton_RaphsonSolver <: Algebra.NonlinearSolver
   ls::Algebra.LinearSolver
   log::ConvergenceLog{Float64}
   linesearch::AbstractLineSearch
+  ctol::Float64
 end
 
-function Newton_RaphsonSolver(ls; maxiter=100, atol=1e-12, rtol=1.e-6,verbose=0, name="Newton-Raphson", linesearch::AbstractLineSearch=LineSearch())
+function Newton_RaphsonSolver(ls; maxiter=100, atol=1e-12, rtol=1.e-6, ctol=1.e-5, verbose=0, name="Newton-Raphson", linesearch::AbstractLineSearch=LineSearch())
   tols = SolverTolerances{Float64}(; maxiter=maxiter, atol=atol, rtol=rtol)
   log = ConvergenceLog(name, tols; verbose=verbose)
-  return Newton_RaphsonSolver(ls, log, linesearch)
+  return Newton_RaphsonSolver(ls, log, linesearch,ctol)
 end
 
 AbstractTrees.children(s::Newton_RaphsonSolver) = [s.ls]
@@ -64,7 +65,7 @@ function _solve_nr!(x, A, b, dx, ns, nls, op)
     solve!(dx, ns, b)
 
     # curvature stopping criterion
-    if abs(b' * dx) < 1e-5
+    if abs(b' * dx) < nls.ctol
       res = 2.22e-22
       done = update!(log, res)
       break