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