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 |
|
|