| 106 |
|
Vector3d pos1 = sd1->getPos(); |
| 107 |
|
Vector3d pos2 = sd2->getPos(); |
| 108 |
|
Vector3d r12 = pos2 - pos1; |
| 109 |
< |
currentSnapshot_->wrapVector(r12); |
| 109 |
> |
if (usePeriodicBoundaryConditions_) |
| 110 |
> |
currentSnapshot_->wrapVector(r12); |
| 111 |
|
|
| 112 |
|
RealType distance = r12.length(); |
| 113 |
|
int whichRBin = distance / deltaR_; |
| 154 |
|
Vector3d pos1 = sd1->getPos(); |
| 155 |
|
Vector3d pos2 = sd2->getPos(); |
| 156 |
|
Vector3d r12 = pos2 - pos1; |
| 157 |
< |
currentSnapshot_->wrapVector(r12); |
| 157 |
> |
|
| 158 |
> |
if (usePeriodicBoundaryConditions_) |
| 159 |
> |
currentSnapshot_->wrapVector(r12); |
| 160 |
> |
|
| 161 |
|
r12.normalize(); |
| 162 |
|
Vector3d dipole = sd1->getElectroFrame().getColumn(2); |
| 163 |
|
dipole.normalize(); |