From aa1f7bdfb2fec8a9d98e3d50d60802a4bd32a5df Mon Sep 17 00:00:00 2001
From: Klaus Mosthaf <klmos@env.dtu.dk>
Date: Fri, 23 Dec 2011 09:38:38 +0000
Subject: [PATCH] Evaluate Neumann conditions only for the equations, which are
of Neumann type
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@7240 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
dumux/boxmodels/common/boxlocalresidual.hh | 8 +++++++-
1 file changed, 7 insertions(+), 1 deletion(-)
diff --git a/dumux/boxmodels/common/boxlocalresidual.hh b/dumux/boxmodels/common/boxlocalresidual.hh
index d64c6cb9de..fb48b16e5a 100644
--- a/dumux/boxmodels/common/boxlocalresidual.hh
+++ b/dumux/boxmodels/common/boxlocalresidual.hh
@@ -442,7 +442,13 @@ protected:
fvElemGeom_().boundaryFace[boundaryFaceIdx].area
* curVolVars_(scvIdx).extrusionFactor();
Valgrind::CheckDefined(values);
- residual_[scvIdx] += values;
+
+ // set the neumann conditions
+ for (int eqIdx = 0; eqIdx < numEq; ++eqIdx) {
+ if (!bcTypes.isNeumann(eqIdx))
+ continue;
+ residual_[scvIdx] += values;
+ }
}
}
--
GitLab