710 |
|
RealType dVdFQ1(0.0); |
711 |
|
RealType dVdFQ2(0.0); |
712 |
|
potVec longRangePotential(0.0); |
713 |
< |
potVec reciprocalPotential(0.0); |
713 |
> |
RealType reciprocalPotential(0.0); |
714 |
|
potVec workPot(0.0); |
715 |
|
potVec exPot(0.0); |
716 |
|
Vector3d eField1(0.0); |
928 |
|
fDecomp_->collectData(); |
929 |
|
if (cutoffMethod_ == EWALD_FULL) { |
930 |
|
interactionMan_->doReciprocalSpaceSum(reciprocalPotential); |
931 |
+ |
|
932 |
+ |
curSnapshot->setReciprocalPotential(reciprocalPotential); |
933 |
|
} |
934 |
|
|
935 |
|
if (info_->requiresSelfCorrection()) { |
943 |
|
fDecomp_->collectSelfData(); |
944 |
|
|
945 |
|
longRangePotential = *(fDecomp_->getEmbeddingPotential()) + |
946 |
< |
*(fDecomp_->getPairwisePotential()) + reciprocalPotential; |
946 |
> |
*(fDecomp_->getPairwisePotential()); |
947 |
|
|
948 |
|
curSnapshot->setLongRangePotential(longRangePotential); |
949 |
|
|