diff --git a/dune/gfe/unitvector.hh b/dune/gfe/unitvector.hh index 960a0a4595d80c4f182161d4a1fe55408aad3792..a0f71d813f2cd482bae662b22b7997ce55ee66a4 100644 --- a/dune/gfe/unitvector.hh +++ b/dune/gfe/unitvector.hh @@ -156,7 +156,7 @@ public: // Compute vector A (see notes) Dune::FieldMatrix<double,1,dim> row; - row[0] = a.globalCoordinates(); + row[0] = b.projectOntoTangentSpace(a.globalCoordinates()); row *= secondDerivativeOfArcCosSquared(sp); Dune::FieldMatrix<double,dim,1> column;