diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh index 9a2d4b989ee67c5f0b6f8f672fd4584131ca54c9..f05cb2c229f96b07985a9d0243b3943aa09618e8 100644 --- a/dune/gfe/rotation.hh +++ b/dune/gfe/rotation.hh @@ -449,7 +449,7 @@ public: Dune::FieldVector<T,3> v = expInv(q); // The derivative of exp at v - Dune::FieldMatrix<T,4,3> A = Dexp(v); + Dune::FieldMatrix<T,4,3> A = Dexp(SkewMatrix<T,3>(v)); // Compute the Moore-Penrose pseudo inverse A^+ = (A^T A)^{-1} A^T Dune::FieldMatrix<T,3,3> ATA;