| 833 |
|
#endif |
| 834 |
|
} |
| 835 |
|
|
| 836 |
< |
inline Vector3d ForceMatrixDecomposition::getIntergroupVector(int cg1, |
| 837 |
< |
int cg2){ |
| 836 |
> |
Vector3d ForceMatrixDecomposition::getIntergroupVector(int cg1, |
| 837 |
> |
int cg2){ |
| 838 |
|
|
| 839 |
|
Vector3d d; |
| 840 |
|
#ifdef IS_MPI |
| 913 |
|
|
| 914 |
|
} |
| 915 |
|
|
| 916 |
< |
inline Vector3d ForceMatrixDecomposition::getInteratomicVector(int atom1, |
| 917 |
< |
int atom2){ |
| 916 |
> |
Vector3d ForceMatrixDecomposition::getInteratomicVector(int atom1, |
| 917 |
> |
int atom2){ |
| 918 |
|
Vector3d d; |
| 919 |
|
|
| 920 |
|
#ifdef IS_MPI |
| 1512 |
|
whichCell.y() = nCells_.y() * scaled.y(); |
| 1513 |
|
whichCell.z() = nCells_.z() * scaled.z(); |
| 1514 |
|
|
| 1515 |
– |
// find single index of this cell: |
| 1516 |
– |
int m1 = Vlinear(whichCell, nCells_); |
| 1517 |
– |
|
| 1515 |
|
for (vector<Vector3i>::iterator os = cellOffsets_.begin(); |
| 1516 |
|
os != cellOffsets_.end(); ++os) { |
| 1517 |
|
|