| 479 |
|
for (unsigned int j = 0; j < toposForAtom[atom1].size(); j++) { |
| 480 |
|
if (toposForAtom[atom1][j] == atom2) |
| 481 |
|
return topoDist[atom1][j]; |
| 482 |
< |
} |
| 482 |
> |
} |
| 483 |
|
return 0; |
| 484 |
|
} |
| 485 |
|
|
| 939 |
|
|
| 940 |
|
|
| 941 |
|
|
| 942 |
< |
int ForceMatrixDecomposition::getNAtomsInRow() { |
| 942 |
> |
int& ForceMatrixDecomposition::getNAtomsInRow() { |
| 943 |
|
#ifdef IS_MPI |
| 944 |
|
return nAtomsInRow_; |
| 945 |
|
#else |
| 950 |
|
/** |
| 951 |
|
* returns the list of atoms belonging to this group. |
| 952 |
|
*/ |
| 953 |
< |
vector<int> ForceMatrixDecomposition::getAtomsInGroupRow(int cg1){ |
| 953 |
> |
vector<int>& ForceMatrixDecomposition::getAtomsInGroupRow(int cg1){ |
| 954 |
|
#ifdef IS_MPI |
| 955 |
|
return groupListRow_[cg1]; |
| 956 |
|
#else |
| 958 |
|
#endif |
| 959 |
|
} |
| 960 |
|
|
| 961 |
< |
vector<int> ForceMatrixDecomposition::getAtomsInGroupColumn(int cg2){ |
| 961 |
> |
vector<int>& ForceMatrixDecomposition::getAtomsInGroupColumn(int cg2){ |
| 962 |
|
#ifdef IS_MPI |
| 963 |
|
return groupListCol_[cg2]; |
| 964 |
|
#else |
| 981 |
|
return d; |
| 982 |
|
} |
| 983 |
|
|
| 984 |
< |
Vector3d ForceMatrixDecomposition::getGroupVelocityColumn(int cg2){ |
| 984 |
> |
Vector3d& ForceMatrixDecomposition::getGroupVelocityColumn(int cg2){ |
| 985 |
|
#ifdef IS_MPI |
| 986 |
|
return cgColData.velocity[cg2]; |
| 987 |
|
#else |
| 989 |
|
#endif |
| 990 |
|
} |
| 991 |
|
|
| 992 |
< |
Vector3d ForceMatrixDecomposition::getAtomVelocityColumn(int atom2){ |
| 992 |
> |
Vector3d& ForceMatrixDecomposition::getAtomVelocityColumn(int atom2){ |
| 993 |
|
#ifdef IS_MPI |
| 994 |
|
return atomColData.velocity[atom2]; |
| 995 |
|
#else |
| 1027 |
|
return d; |
| 1028 |
|
} |
| 1029 |
|
|
| 1030 |
< |
RealType ForceMatrixDecomposition::getMassFactorRow(int atom1) { |
| 1030 |
> |
RealType& ForceMatrixDecomposition::getMassFactorRow(int atom1) { |
| 1031 |
|
#ifdef IS_MPI |
| 1032 |
|
return massFactorsRow[atom1]; |
| 1033 |
|
#else |
| 1035 |
|
#endif |
| 1036 |
|
} |
| 1037 |
|
|
| 1038 |
< |
RealType ForceMatrixDecomposition::getMassFactorColumn(int atom2) { |
| 1038 |
> |
RealType& ForceMatrixDecomposition::getMassFactorColumn(int atom2) { |
| 1039 |
|
#ifdef IS_MPI |
| 1040 |
|
return massFactorsCol[atom2]; |
| 1041 |
|
#else |
| 1058 |
|
return d; |
| 1059 |
|
} |
| 1060 |
|
|
| 1061 |
< |
vector<int> ForceMatrixDecomposition::getExcludesForAtom(int atom1) { |
| 1061 |
> |
vector<int>& ForceMatrixDecomposition::getExcludesForAtom(int atom1) { |
| 1062 |
|
return excludesForAtom[atom1]; |
| 1063 |
|
} |
| 1064 |
|
|