diff --git a/ProcessLib/HydroMechanics/HydroMechanicsFEM.h b/ProcessLib/HydroMechanics/HydroMechanicsFEM.h index c08978b8627c7b30845b20da3a11dd312b1f1cce..fb84183887b128d3dab10a69c2024a378ceb293f 100644 --- a/ProcessLib/HydroMechanics/HydroMechanicsFEM.h +++ b/ProcessLib/HydroMechanics/HydroMechanicsFEM.h @@ -379,19 +379,19 @@ public: local_Jac .template block<displacement_size, pressure_size>( displacement_index, pressure_index) - .noalias() -= Kup; + .noalias() = -Kup; // pressure equation, pressure part. local_Jac .template block<pressure_size, pressure_size>(pressure_index, pressure_index) - .noalias() += laplace_p + storage_p / dt; + .noalias() = laplace_p + storage_p / dt; // pressure equation, displacement part. local_Jac .template block<pressure_size, displacement_size>( pressure_index, displacement_index) - .noalias() += Kup.transpose() / dt; + .noalias() = Kup.transpose() / dt; // pressure equation local_rhs.template segment<pressure_size>(pressure_index).noalias() -=