79 |
|
Vector3d COM; /**< location of center of mass */ |
80 |
|
Vector3d COMvel; /**< system center of mass velocity */ |
81 |
|
Vector3d COMw; /**< system center of mass angular velocity */ |
82 |
< |
Mat3x3d tau; /**< stress tensor */ |
82 |
> |
Mat3x3d stressTensor; /**< stress tensor */ |
83 |
|
Mat3x3d pressureTensor; /**< pressure tensor */ |
84 |
|
Vector3d systemDipole; /**< total system dipole moment */ |
85 |
+ |
Vector3d conductiveHeatFlux; /**< heat flux vector (conductive only) */ |
86 |
|
}; |
87 |
|
|
88 |
|
|
122 |
|
frameData.COM = V3Zero; |
123 |
|
frameData.COMvel = V3Zero; |
124 |
|
frameData.COMw = V3Zero; |
125 |
< |
frameData.tau = Mat3x3d(0.0); |
125 |
> |
frameData.stressTensor = Mat3x3d(0.0); |
126 |
|
frameData.pressureTensor = Mat3x3d(0.0); |
127 |
|
frameData.systemDipole = V3Zero; |
128 |
+ |
frameData.conductiveHeatFlux = V3Zero; |
129 |
|
} |
130 |
|
|
131 |
|
Snapshot(int nAtoms, int nRigidbodies, int nCutoffGroups, |
155 |
|
frameData.COM = V3Zero; |
156 |
|
frameData.COMvel = V3Zero; |
157 |
|
frameData.COMw = V3Zero; |
158 |
< |
frameData.tau = Mat3x3d(0.0); |
158 |
> |
frameData.stressTensor = Mat3x3d(0.0); |
159 |
|
frameData.pressureTensor = Mat3x3d(0.0); |
160 |
|
frameData.systemDipole = V3Zero; |
161 |
+ |
frameData.conductiveHeatFlux = V3Zero; |
162 |
|
} |
163 |
|
|
164 |
|
/** Returns the id of this Snapshot */ |
285 |
|
frameData.eta = eta; |
286 |
|
} |
287 |
|
|
288 |
< |
Mat3x3d getTau() { |
289 |
< |
return frameData.tau; |
288 |
> |
Mat3x3d getStressTensor() { |
289 |
> |
return frameData.stressTensor; |
290 |
|
} |
291 |
|
|
292 |
< |
void setTau(const Mat3x3d& tau) { |
293 |
< |
frameData.tau = tau; |
292 |
> |
void setStressTensor(const Mat3x3d& stressTensor) { |
293 |
> |
frameData.stressTensor = stressTensor; |
294 |
|
} |
295 |
|
|
296 |
+ |
Vector3d getConductiveHeatFlux() { |
297 |
+ |
return frameData.conductiveHeatFlux; |
298 |
+ |
} |
299 |
+ |
|
300 |
+ |
void setConductiveHeatFlux(const Vector3d& heatFlux) { |
301 |
+ |
frameData.conductiveHeatFlux = heatFlux; |
302 |
+ |
} |
303 |
+ |
|
304 |
|
bool hasCOM() { |
305 |
|
return hasCOM_; |
306 |
|
} |