| 98 |
|
thermo.saveStat(); |
| 99 |
|
saveConservedQuantity(); |
| 100 |
|
if (simParams->getUseRNEMD()) |
| 101 |
< |
rnemd_->getStatus(); |
| 101 |
> |
rnemd_->getStarted(); |
| 102 |
|
|
| 103 |
|
statWriter->writeStat(currentSnapshot_->statData); |
| 104 |
|
|
| 109 |
|
currReset = resetTime + currentSnapshot_->getTime(); |
| 110 |
|
} |
| 111 |
|
if (simParams->getUseRNEMD()){ |
| 112 |
< |
currRNEMD = RNEMD_swapTime + currentSnapshot_->getTime(); |
| 112 |
> |
currRNEMD = RNEMD_exchangeTime + currentSnapshot_->getTime(); |
| 113 |
|
} |
| 114 |
|
needPotential = false; |
| 115 |
|
needStress = false; |
| 160 |
|
} |
| 161 |
|
} |
| 162 |
|
if (useRNEMD) { |
| 163 |
+ |
rnemd_->collectData(); |
| 164 |
|
if (currentSnapshot_->getTime() >= currRNEMD) { |
| 165 |
< |
rnemd_->doSwap(); |
| 166 |
< |
currRNEMD += RNEMD_swapTime; |
| 165 |
> |
rnemd_->doRNEMD(); |
| 166 |
> |
currRNEMD += RNEMD_exchangeTime; |
| 167 |
|
} |
| 168 |
|
} |
| 169 |
|
|
| 178 |
|
thermo.saveStat(); |
| 179 |
|
saveConservedQuantity(); |
| 180 |
|
|
| 181 |
< |
if (simParams->getUseRNEMD()) |
| 181 |
> |
if (simParams->getUseRNEMD()) { |
| 182 |
|
rnemd_->getStatus(); |
| 183 |
+ |
} |
| 184 |
|
|
| 185 |
|
statWriter->writeStat(currentSnapshot_->statData); |
| 186 |
|
|
| 268 |
|
} |
| 269 |
|
|
| 270 |
|
if (simParams->getUseRNEMD()) { |
| 271 |
< |
mask.set(Stats::RNEMD_SWAP_TOTAL); |
| 271 |
> |
mask.set(Stats::RNEMD_EXCHANGE_TOTAL); |
| 272 |
|
} |
| 273 |
< |
|
| 273 |
> |
|
| 274 |
|
|
| 275 |
< |
return new StatWriter(info_->getStatFileName(), mask); |
| 275 |
> |
return new StatWriter(info_->getStatFileName(), mask); |
| 276 |
|
} |
| 277 |
|
|
| 278 |
|
|