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