100 |
|
//save statistics, before writeStat, we must save statistics |
101 |
|
thermo.saveStat(); |
102 |
|
saveConservedQuantity(); |
103 |
< |
if (simParams->getUseRNEMD()) |
103 |
> |
if (simParams->getRNEMDParameters()->getUseRNEMD()) |
104 |
|
rnemd_->getStarted(); |
105 |
|
|
106 |
|
statWriter->writeStat(currentSnapshot_->statData); |
111 |
|
if (needReset) { |
112 |
|
currReset = resetTime + currentSnapshot_->getTime(); |
113 |
|
} |
114 |
< |
if (simParams->getUseRNEMD()){ |
114 |
> |
if (simParams->getRNEMDParameters()->getUseRNEMD()){ |
115 |
|
currRNEMD = RNEMD_exchangeTime + currentSnapshot_->getTime(); |
116 |
|
} |
117 |
|
needPotential = false; |
174 |
|
thermo.saveStat(); |
175 |
|
saveConservedQuantity(); |
176 |
|
|
177 |
< |
if (simParams->getUseRNEMD()) { |
177 |
> |
if (simParams->getRNEMDParameters()->getUseRNEMD()) { |
178 |
|
rnemd_->getStatus(); |
179 |
|
} |
180 |
|
|
232 |
|
|
233 |
|
|
234 |
|
if (simParams->getUseThermodynamicIntegration()) |
235 |
< |
mask.set(Stats::VRAW); |
235 |
> |
mask.set(Stats::RAW_POTENTIAL); |
236 |
|
|
237 |
|
// if we've got restraints turned on, we'll also want a report of the |
238 |
|
// total harmonic restraints |
239 |
|
if (simParams->getUseRestraints()){ |
240 |
< |
mask.set(Stats::VHARM); |
240 |
> |
mask.set(Stats::RESTRAINT_POTENTIAL); |
241 |
|
} |
242 |
|
|
243 |
|
if (simParams->havePrintPressureTensor() && |
271 |
|
} |
272 |
|
} |
273 |
|
|
274 |
< |
if (simParams->getUseRNEMD()) { |
274 |
> |
if (simParams->getRNEMDParameters()->getUseRNEMD()) { |
275 |
|
mask.set(Stats::RNEMD_EXCHANGE_TOTAL); |
276 |
|
} |
277 |
|
|