Delete m_gradJ from FEMaterialPoint

It is not needed.
febiosoftware · Dec 12, 2024 · f9bade4 · f9bade4
1 parent 8940f9c
commit f9bade4
Show file tree

Hide file tree

Showing 4 changed files with 3 additions and 6 deletions.
diff --git a/FEBioFluid/FEBiphasicFSIDomain3D.cpp b/FEBioFluid/FEBiphasicFSIDomain3D.cpp
@@ -990,7 +990,6 @@ void FEBiphasicFSIDomain3D::UpdateElementStress(int iel, const FETimeInfo& tp)
 
         //calculate gradJ gradphif
         bt.m_gradJ = Fi.transpose()*(GradJ*alphaf + GradJp*(1-alphaf));
-        ept.m_gradJ = bt.m_gradJ;
 
         // FSI material point data
         ft.m_w = el.Evaluate(w, n);

diff --git a/FEBioFluid/FEMultiphasicFSIDomain3D.cpp b/FEBioFluid/FEMultiphasicFSIDomain3D.cpp
@@ -1596,7 +1596,6 @@ void FEMultiphasicFSIDomain3D::UpdateElementStress(int iel, const FETimeInfo& tp
 
         //calculate gradJ gradphif
         bt.m_gradJ = Fi.transpose()*(GradJ*alphaf + GradJp*(1-alphaf));
-        ept.m_gradJ = bt.m_gradJ;
 
         // FSI material point data
         ft.m_w = el.Evaluate(w, n);

diff --git a/FEBioMech/FEElasticMaterialPoint.cpp b/FEBioMech/FEElasticMaterialPoint.cpp
@@ -35,7 +35,7 @@ FEElasticMaterialPoint::FEElasticMaterialPoint(FEMaterialPointData* mp) : FEMate
 	m_F.unit();
 	m_J = 1;
 	m_s.zero();
-    m_v = m_a = m_gradJ = vec3d(0, 0, 0);
+    m_v = m_a = vec3d(0, 0, 0);
     m_buncoupled = false;
     m_Wt = m_Wp = 0;
     m_p = 0;
@@ -58,7 +58,7 @@ void FEElasticMaterialPoint::Init()
 
 	m_s.zero();
 
-    m_v = m_a = m_gradJ = vec3d(0, 0, 0);
+    m_v = m_a = vec3d(0, 0, 0);
     m_L.zero();
 
     m_Wt = m_Wp = 0;
@@ -73,7 +73,7 @@ void FEElasticMaterialPoint::Init()
 void FEElasticMaterialPoint::Serialize(DumpStream& ar)
 {
 	FEMaterialPointData::Serialize(ar);
-    ar & m_F & m_J & m_s & m_v & m_a & m_gradJ & m_L & m_Wt & m_Wp & m_p;
+    ar & m_F & m_J & m_s & m_v & m_a & m_L & m_Wt & m_Wp & m_p;
 	ar & m_buncoupled;
 }
 

diff --git a/FEBioMech/FEElasticMaterialPoint.h b/FEBioMech/FEElasticMaterialPoint.h
@@ -84,7 +84,6 @@ class FEBIOMECH_API FEElasticMaterialPoint : public FEMaterialPointData
 	// deformation data at intermediate time
 	mat3d	m_F;	//!< deformation gradient
 	double	m_J;	//!< determinant of F
-    vec3d   m_gradJ;  //!< gradient of J
     vec3d   m_v;    //!< velocity
     vec3d   m_a;    //!< acceleration
     mat3d   m_L;    //!< spatial velocity gradient