Commit 550dcd17 authored by Bernd Flemisch's avatar Bernd Flemisch
Browse files

[box] simplify calculation of the boundary face normal

For ages, the method centerUnitOuterNormal exists which does not require
a local coordinate on the intersection. Therefore, we can also get rid
of the assumption.

Reviewed by gruenich.

git-svn-id: svn:// 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 60a7e90d
......@@ -848,9 +848,7 @@ public:
bFace.i = scvIdx;
bFace.j = scvIdx;
// ASSUME constant normal
Dune::FieldVector<CoordScalar, dim-1> localDimM1(0);
bFace.normal = isIt->unitOuterNormal(localDimM1);
bFace.normal = isIt->centerUnitOuterNormal();
bFace.normal *= bFace.area;
typedef Dune::FieldVector< Scalar, 1 > ShapeValue;
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