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