From 7c1c5d3a724f53838907c3ef82d42777b0840318 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 14 Jan 2011 14:57:56 +0000
Subject: [PATCH] bugfix: 1e-10 instead of 1e10
[[Imported from SVN: r6756]]
---
dune/gfe/rotation.hh | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh
index c482677d..22290429 100644
--- a/dune/gfe/rotation.hh
+++ b/dune/gfe/rotation.hh
@@ -270,8 +270,7 @@ public:
for (int m=0; m<3; m++) {
-#warning There is a bug in the following line!
- result[m][i] = (norm<1e10)
+ result[m][i] = (norm<1e-10)
/** \todo Isn't there a better way to implement this stably? */
? 0.5 * (i==m)
: 0.5 * std::cos(norm/2) * v[i] * v[m] / (norm*norm) + sincHalf(norm) * ( (i==m) - v[i]*v[m]/(norm*norm));
--
GitLab