diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh
index 6684ee404e4a3fe13ea55f9ffd15c739516ada8b..b4e2fce9b25d8e8706195575e0dae1aa2a061af0 100644
--- a/dune/gfe/rotation.hh
+++ b/dune/gfe/rotation.hh
@@ -971,7 +971,7 @@ public:
     {
       Tensor3<T,4,3,3> result;
 
-      Dune::FieldVector<T,3> p;
+      Dune::FieldVector<T,4> p;
 
         // The following equations for the derivation of a unit quaternion from a rotation
         // matrix comes from 'E. Salamin, Application of Quaternions to Computation with