From 09d3f6b1bbba3490facd7a103b583da45300c2cc Mon Sep 17 00:00:00 2001
From: Natalie Schroeder <natalies.schroeder@iws.uni-stuttgart.de>
Date: Thu, 15 Jan 2015 13:04:09 +0000
Subject: [PATCH] changes in implicit velocity output: works now also with
branched grids and withgrids where dim != dimworld
reviewed by Bernd
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@14065 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
dumux/implicit/common/implicitvelocityoutput.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dumux/implicit/common/implicitvelocityoutput.hh b/dumux/implicit/common/implicitvelocityoutput.hh
index 5466be50cd..5de58cd233 100644
--- a/dumux/implicit/common/implicitvelocityoutput.hh
+++ b/dumux/implicit/common/implicitvelocityoutput.hh
@@ -218,7 +218,7 @@ public:
elemVolVars);
Scalar flux = fluxVars.volumeFlux(phaseIdx);
- scvVelocities[fIdx] = flux;
+ scvVelocities[fIdx] += flux;
fIdxInner++;
}
@@ -239,7 +239,7 @@ public:
for (int i = 0; i < dim; i++)
refVelocity[i] = 0.5 * (scvVelocities[2*i + 1] - scvVelocities[2*i]);
- Dune::FieldVector<Scalar, dim> scvVelocity(0);
+ Dune::FieldVector<Scalar, dimWorld> scvVelocity(0);
jacobianT2.mtv(refVelocity, scvVelocity);
scvVelocity /= element.geometry().integrationElement(localPos);
--
GitLab