| 106 |  | Vector3d r3 =sd3->getPos(); | 
| 107 |  | Vector3d r1 = sd1->getPos(); | 
| 108 |  | Vector3d v1 =  r3 - r1; | 
| 109 | < | info_->getSnapshotManager()->getCurrentSnapshot()->wrapVector(v1); | 
| 109 | > | if (usePeriodicBoundaryConditions_) | 
| 110 | > | info_->getSnapshotManager()->getCurrentSnapshot()->wrapVector(v1); | 
| 111 |  | Vector3d zaxis = sd1->getElectroFrame().getColumn(2); | 
| 112 |  | Vector3d xaxis = cross(v1, zaxis); | 
| 113 |  | Vector3d yaxis = cross(zaxis, xaxis); | 
| 131 |  | Vector3d pos1 = sd1->getPos(); | 
| 132 |  | Vector3d pos2 = sd2->getPos(); | 
| 133 |  | Vector3d r12 = pos2 - pos1; | 
| 134 | < | currentSnapshot_->wrapVector(r12); | 
| 134 | > | if (usePeriodicBoundaryConditions_) | 
| 135 | > | currentSnapshot_->wrapVector(r12); | 
| 136 |  |  | 
| 137 |  | std::map<int, RotMat3x3d>::iterator i = rotMats_.find(sd1->getGlobalIndex()); | 
| 138 |  | assert(i != rotMats_.end()); |