35 |
|
* |
36 |
|
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
37 |
|
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
38 |
< |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). |
38 |
> |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). |
39 |
|
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
40 |
|
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
41 |
|
*/ |
94 |
|
frameData.id = -1; |
95 |
|
frameData.currentTime = 0; |
96 |
|
frameData.hmat = Mat3x3d(0.0); |
97 |
< |
frameData.invHmat = Mat3x3d(0.0); |
97 |
> |
frameData.invHmat = Mat3x3d(0.0); |
98 |
> |
frameData.bBox = Mat3x3d(0.0); |
99 |
> |
frameData.invBbox = Mat3x3d(0.0); |
100 |
|
frameData.orthoRhombic = false; |
101 |
|
frameData.bondPotential = 0.0; |
102 |
|
frameData.bendPotential = 0.0; |
133 |
|
frameData.electronicTemperature = 0.0; |
134 |
|
frameData.COM = V3Zero; |
135 |
|
frameData.COMvel = V3Zero; |
136 |
< |
frameData.COMw = V3Zero; |
136 |
> |
frameData.COMw = V3Zero; |
137 |
|
|
138 |
|
hasTotalEnergy = false; |
139 |
|
hasTranslationalKineticEnergy = false; |
254 |
|
/** Returns the inverse H-Matrix */ |
255 |
|
Mat3x3d Snapshot::getInvHmat() { |
256 |
|
return frameData.invHmat; |
257 |
+ |
} |
258 |
+ |
|
259 |
+ |
/** Returns the Bounding Box */ |
260 |
+ |
Mat3x3d Snapshot::getBoundingBox() { |
261 |
+ |
return frameData.bBox; |
262 |
|
} |
263 |
|
|
264 |
+ |
/** Sets the Bounding Box */ |
265 |
+ |
void Snapshot::setBoundingBox(const Mat3x3d& m) { |
266 |
+ |
frameData.bBox = m; |
267 |
+ |
frameData.invBbox = frameData.bBox.inverse(); |
268 |
+ |
} |
269 |
+ |
|
270 |
+ |
/** Returns the inverse Bounding Box */ |
271 |
+ |
Mat3x3d Snapshot::getInvBoundingBox() { |
272 |
+ |
return frameData.invBbox; |
273 |
+ |
} |
274 |
+ |
|
275 |
|
RealType Snapshot::getXYarea() { |
276 |
|
if (!hasXYarea) { |
277 |
|
Vector3d x = frameData.hmat.getColumn(0); |