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