| 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); |