| 709 |  | RealType dVdFQ1(0.0); | 
| 710 |  | RealType dVdFQ2(0.0); | 
| 711 |  | potVec longRangePotential(0.0); | 
| 712 | + | potVec reciprocalPotential(0.0); | 
| 713 |  | potVec workPot(0.0); | 
| 714 |  | potVec exPot(0.0); | 
| 715 |  | Vector3d eField1(0.0); | 
| 927 |  | // collects pairwise information | 
| 928 |  | fDecomp_->collectData(); | 
| 929 |  | if (cutoffMethod_ == EWALD_FULL) { | 
| 930 | < | interactionMan_->doReciprocalSpaceSum(); | 
| 930 | > | interactionMan_->doReciprocalSpaceSum(reciprocalPotential); | 
| 931 |  | } | 
| 932 |  |  | 
| 933 |  | if (info_->requiresSelfCorrection()) { | 
| 941 |  | fDecomp_->collectSelfData(); | 
| 942 |  |  | 
| 943 |  | longRangePotential = *(fDecomp_->getEmbeddingPotential()) + | 
| 944 | < | *(fDecomp_->getPairwisePotential()); | 
| 944 | > | *(fDecomp_->getPairwisePotential()) + reciprocalPotential; | 
| 945 |  |  | 
| 946 |  | curSnapshot->setLongRangePotential(longRangePotential); | 
| 947 |  |  | 
| 950 |  |  | 
| 951 |  | } | 
| 952 |  |  | 
| 952 | – |  | 
| 953 |  | void ForceManager::postCalculation() { | 
| 954 |  |  | 
| 955 |  | vector<Perturbation*>::iterator pi; |