| # | Line 72 | Line 72 | namespace oopse { | |
|---|---|---|
| 72 | bondType_->calcForce(len, potential_, dvdr); | |
| 73 | ||
| 74 | force = r12 * (-dvdr / len); | |
| 75 | < | |
| 75 | > | |
| 76 | atom1_->addFrc(-force); | |
| 77 | atom2_->addFrc(force); | |
| 78 | + | |
| 79 | + | atom1_->addParticlePot(potential_); |
| 80 | + | atom2_->addParticlePot(potential_); |
| 81 | } | |
| 82 | ||
| 83 | RealType getPotential() { | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |