From 487d0c6806ec0ce77e12a3c22f45b5af80e26ba4 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 4 Jan 2012 11:02:00 +0000
Subject: [PATCH] The Dexp method wants a SkewMatrix instead of a vector now
[[Imported from SVN: r8346]]
---
dune/gfe/rodlocalstiffness.hh | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/dune/gfe/rodlocalstiffness.hh b/dune/gfe/rodlocalstiffness.hh
index cf494ec9..f518bad2 100644
--- a/dune/gfe/rodlocalstiffness.hh
+++ b/dune/gfe/rodlocalstiffness.hh
@@ -242,7 +242,7 @@ interpolationDerivative(const Rotation<RT,3>& q0, const Rotation<RT,3>& q1, doub
Dune::FieldVector<RT,3> v = Rotation<RT,3>::expInv(q1InvQ0);
v *= (1-s);
- Dune::FieldMatrix<RT,4,3> dExp_v = Rotation<RT,3>::Dexp(v);
+ Dune::FieldMatrix<RT,4,3> dExp_v = Rotation<RT,3>::Dexp(SkewMatrix<RT,3>(v));
Dune::FieldMatrix<RT,3,4> dExpInv = Rotation<RT,3>::DexpInv(q1InvQ0);
@@ -279,7 +279,7 @@ interpolationDerivative(const Rotation<RT,3>& q0, const Rotation<RT,3>& q1, doub
Dune::FieldVector<RT,3> v = Rotation<RT,3>::expInv(q0InvQ1);
v *= s;
- Dune::FieldMatrix<RT,4,3> dExp_v = Rotation<RT,3>::Dexp(v);
+ Dune::FieldMatrix<RT,4,3> dExp_v = Rotation<RT,3>::Dexp(SkewMatrix<RT,3>(v));
Dune::FieldMatrix<RT,3,4> dExpInv = Rotation<RT,3>::DexpInv(q0InvQ1);
@@ -324,7 +324,7 @@ interpolationVelocityDerivative(const Rotation<RT,3>& q0, const Rotation<RT,3>&
Dune::FieldVector<RT,3> v = Rotation<RT,3>::expInv(q0Inv.mult(q1));
v *= s/intervalLength;
- Dune::FieldMatrix<RT,4,3> dExp_v = Rotation<RT,3>::Dexp(v);
+ Dune::FieldMatrix<RT,4,3> dExp_v = Rotation<RT,3>::Dexp(SkewMatrix<RT,3>(v));
Dune::array<Dune::FieldMatrix<RT,3,3>, 4> ddExp;
Rotation<RT,3>::DDexp(v, ddExp);
--
GitLab