[FEM.Elastic] introduce virtual method called before force derivatives (sofa-framework#6139)

alxbilger · hugtalbot · web-flow · commit 65ccbee4967c · 2026-06-05T13:19:16.000+02:00
introduce function before force deriv

Co-authored-by: Hugo &lt;hugo.talbot@sofa-framework.org&gt;
diff --git a/Sofa/Component/SolidMechanics/FEM/Elastic/src/sofa/component/solidmechanics/fem/elastic/FEMForceField.h b/Sofa/Component/SolidMechanics/FEM/Elastic/src/sofa/component/solidmechanics/fem/elastic/FEMForceField.h
@@ -106,6 +106,8 @@ class FEMForceField :
         sofa::type::vector<ElementForce>& df,
         const sofa::VecDeriv_t<DataTypes>& dx);
 
+    virtual void beforeElementForceDeriv(const sofa::core::MechanicalParams* mparams) {}
+
     virtual void computeElementsForcesDeriv(
         const sofa::simulation::Range<std::size_t>& range,
         const sofa::core::MechanicalParams* mparams,
diff --git a/Sofa/Component/SolidMechanics/FEM/Elastic/src/sofa/component/solidmechanics/fem/elastic/FEMForceField.inl b/Sofa/Component/SolidMechanics/FEM/Elastic/src/sofa/component/solidmechanics/fem/elastic/FEMForceField.inl
@@ -164,6 +164,8 @@ void FEMForceField<DataTypes, ElementType>::computeElementsForcesDeriv(
 {
     SCOPED_TIMER("ElementForcesDeriv");
 
+    this->beforeElementForceDeriv(mparams);
+
     const auto& elements = trait::FiniteElement::getElementSequence(*this->l_topology);
 
     sofa::simulation::forEachRange(getExecutionPolicy(d_computeForceDerivStrategy), *this->m_taskScheduler,

Original file line number	Diff line number	Diff line change
`@@ -164,6 +164,8 @@ void FEMForceField<DataTypes, ElementType>::computeElementsForcesDeriv(`
`164`	`164`	`{`
`165`	`165`	`SCOPED_TIMER("ElementForcesDeriv");`
`166`	`166`
	`167`	`+ this->beforeElementForceDeriv(mparams);`
	`168`	`+`
`167`	`169`	`const auto& elements = trait::FiniteElement::getElementSequence(*this->l_topology);`
`168`	`170`
`169`	`171`	`sofa::simulation::forEachRange(getExecutionPolicy(d_computeForceDerivStrategy), *this->m_taskScheduler,`