Adds curvature stopping criterion to NR solver

[jmartfrut]

Implements a curvature stopping criterion in the Newton-Raphson solver. This criterion checks if the approximate solution has converged based on the change in the residual. This prevents unnecessary iterations.

[github-actions]

Benchmark Results (Julia v1)

Time benchmarks

	main	fac0aee...	main / fac0aee...
Constitutive models/Visco-elastic Ψ	0.0496 ± 0.027 ms	0.0489 ± 0.027 ms	1.01 ± 0.78
Constitutive models/Visco-elastic ∂Ψu	0.0651 ± 0.0057 ms	0.0648 ± 0.013 ms	1 ± 0.22
Constitutive models/Visco-elastic ∂Ψuu	0.114 ± 0.015 ms	0.114 ± 0.014 ms	1 ± 0.18
Simulations/StaticMechanicalDirichlet	0.134 ± 0.0073 s	0.133 ± 0.0067 s	1.01 ± 0.075
Simulations/StaticMechanicalNeumann	0.129 ± 0.012 s	0.13 ± 0.013 s	0.995 ± 0.13
Simulations/ViscoElastic	20.7 s	20.8 s	0.995
Tensor algebra/Cofactor	0.07 ± 0.01 μs	0.07 ± 0.001 μs	1 ± 0.14
Tensor algebra/Det(A)Inv(A')	0.13 ± 0.001 μs	0.13 ± 0.01 μs	1 ± 0.077
Tensor algebra/δδ_λ_2d	30 ± 0 ns	30 ± 0 ns	1 ± 0
Tensor algebra/δδ_μ_2d	30 ± 0 ns	30 ± 0 ns	1 ± 0
time_to_load	2.36 ± 0.019 s	2.36 ± 0.024 s	0.999 ± 0.013

Memory benchmarks

	main	fac0aee...	main / fac0aee...
Constitutive models/Visco-elastic Ψ	1.24 k allocs: 0.0777 MB	1.24 k allocs: 0.0777 MB	1
Constitutive models/Visco-elastic ∂Ψu	1.36 k allocs: 0.0898 MB	1.36 k allocs: 0.0898 MB	1
Constitutive models/Visco-elastic ∂Ψuu	2.32 k allocs: 0.155 MB	2.32 k allocs: 0.155 MB	1
Simulations/StaticMechanicalDirichlet	1.51 M allocs: 0.109 GB	1.51 M allocs: 0.109 GB	1
Simulations/StaticMechanicalNeumann	1.48 M allocs: 0.0933 GB	1.48 M allocs: 0.0933 GB	1
Simulations/ViscoElastic	0.343 G allocs: 22 GB	0.343 G allocs: 22 GB	1
Tensor algebra/Cofactor	1 allocs: 0.0781 kB	1 allocs: 0.0781 kB	1
Tensor algebra/Det(A)Inv(A')	4 allocs: 0.25 kB	4 allocs: 0.25 kB	1
Tensor algebra/δδ_λ_2d	0 allocs: 0 B	0 allocs: 0 B
Tensor algebra/δδ_μ_2d	0 allocs: 0 B	0 allocs: 0 B
time_to_load	0.149 k allocs: 11.1 kB	0.149 k allocs: 11.1 kB	1