| 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() && |
| 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()) { |
| 269 |
|
if (simParams->getPrintTaggedPairDistance()) { |
| 271 |
|
} |
| 272 |
|
} |
| 273 |
|
|
| 274 |
< |
if (simParams->getUseRNEMD()) { |
| 274 |
> |
if (simParams->getRNEMDParameters()->getUseRNEMD()) { |
| 275 |
|
mask.set(Stats::RNEMD_EXCHANGE_TOTAL); |
| 276 |
|
} |
| 277 |
|
|