[HyperElastic] Enable the SelfAdjointEigenSolver in Ogden - with bool parameter

th-skam · th-skam · commit 4cba58b85d04 · 2026-02-13T15:42:24.000+01:00
diff --git a/Sofa/Component/SolidMechanics/FEM/HyperElastic/src/sofa/component/solidmechanics/fem/hyperelastic/material/Ogden.h b/Sofa/Component/SolidMechanics/FEM/HyperElastic/src/sofa/component/solidmechanics/fem/hyperelastic/material/Ogden.h
@@ -81,7 +81,8 @@ class Ogden: public HyperelasticMaterial<DataTypes>
 
         // 17/11/2025: Disable /*Eigen::SelfAdjointEigenSolver<EigenMatrix>*/
         // due to incorrect eigenvector computation for 3x3 matrices.
-        Eigen::EigenSolver<Eigen::Matrix<Real, 3, 3> > EigenProblemSolver(CEigen, true);
+        //Eigen::EigenSolver<Eigen::Matrix<Real, 3, 3> > EigenProblemSolver(CEigen, true);
+        Eigen::SelfAdjointEigenSolver<Eigen::Matrix<Real, 3, 3> > EigenProblemSolver(CEigen, true);
         if (EigenProblemSolver.info() != Eigen::Success)
         {
             dmsg_warning("Ogden") << "EigenSolver iterations failed to converge";

Original file line number	Diff line number	Diff line change
`@@ -81,7 +81,8 @@ class Ogden: public HyperelasticMaterial<DataTypes>`
`81`	`81`
`82`	`82`	`// 17/11/2025: Disable /Eigen::SelfAdjointEigenSolver<EigenMatrix>/`
`83`	`83`	`// due to incorrect eigenvector computation for 3x3 matrices.`
`84`		`- Eigen::EigenSolver<Eigen::Matrix<Real, 3, 3> > EigenProblemSolver(CEigen, true);`
	`84`	`+ //Eigen::EigenSolver<Eigen::Matrix<Real, 3, 3> > EigenProblemSolver(CEigen, true);`
	`85`	`+ Eigen::SelfAdjointEigenSolver<Eigen::Matrix<Real, 3, 3> > EigenProblemSolver(CEigen, true);`
`85`	`86`	`if (EigenProblemSolver.info() != Eigen::Success)`
`86`	`87`	`{`
`87`	`88`	`dmsg_warning("Ogden") << "EigenSolver iterations failed to converge";`