diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh index 42e544798301cd15e2271b0dbfeae35e1d24f7e7..51c9231b436f8283c12218ebc18762c85f4709a6 100644 --- a/dune/gfe/rotation.hh +++ b/dune/gfe/rotation.hh @@ -244,6 +244,8 @@ public: */ static Rotation<3,T> exp(const Rotation<3,T>& p, const EmbeddedTangentVector& v) { + assert( std::fabs(p*v) < 1e-8 ); + // The vector v as a quaternion Quaternion<T> vQuat(v);