From c3efcc6b8031f13fd4a661eee879eb0cf3dba374 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 29 Mar 2011 16:59:41 +0000
Subject: [PATCH] bugfix: basis vectors of the orthonormal frame were not
normalized
[[Imported from SVN: r7033]]
---
dune/gfe/unitvector.hh | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/dune/gfe/unitvector.hh b/dune/gfe/unitvector.hh
index d71aaf09..4251dae8 100644
--- a/dune/gfe/unitvector.hh
+++ b/dune/gfe/unitvector.hh
@@ -325,6 +325,10 @@ public:
if (data_[N-1] > 0)
for (size_t j=0; j<N-1; j++)
result[j][N-1] *= -1;
+
+ // normalize the rows to make the orthogonal basis orthonormal
+ for (size_t i=0; i<N-1; i++)
+ result[i] /= result[i].two_norm();
return result;
}
--
GitLab