From 6191ec45842ec5b9f523894ff2553d89a9354e52 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Tue, 11 Jul 2017 10:32:50 +0200
Subject: [PATCH] [PL] HM: Use assignment, not += for matrices.
---
ProcessLib/HydroMechanics/HydroMechanicsFEM.h | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/ProcessLib/HydroMechanics/HydroMechanicsFEM.h b/ProcessLib/HydroMechanics/HydroMechanicsFEM.h
index c08978b8627..fb84183887b 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() -=
--
GitLab