Merge branch 'fix_LIE' into 'master'

Fixed a bug in LIE

Closes #3173

See merge request !3834

ProcessLib/LIE/HydroMechanics/LocalAssembler/HydroMechanicsLocalAssemblerFracture-impl.h

@@ -196,8 +196,8 @@ void HydroMechanicsLocalAssemblerFracture<ShapeFunctionDisplacement,
     Eigen::MatrixXd const global2local_rotation =
         R.template topLeftCorner<ShapeFunctionPressure::DIM, GlobalDim>();
 
-    auto const& gravity_vec =
-        (global2local_rotation * _process_data.specific_body_force).eval();
+    DimVectorType const gravity_vec =
+        global2local_rotation * _process_data.specific_body_force;
 
     ParameterLib::SpatialPosition x_position;
     x_position.setElementID(_element.getID());
@@ -294,7 +294,7 @@ void HydroMechanicsLocalAssemblerFracture<ShapeFunctionDisplacement,
         Eigen::Matrix<double, 1, displacement_size> const mT_R_Hg =
             identity2.transpose() * R * H_g;
         // velocity in global coordinates
-        ip_data.darcy_velocity.head(GlobalDim).noalias() =
+        ip_data.darcy_velocity =
             -global2local_rotation.transpose() * k * grad_head_over_mu;
         J_pg.noalias() += N_p.transpose() * S * N_p * p_dot * mT_R_Hg * ip_w;
         J_pg.noalias() +=

ProcessLib/LIE/HydroMechanics/LocalAssembler/HydroMechanicsLocalAssemblerFracture.h

@@ -80,14 +80,16 @@ private:
 
     // Types for displacement.
     using ShapeMatricesTypeDisplacement =
-        ShapeMatrixPolicyType<ShapeFunctionDisplacement, GlobalDim>;
+        ShapeMatrixPolicyType<ShapeFunctionDisplacement,
+                              ShapeFunctionDisplacement::DIM>;
     using HMatricesType =
         HMatrixPolicyType<ShapeFunctionDisplacement, GlobalDim>;
     using HMatrixType = typename HMatricesType::HMatrixType;
 
     // Types for pressure.
     using ShapeMatricesTypePressure =
-        ShapeMatrixPolicyType<ShapeFunctionPressure, GlobalDim>;
+        ShapeMatrixPolicyType<ShapeFunctionPressure,
+                              ShapeFunctionPressure::DIM>;
 
     // Types for the integration point data
     using IntegrationPointDataType =
@@ -97,6 +99,7 @@ private:
                                      GlobalDim>;
 
     using GlobalDimVector = Eigen::Matrix<double, GlobalDim, 1>;
+    using DimVectorType = typename ShapeMatricesTypePressure::DimVectorType;
 
     HydroMechanicsProcessData<GlobalDim>& _process_data;
 

ProcessLib/LIE/HydroMechanics/Tests.cmake

@@ -2,7 +2,7 @@
 if (NOT (OGS_USE_LIS OR OGS_USE_MPI))
     # Temporary disable test until resolution of the issue
     # https://gitlab.opengeosys.org/ogs/ogs/-/issues/3173
-    # OgsTest(PROJECTFILE LIE/HydroMechanics/single_fracture_constK.prj RUNTIME 17)
+    OgsTest(PROJECTFILE LIE/HydroMechanics/single_fracture_constK.prj RUNTIME 17)
     OgsTest(PROJECTFILE LIE/HydroMechanics/single_fracture_CZ_kf_const.prj RUNTIME 250)
     OgsTest(PROJECTFILE LIE/HydroMechanics/single_fracture_CZ_kf_cubic.prj RUNTIME 120)
     OgsTest(PROJECTFILE LIE/HydroMechanics/single_fracture_LE_kf_const.prj RUNTIME 90)