| 35 |
|
* |
| 36 |
|
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
| 37 |
|
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
| 38 |
< |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). |
| 38 |
> |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). |
| 39 |
|
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
| 40 |
|
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
| 41 |
|
*/ |
| 81 |
|
} |
| 82 |
|
|
| 83 |
|
|
| 84 |
< |
void GofXyz::initalizeHistogram() { |
| 84 |
> |
void GofXyz::initializeHistogram() { |
| 85 |
|
//calculate the center of mass of the molecule of selected stuntdouble in selection1 |
| 86 |
|
|
| 87 |
|
if (!evaluator3_.isDynamic()) { |
| 151 |
|
|
| 152 |
|
Vector3d newR12 = i->second * r12; |
| 153 |
|
// x, y and z's possible values range -halfLen_ to halfLen_ |
| 154 |
< |
int xbin = (newR12.x() + halfLen_) / deltaR_; |
| 155 |
< |
int ybin = (newR12.y() + halfLen_) / deltaR_; |
| 156 |
< |
int zbin = (newR12.z() + halfLen_) / deltaR_; |
| 154 |
> |
int xbin = int( (newR12.x() + halfLen_) / deltaR_); |
| 155 |
> |
int ybin = int( (newR12.y() + halfLen_) / deltaR_); |
| 156 |
> |
int zbin = int( (newR12.z() + halfLen_) / deltaR_); |
| 157 |
|
|
| 158 |
|
if (xbin < nRBins_ && xbin >=0 && |
| 159 |
|
ybin < nRBins_ && ybin >= 0 && |