diff --git a/ProcessLib/LIE/HydroMechanics/LocalAssembler/HydroMechanicsLocalAssemblerMatrix-impl.h b/ProcessLib/LIE/HydroMechanics/LocalAssembler/HydroMechanicsLocalAssemblerMatrix-impl.h
index 253c7ff3b6a5d50352d7695933ad90bd189b3b79..4cdda52c7644bd3dfd29aa8e170211e25fcf93b1 100644
--- a/ProcessLib/LIE/HydroMechanics/LocalAssembler/HydroMechanicsLocalAssemblerMatrix-impl.h
+++ b/ProcessLib/LIE/HydroMechanics/LocalAssembler/HydroMechanicsLocalAssemblerMatrix-impl.h
@@ -249,16 +249,14 @@ void HydroMechanicsLocalAssemblerMatrix<ShapeFunctionDisplacement,
         rhs_u.noalias() -= -H_u.transpose() * rho * gravity_vec * ip_w;
 
         //
-        // displacement equation, pressure part
-        //
-        Kup.noalias() += B.transpose() * alpha * identity2 * N_p * ip_w;
-
-        //
-        // pressure equation, pressure part.
+        // pressure equation, pressure part and displacement equation, pressure
+        // part
         //
         if (!_process_data.deactivate_matrix_in_flow)  // Only for hydraulically
                                                        // active matrix
         {
+            Kup.noalias() += B.transpose() * alpha * identity2 * N_p * ip_w;
+
             double const k_over_mu =
                 _process_data.intrinsic_permeability(t, x_position)[0] /
                 _process_data.fluid_viscosity(t, x_position)[0];