From a506a0e797d5273499478e150556cc863206c405 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 31 Jan 2011 08:55:45 +0000
Subject: [PATCH] fix the canonical --> local force transformation
[[Imported from SVN: r6903]]
---
dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh | 1 +
1 file changed, 1 insertion(+)
diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
index eb70e08e..3a4952e8 100644
--- a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
+++ b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
@@ -664,6 +664,7 @@ linearizedRodNeumannToDirichletMap(const std::string& rodName,
/** \todo Why don't we have to transform the force as well? */
Dune::FieldVector<double,3> localForce, localTorque;
orientationMatrix.mv(canonicalTorque,localTorque);
+ localForce = canonicalForce;
RigidBodyMotion<3>::TangentVector localForceTorque;
for (int i=0; i<3; i++) {
--
GitLab