| 93 |
|
Vector3d pos1 = sd1->getPos(); |
| 94 |
|
Vector3d pos2 = sd2->getPos(); |
| 95 |
|
Vector3d r12 = pos1 - pos2; |
| 96 |
< |
currentSnapshot_->wrapVector(r12); |
| 96 |
> |
if (usePeriodicBoundaryConditions_) |
| 97 |
> |
currentSnapshot_->wrapVector(r12); |
| 98 |
|
Vector3d dipole1 = sd1->getElectroFrame().getColumn(2); |
| 99 |
|
Vector3d dipole2 = sd2->getElectroFrame().getColumn(2); |
| 100 |
|
|
| 108 |
|
|
| 109 |
|
RealType halfBin = (nAngleBins_ - 1) * 0.5; |
| 110 |
|
int angleBin1 = halfBin * (cosAngle1 + 1.0); |
| 111 |
< |
int angleBin2 = halfBin * (cosAngle1 + 1.0); |
| 111 |
> |
int angleBin2 = halfBin * (cosAngle2 + 1.0); |
| 112 |
|
|
| 113 |
< |
++histogram_[angleBin1][angleBin1]; |
| 113 |
> |
++histogram_[angleBin1][angleBin2]; |
| 114 |
|
++npairs_; |
| 115 |
|
} |
| 116 |
|
|