diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh index feccfe5562c3954cc18e903bd1cc5d31acceb6f2..a53c959607ad1b5a208f35c7d766f0280c30150e 100644 --- a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh +++ b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh @@ -992,7 +992,7 @@ iterate(std::map<std::pair<std::string,std::string>, RigidBodyMotion<3> >& lambd RigidBodyMotion<3>::TangentVector& correction = interfaceCorrection[interfaceName]; for (int j=0; j<6; j++) correction[j] = (alpha_[0] * continuumCorrection[interfaceName][j] + alpha_[1] * rodCorrection[interfaceName][j]) - / alpha_[0] + alpha_[1]; + / (alpha_[0] + alpha_[1]; } @@ -1409,7 +1409,7 @@ iterateWithContact(std::map<std::pair<std::string,std::string>, RigidBodyMotion< RigidBodyMotion<3>::TangentVector& correction = interfaceCorrection[interfaceName]; for (int j=0; j<6; j++) correction[j] = (alpha_[0] * continuumCorrection[interfaceName][j] + alpha_[1] * rodCorrection[interfaceName][j]) - / alpha_[0] + alpha_[1]; + / (alpha_[0] + alpha_[1]); }