diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh index 31a65c6dce564c2a22b0eaae67a22ad65e48d6c8..9f89f6ac861e90035ed987e150f552b0c985c448 100644 --- a/dune/gfe/rotation.hh +++ b/dune/gfe/rotation.hh @@ -454,7 +454,8 @@ public: Dune::FieldMatrix<T,3,3> mat = s.toMatrix(); mat *= 0.5; - Dune::FieldMatrix<T,3,3> skewSquare = mat.rightmultiply(mat); + Dune::FieldMatrix<T,3,3> skewSquare = mat; + skewSquare.rightmultiply(mat); mat += skewSquare; mat *= 2/norm;