From 3e3f3aa60d55d1b54def9173bb0d73651cd3cf3d Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Dennis=20Gl=C3=A4ser?= <dennis.glaeser@iws.uni-stuttgart.de>
Date: Mon, 19 Nov 2018 16:36:17 +0100
Subject: [PATCH] [mpfa][ivlocalassembler][gravity] remove obsolete reset to
zero
---
dumux/discretization/cellcentered/mpfa/localassembler.hh | 8 ++------
1 file changed, 2 insertions(+), 6 deletions(-)
diff --git a/dumux/discretization/cellcentered/mpfa/localassembler.hh b/dumux/discretization/cellcentered/mpfa/localassembler.hh
index c3247c6b44..8c644aecb5 100644
--- a/dumux/discretization/cellcentered/mpfa/localassembler.hh
+++ b/dumux/discretization/cellcentered/mpfa/localassembler.hh
@@ -182,7 +182,6 @@ class InteractionVolumeAssemblerBase
FaceVector sum_alphas;
resizeVector_(sum_alphas, iv.numUnknowns());
sum_alphas = 0.0;
- std::fill(g.begin(), g.end(), 0.0);
for (LocalIndexType faceIdx = 0; faceIdx < iv.numFaces(); ++faceIdx)
{
// gravitational acceleration on this face
@@ -207,10 +206,7 @@ class InteractionVolumeAssemblerBase
Scalar rho;
if (isSurfaceGrid)
- {
resizeVector_(outsideG[faceIdx], numOutsideFaces);
- std::fill(outsideG[faceIdx].begin(), outsideG[faceIdx].end(), 0.0);
- }
if (!curLocalScvf.isDirichlet())
{
@@ -249,13 +245,13 @@ class InteractionVolumeAssemblerBase
rho = getRho(elemVolVars()[curGlobalScvf.outsideScvIdx()]);
// add "inside" & "outside" alphas to gravity containers
- g[faceIdx] += alpha_inside*rho*curGlobalScvf.area();
+ g[faceIdx] = alpha_inside*rho*curGlobalScvf.area();
if (isSurfaceGrid)
{
unsigned int i = 0;
for (const auto& alpha : alpha_outside)
- outsideG[faceIdx][i++] += alpha*rho*curGlobalScvf.area();
+ outsideG[faceIdx][i++] = alpha*rho*curGlobalScvf.area();
}
}
--
GitLab