removed CalculateMassMatrix

KratosMultiphysics · Apr 19, 2024 · ea4cc17 · ea4cc17
1 parent 50f3582
commit ea4cc17
Show file tree

Hide file tree

Showing 4 changed files with 2 additions and 37 deletions.
diff --git a/applications/GeoMechanicsApplication/custom_elements/transient_Pw_element.cpp b/applications/GeoMechanicsApplication/custom_elements/transient_Pw_element.cpp
@@ -52,22 +52,6 @@ void TransientPwElement<TDim, TNumNodes>::EquationIdVector(EquationIdVectorType&
     rResult = Geo::DofUtilities::ExtractEquationIdsFrom(GetDofs());
 }
 
-//----------------------------------------------------------------------------------------
-template <unsigned int TDim, unsigned int TNumNodes>
-void TransientPwElement<TDim, TNumNodes>::CalculateMassMatrix(MatrixType& rMassMatrix,
-                                                              const ProcessInfo& rCurrentProcessInfo)
-{
-    KRATOS_TRY
-
-    const unsigned int N_DOF = this->GetNumberOfDOF();
-
-    // Resizing mass matrix
-    if (rMassMatrix.size1() != N_DOF) rMassMatrix.resize(N_DOF, N_DOF, false);
-    noalias(rMassMatrix) = ZeroMatrix(N_DOF, N_DOF);
-
-    KRATOS_CATCH("")
-}
-
 //----------------------------------------------------------------------------------------
 template <unsigned int TDim, unsigned int TNumNodes>
 void TransientPwElement<TDim, TNumNodes>::CalculateDampingMatrix(MatrixType& rDampingMatrix,

diff --git a/applications/GeoMechanicsApplication/custom_elements/transient_Pw_element.hpp b/applications/GeoMechanicsApplication/custom_elements/transient_Pw_element.hpp
@@ -103,9 +103,7 @@ class KRATOS_API(GEO_MECHANICS_APPLICATION) TransientPwElement : public UPwSmall
     void GetDofList(DofsVectorType& rElementalDofList, const ProcessInfo&) const override;
 
     void EquationIdVector(EquationIdVectorType& rResult, const ProcessInfo&) const override;
-
-    void CalculateMassMatrix(MatrixType& rMassMatrix, const ProcessInfo& rCurrentProcessInfo) override;
-
+
     void CalculateDampingMatrix(MatrixType& rDampingMatrix, const ProcessInfo& rCurrentProcessInfo) override;
 
     void GetValuesVector(Vector& rValues, int Step = 0) const override;

diff --git a/applications/GeoMechanicsApplication/custom_elements/transient_Pw_interface_element.cpp b/applications/GeoMechanicsApplication/custom_elements/transient_Pw_interface_element.cpp
@@ -127,21 +127,6 @@ void TransientPwInterfaceElement<TDim, TNumNodes>::Initialize(const ProcessInfo&
     KRATOS_CATCH("")
 }
 
-template <unsigned int TDim, unsigned int TNumNodes>
-void TransientPwInterfaceElement<TDim, TNumNodes>::CalculateMassMatrix(MatrixType& rMassMatrix,
-                                                                       const ProcessInfo& rCurrentProcessInfo)
-{
-    KRATOS_TRY
-
-    const unsigned int N_DOF = this->GetNumberOfDOF();
-
-    // Resizing mass matrix
-    if (rMassMatrix.size1() != N_DOF) rMassMatrix.resize(N_DOF, N_DOF, false);
-    noalias(rMassMatrix) = ZeroMatrix(N_DOF, N_DOF);
-
-    KRATOS_CATCH("")
-}
-
 template <unsigned int TDim, unsigned int TNumNodes>
 void TransientPwInterfaceElement<TDim, TNumNodes>::InitializeSolutionStep(const ProcessInfo& rCurrentProcessInfo)
 {

diff --git a/applications/GeoMechanicsApplication/custom_elements/transient_Pw_interface_element.hpp b/applications/GeoMechanicsApplication/custom_elements/transient_Pw_interface_element.hpp
@@ -106,9 +106,7 @@ class KRATOS_API(GEO_MECHANICS_APPLICATION) TransientPwInterfaceElement
     void GetSecondDerivativesVector(Vector& rValues, int Step = 0) const override;
 
     void Initialize(const ProcessInfo& rCurrentProcessInfo) override;
-
-    void CalculateMassMatrix(MatrixType& rMassMatrix, const ProcessInfo& rCurrentProcessInfo) override;
-
+
     void InitializeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
     void FinalizeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;