| 222 |
|
} |
| 223 |
|
|
| 224 |
|
StatWriter* VelocityVerletIntegrator::createStatWriter() { |
| 225 |
+ |
|
| 226 |
+ |
std::string statFileFormatString = simParams->getStatFileFormat(); |
| 227 |
+ |
StatsBitSet mask = parseStatFileFormat(statFileFormatString); |
| 228 |
+ |
|
| 229 |
|
// if solidThermInt is true, add extra information to the statfile |
| 230 |
|
if (simParams->getUseSolidThermInt()){ |
| 227 |
– |
StatsBitSet mask; |
| 228 |
– |
mask.set(Stats::TIME); |
| 229 |
– |
mask.set(Stats::TOTAL_ENERGY); |
| 230 |
– |
mask.set(Stats::POTENTIAL_ENERGY); |
| 231 |
– |
mask.set(Stats::KINETIC_ENERGY); |
| 232 |
– |
mask.set(Stats::TEMPERATURE); |
| 233 |
– |
mask.set(Stats::PRESSURE); |
| 234 |
– |
mask.set(Stats::VOLUME); |
| 235 |
– |
mask.set(Stats::CONSERVED_QUANTITY); |
| 231 |
|
mask.set(Stats::VRAW); |
| 232 |
|
mask.set(Stats::VHARM); |
| 238 |
– |
return new StatWriter(info_->getStatFileName(), mask); |
| 233 |
|
} |
| 234 |
|
|
| 235 |
|
if (simParams->havePrintPressureTensor() && simParams->getPrintPressureTensor()){ |
| 242 |
– |
StatsBitSet mask; |
| 243 |
– |
mask.set(Stats::TIME); |
| 244 |
– |
mask.set(Stats::TOTAL_ENERGY); |
| 245 |
– |
mask.set(Stats::POTENTIAL_ENERGY); |
| 246 |
– |
mask.set(Stats::KINETIC_ENERGY); |
| 247 |
– |
mask.set(Stats::TEMPERATURE); |
| 248 |
– |
mask.set(Stats::PRESSURE); |
| 249 |
– |
mask.set(Stats::VOLUME); |
| 250 |
– |
mask.set(Stats::CONSERVED_QUANTITY); |
| 236 |
|
mask.set(Stats::PRESSURE_TENSOR_X); |
| 237 |
|
mask.set(Stats::PRESSURE_TENSOR_Y); |
| 238 |
|
mask.set(Stats::PRESSURE_TENSOR_Z); |
| 254 |
– |
return new StatWriter(info_->getStatFileName(), mask); |
| 239 |
|
} |
| 240 |
|
|
| 241 |
< |
return new StatWriter(info_->getStatFileName()); |
| 241 |
> |
return new StatWriter(info_->getStatFileName(), mask); |
| 242 |
|
} |
| 243 |
|
|
| 244 |
|
RestWriter* VelocityVerletIntegrator::createRestWriter(){ |
