From d5e5afd21da82294d61e072c026562d2b9942a83 Mon Sep 17 00:00:00 2001
From: Jonathan Youett <youett@mi.fu-berlin.de>
Date: Thu, 3 Nov 2011 17:52:47 +0000
Subject: [PATCH] missed a line
[[Imported from SVN: r8111]]
---
dune/gfe/rotation.hh | 2 ++
1 file changed, 2 insertions(+)
diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh
index b1d8ec3a..dd038f07 100644
--- a/dune/gfe/rotation.hh
+++ b/dune/gfe/rotation.hh
@@ -531,6 +531,8 @@ public:
// Make sure we do the right thing if a and b are not in the same sheet
// of the double covering of the unit quaternions over SO(3)
+ T dist = 2*std::acos( std::min(diff[3],1.0) );
+
if (dist>=M_PI)
diff *= -1;
--
GitLab