--- branches/development/src/primitives/GhostBend.cpp 2011/11/22 20:38:56 1665 +++ branches/development/src/primitives/GhostBend.cpp 2012/05/19 13:30:21 1712 @@ -45,7 +45,7 @@ namespace OpenMD { namespace OpenMD { /**@todo still a lot left to improve*/ - void GhostBend::calcForce(RealType& angle) { + void GhostBend::calcForce(RealType& angle, bool doParticlePot) { DirectionalAtom* ghostAtom = static_cast(atom2_); Vector3d pos1 = atom1_->getPos(); @@ -95,10 +95,11 @@ namespace OpenMD { ghostAtom->addFrc(-force1); ghostAtom->addTrq( cross(r23, force3) ); + if(doParticlePot) { + atom1_->addParticlePot(potential_); + ghostAtom->addParticlePot(potential_); + } - atom1_->addParticlePot(potential_); - ghostAtom->addParticlePot(potential_); - angle = theta /M_PI * 180.0; }