diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh index 51c9231b436f8283c12218ebc18762c85f4709a6..e002b6cac1a2ccc970ee5d26d4449d27aa3b4876 100644 --- a/dune/gfe/rotation.hh +++ b/dune/gfe/rotation.hh @@ -125,7 +125,7 @@ class Rotation<3,T> : public Quaternion<T> /** \brief Computes sin(x/2) / x without getting unstable for small x */ static T sincHalf(const T& x) { - return (x < 1e-4) ? 0.5 + (x*x/48) : std::sin(x/2)/x; + return (x < 1e-4) ? 0.5 - (x*x/48) : std::sin(x/2)/x; } /** \brief Compute the derivative of arccos^2 without getting unstable for x close to 1 */