change NR settings

Hugo Verhelst · Hugo Verhelst · commit e26bfa2d6322 · 2025-01-15T11:53:09.000+01:00
diff --git a/src/gsMaterialMatrixTFT.hpp b/src/gsMaterialMatrixTFT.hpp
@@ -486,7 +486,7 @@ gsMatrix<T> gsMaterialMatrixTFT<dim,T,linear>::eval_theta(const gsMatrix<T> & Cs
                 converged = obj.findRootBisection(f,theta_interval(0),theta_interval(1),theta,1e-6);
                 if (!converged)
                 {
-                    obj.newtonRaphson(zeros,arg,false,1e-12,100,1);
+                    obj.newtonRaphson(zeros,arg,false,1e-8,1000,1);
                     theta = arg(0);
                     f = obj.eval(theta);
                 }
@@ -522,7 +522,7 @@ gsMatrix<T> gsMaterialMatrixTFT<dim,T,linear>::eval_theta(const gsMatrix<T> & Cs
                     converged = obj.findRootBisection(f,intervals[i].first,intervals[i].second,theta,1e-6);
                     if (!converged)
                     {
-                        obj.newtonRaphson(zeros,arg,false,1e-12,100,1);
+                        obj.newtonRaphson(zeros,arg,false,1e-8,1000,1);
                         theta = arg(0);
                         f = obj.eval(theta);
                     }

Original file line number	Diff line number	Diff line change
`@@ -486,7 +486,7 @@ gsMatrix<T> gsMaterialMatrixTFT<dim,T,linear>::eval_theta(const gsMatrix<T> & Cs`
`486`	`486`	`converged = obj.findRootBisection(f,theta_interval(0),theta_interval(1),theta,1e-6);`
`487`	`487`	`if (!converged)`
`488`	`488`	`{`
`489`		`- obj.newtonRaphson(zeros,arg,false,1e-12,100,1);`
	`489`	`+ obj.newtonRaphson(zeros,arg,false,1e-8,1000,1);`
`490`	`490`	`theta = arg(0);`
`491`	`491`	`f = obj.eval(theta);`
`492`	`492`	`}`
`@@ -522,7 +522,7 @@ gsMatrix<T> gsMaterialMatrixTFT<dim,T,linear>::eval_theta(const gsMatrix<T> & Cs`
`522`	`522`	`converged = obj.findRootBisection(f,intervals[i].first,intervals[i].second,theta,1e-6);`
`523`	`523`	`if (!converged)`
`524`	`524`	`{`
`525`		`- obj.newtonRaphson(zeros,arg,false,1e-12,100,1);`
	`525`	`+ obj.newtonRaphson(zeros,arg,false,1e-8,1000,1);`
`526`	`526`	`theta = arg(0);`
`527`	`527`	`f = obj.eval(theta);`
`528`	`528`	`}`