From f6970a6998e1bec41d9fd8bc3ce442ffe606ce56 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 1 Jan 2012 17:43:28 +0000
Subject: [PATCH] bugfix: Dexp wants a SkewMatrix as its argument, not a
FieldVector
[[Imported from SVN: r8333]]
---
dune/gfe/rotation.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh
index 9a2d4b98..f05cb2c2 100644
--- a/dune/gfe/rotation.hh
+++ b/dune/gfe/rotation.hh
@@ -449,7 +449,7 @@ public:
Dune::FieldVector<T,3> v = expInv(q);
// The derivative of exp at v
- Dune::FieldMatrix<T,4,3> A = Dexp(v);
+ Dune::FieldMatrix<T,4,3> A = Dexp(SkewMatrix<T,3>(v));
// Compute the Moore-Penrose pseudo inverse A^+ = (A^T A)^{-1} A^T
Dune::FieldMatrix<T,3,3> ATA;
--
GitLab