Commit 09d3f6b1 authored by Natalie Schroeder's avatar Natalie Schroeder
Browse files

changes in implicit velocity output: works now also with branched grids and...

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
parent 526e0456
......@@ -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);
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment