| 488 |
|
pairwisePot = 0.0; |
| 489 |
|
embeddingPot = 0.0; |
| 490 |
|
excludedPot = 0.0; |
| 491 |
+ |
excludedSelfPot = 0.0; |
| 492 |
|
|
| 493 |
|
#ifdef IS_MPI |
| 494 |
|
if (storageLayout_ & DataStorage::dslForce) { |
| 908 |
|
MPI::COMM_WORLD.Allreduce(&ploc1, &ploc2, 1, MPI::REALTYPE, MPI::SUM); |
| 909 |
|
embeddingPot[ii] = ploc2; |
| 910 |
|
} |
| 911 |
+ |
for (int ii = 0; ii < N_INTERACTION_FAMILIES; ii++) { |
| 912 |
+ |
RealType ploc1 = excludedSelfPot[ii]; |
| 913 |
+ |
RealType ploc2 = 0.0; |
| 914 |
+ |
MPI::COMM_WORLD.Allreduce(&ploc1, &ploc2, 1, MPI::REALTYPE, MPI::SUM); |
| 915 |
+ |
excludedSelfPot[ii] = ploc2; |
| 916 |
+ |
} |
| 917 |
|
#endif |
| 918 |
|
|
| 919 |
|
} |