209 |
|
RealType volume = this->getVolume(); |
210 |
|
Snapshot* curSnapshot = info_->getSnapshotManager()->getCurrentSnapshot(); |
211 |
|
Mat3x3d tau = curSnapshot->statData.getTau(); |
212 |
< |
|
212 |
> |
|
213 |
|
pressureTensor = (p_global + OOPSEConstant::energyConvert* tau)/volume; |
214 |
|
|
215 |
|
return pressureTensor; |
228 |
|
stat[Stats::VOLUME] = getVolume(); |
229 |
|
|
230 |
|
Mat3x3d tensor =getPressureTensor(); |
231 |
< |
stat[Stats::PRESSURE_TENSOR_X] = tensor(0, 0); |
232 |
< |
stat[Stats::PRESSURE_TENSOR_Y] = tensor(1, 1); |
233 |
< |
stat[Stats::PRESSURE_TENSOR_Z] = tensor(2, 2); |
231 |
> |
stat[Stats::PRESSURE_TENSOR_XX] = tensor(0, 0); |
232 |
> |
stat[Stats::PRESSURE_TENSOR_XY] = tensor(0, 1); |
233 |
> |
stat[Stats::PRESSURE_TENSOR_XZ] = tensor(0, 2); |
234 |
> |
stat[Stats::PRESSURE_TENSOR_YX] = tensor(1, 0); |
235 |
> |
stat[Stats::PRESSURE_TENSOR_YY] = tensor(1, 1); |
236 |
> |
stat[Stats::PRESSURE_TENSOR_YZ] = tensor(1, 2); |
237 |
> |
stat[Stats::PRESSURE_TENSOR_ZX] = tensor(2, 0); |
238 |
> |
stat[Stats::PRESSURE_TENSOR_ZY] = tensor(2, 1); |
239 |
> |
stat[Stats::PRESSURE_TENSOR_ZZ] = tensor(2, 2); |
240 |
|
|
241 |
|
|
242 |
|
/**@todo need refactorying*/ |