From a9ac6d246f819d42ab7f2012752fdbce6c034246 Mon Sep 17 00:00:00 2001
From: Kilian Weishaupt <kilian.weishaupt@iws.uni-stuttgart.de>
Date: Mon, 13 Nov 2017 09:56:00 +0100
Subject: [PATCH] [staggered][localAssembler] Bugfix in calculation of residual
* take the full residual vector, not just its first component
---
dumux/assembly/staggeredlocalassembler.hh | 8 ++------
1 file changed, 2 insertions(+), 6 deletions(-)
diff --git a/dumux/assembly/staggeredlocalassembler.hh b/dumux/assembly/staggeredlocalassembler.hh
index 4cfec4176d..a60f8dcdf0 100644
--- a/dumux/assembly/staggeredlocalassembler.hh
+++ b/dumux/assembly/staggeredlocalassembler.hh
@@ -139,11 +139,10 @@ public:
prevElemeFaceVars,
curElemeFaceVars,
elemBcTypes,
- elemFluxVarsCache)[0];
+ elemFluxVarsCache);
// treat the local residua of the face dofs:
// create a cache to reuse some results for the calculation of the derivatives
-
FaceSolutionVector faceResidualCache;
faceResidualCache.resize(fvGeometry.numScvf());
faceResidualCache = 0.0;
@@ -159,7 +158,7 @@ public:
prevElemeFaceVars,
curElemeFaceVars,
elemBcTypes,
- elemFluxVarsCache)[0];
+ elemFluxVarsCache);
res[faceIdx][scvf.dofIndex()] += faceResidualCache[scvf.localFaceIdx()] ;
}
@@ -179,9 +178,6 @@ public:
res[cellCenterIdx][cellCenterGlobalI],
faceResidualCache);
-
-
-
}
/*!
--
GitLab