Skip to content
Snippets Groups Projects
Commit 487d0c68 authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

The Dexp method wants a SkewMatrix instead of a vector now 

[[Imported from SVN: r8346]]
parent c8b7d611
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
dune/gfe/rodlocalstiffness.hh
+ 3
− 3
View file @ 487d0c68
......@@ -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);
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment