| # | Line 213 | Line 213 | namespace OpenMD { | |
|---|---|---|
| 213 | Vector3d vec = j->first->getPos() - j->second->getPos(); | |
| 214 | if (usePeriodicBoundaryConditions_) | |
| 215 | currentSnapshot_->wrapVector(vec); | |
| 216 | < | vec.normalize(); |
| 217 | < | |
| 216 | > | vec.normalize(); |
| 217 | angle += acos(dot(vec, director)) ; | |
| 218 | } | |
| 219 | angle = angle / (sdPairs_.size() * NumericConstant::PI) * 180.0; | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |