diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh index 3c002e46d5c7388785e857980cf6f7ba2c8c02e6..f60d0c04eaa33b26996c28cc5c08d1f938285294 100644 --- a/dune/gfe/rotation.hh +++ b/dune/gfe/rotation.hh @@ -281,13 +281,14 @@ public: vAtIdentity[2] = 0.5*v.axial()[2]; // multiply with base point to get real embedded tangent vector - Quaternion<T> vQuat = p.mult(vAtIdentity); + Quaternion<T> vQuat = ((Quaternion<T>) p).mult(vAtIdentity); //get basis of the tangent space Dune::FieldMatrix<T,3,4> basis = p.orthonormalFrame(); // transform coordinates TangentVector tang; + tang = 0; basis.mv(vQuat,tang); return tang;