Merge branch 'PhaseField' into 'master'

A bug fix in PhaseField parallel example

Closes #3224

See merge request ogs/ogs!3898

ProcessLib/BoundaryConditionAndSourceTerm/PhaseFieldIrreversibleDamageOracleBoundaryCondition.cpp

@@ -57,7 +57,7 @@ void PhaseFieldIrreversibleDamageOracleBoundaryCondition::preTimestep(
         const auto g_idx =
             _dof_table.getGlobalIndex(l, _variable_id, _component_id);
 
-        if ((*x[process_id])[node_id] <= irreversibleDamage)
+        if ((*x[process_id])[g_idx] <= irreversibleDamage)
         {
             _bc_values.ids.emplace_back(g_idx);
             _bc_values.values.emplace_back(0.0);

ProcessLib/PhaseField/Tests.cmake

@@ -42,3 +42,20 @@ AddTest(
         expected_AT1_vd_tension_ts_10_t_1_000000_0.vtu AT1_vd_tension_ts_10_t_1_000000_0.vtu displacement displacement 1e-5 0
         expected_AT1_vd_tension_ts_10_t_1_000000_0.vtu AT1_vd_tension_ts_10_t_1_000000_0.vtu phasefield phasefield 1e-6 0
 )
+
+AddTest(
+    NAME PhaseField_3D_beam_tens_AT1_vd_2core
+    PATH PhaseField/beam
+    EXECUTABLE ogs
+    EXECUTABLE_ARGS AT1_vd_tensile.prj
+    WRAPPER mpirun
+    WRAPPER_ARGS -np 2
+    TESTER vtkdiff
+    REQUIREMENTS OGS_USE_MPI
+    RUNTIME 18
+    DIFF_DATA
+        expected_AT1_vd_tension_2core_ts_10_t_1_000000_0.vtu AT1_vd_tension_ts_10_t_1_000000_0.vtu displacement displacement 1e-5 0
+        expected_AT1_vd_tension_2core_ts_10_t_1_000000_1.vtu AT1_vd_tension_ts_10_t_1_000000_1.vtu displacement displacement 1e-5 0
+        expected_AT1_vd_tension_2core_ts_10_t_1_000000_0.vtu AT1_vd_tension_ts_10_t_1_000000_0.vtu phasefield phasefield 1e-6 0
+        expected_AT1_vd_tension_2core_ts_10_t_1_000000_1.vtu AT1_vd_tension_ts_10_t_1_000000_1.vtu phasefield phasefield 1e-6 0
+)