diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh index eb70e08e078e1ecbf0edd08724c270cd94bac46f..3a4952e89f1c285c666239a43c1b26e24d4ca645 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++) {