diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh
index 055de5ff44420819c42f739176f7384c1367160c..8396bf311ecf66cae5f73563550ae49cbc0ba29c 100644
--- a/dune/gfe/rotation.hh
+++ b/dune/gfe/rotation.hh
@@ -187,15 +187,11 @@ public:
     {
         for (int i=0; i<4; i++)
             (*this)[i] = c[i];
-
-        *this /= this->two_norm();
     }
 
     explicit Rotation<T,3>(const Dune::FieldVector<T,4>& c)
         : Quaternion<T>(c)
-    {
-        *this /= this->two_norm();
-    }
+    {}
 
     Rotation<T,3>(Dune::FieldVector<T,3> axis, T angle)
     {
diff --git a/dune/gfe/unitvector.hh b/dune/gfe/unitvector.hh
index a64902d1583a26c1adbc641b84c599d94c37e0d8..3e3bf8ef51545dc70bef6982d73475a46377e2bc 100644
--- a/dune/gfe/unitvector.hh
+++ b/dune/gfe/unitvector.hh
@@ -89,19 +89,16 @@ public:
     UnitVector()
     {}
 
-    /** \brief Constructor from a vector.  The vector gets normalized */
+    /** \brief Constructor from a vector.  The vector does not get normalized! */
     UnitVector(const Dune::FieldVector<T,N>& vector)
         : data_(vector)
-    {
-        data_ /= data_.two_norm();
-    }
+    {}
 
-    /** \brief Constructor from an array.  The array gets normalized */
+    /** \brief Constructor from an array.  The array does not get normalized! */
     UnitVector(const Dune::array<T,N>& vector)
     {
         for (int i=0; i<N; i++)
             data_[i] = vector[i];
-        data_ /= data_.two_norm();
     }
 
     /** \brief Assigment from UnitVector with different type -- used for automatic differentiation with ADOL-C */
@@ -124,7 +121,6 @@ public:
     UnitVector<T,N>& operator=(const Dune::FieldVector<T,N>& vector)
     {
         data_ = vector;
-        data_ /= data_.two_norm();
         return *this;
     }