| 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 |
|
|
| 257 |
|
mask.set(Stats::BOX_DIPOLE_X); |
| 258 |
|
mask.set(Stats::BOX_DIPOLE_Y); |
| 259 |
|
mask.set(Stats::BOX_DIPOLE_Z); |
| 260 |
+ |
} |
| 261 |
+ |
|
| 262 |
+ |
if (simParams->getPrintHeatFlux()) { |
| 263 |
+ |
mask.set(Stats::HEATFLUX_X); |
| 264 |
+ |
mask.set(Stats::HEATFLUX_Y); |
| 265 |
+ |
mask.set(Stats::HEATFLUX_Z); |
| 266 |
|
} |
| 267 |
|
|
| 268 |
|
if (simParams->haveTaggedAtomPair() && simParams->havePrintTaggedPairDistance()) { |
| 271 |
|
} |
| 272 |
|
} |
| 273 |
|
|
| 274 |
< |
if (simParams->getUseRNEMD()) { |
| 274 |
> |
if (simParams->getRNEMDParameters()->getUseRNEMD()) { |
| 275 |
|
mask.set(Stats::RNEMD_EXCHANGE_TOTAL); |
| 276 |
|
} |
| 277 |
|
|