Commit dfc9c902 authored by Philipp Nuske's avatar Philipp Nuske
Browse files

velocities should be of precision Scalar and not the one of the grid. This...

velocities should be of precision Scalar and not the one of the grid. This makes it possible to use quad in stead of double precision. 
(dunecontrol with  --enable-quad and 
SET_TYPE_PROP(CombustionProblem, Scalar, double ); // quad / double  )

reviewed by Bernd

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@11971 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 3561de44
......@@ -229,11 +229,11 @@ public:
}
}
Dune::FieldVector < CoordScalar, dim > refVelocity(0);
Dune::FieldVector < Scalar, dim > refVelocity(0);
for (int i = 0; i < dim; i++)
refVelocity[i] = 0.5 * (scvVelocities[2*i + 1] - scvVelocities[2*i]);
Dune::FieldVector<CoordScalar, dim> scvVelocity(0);
Dune::FieldVector<Scalar, dim> 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