215 |
|
|
216 |
|
void VelocityVerletIntegrator::calcForce() { |
217 |
|
forceMan_->calcForces(); |
218 |
+ |
flucQ_->applyConstraints(); |
219 |
|
} |
220 |
|
|
221 |
|
DumpWriter* VelocityVerletIntegrator::createDumpWriter() { |
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()) { |