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