| 72 |
|
} |
| 73 |
|
} |
| 74 |
|
|
| 75 |
< |
void GofAngle2::initalizeHistogram() { |
| 75 |
> |
void GofAngle2::initializeHistogram() { |
| 76 |
|
npairs_ = 0; |
| 77 |
|
for (unsigned int i = 0; i < histogram_.size(); ++i) |
| 78 |
|
std::fill(histogram_[i].begin(), histogram_[i].end(), 0); |
| 108 |
|
RealType cosAngle2 = dot(dipole1, dipole2); |
| 109 |
|
|
| 110 |
|
RealType halfBin = (nAngleBins_ - 1) * 0.5; |
| 111 |
< |
int angleBin1 = halfBin * (cosAngle1 + 1.0); |
| 112 |
< |
int angleBin2 = halfBin * (cosAngle2 + 1.0); |
| 111 |
> |
int angleBin1 = int(halfBin * (cosAngle1 + 1.0)); |
| 112 |
> |
int angleBin2 = int(halfBin * (cosAngle2 + 1.0)); |
| 113 |
|
|
| 114 |
|
++histogram_[angleBin1][angleBin2]; |
| 115 |
|
++npairs_; |