| 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(longRangePotential); |
| 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 |
|
|