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). |
39 |
< |
* [4] Vardeman & Gezelter, in progress (2009). |
39 |
> |
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
40 |
> |
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
41 |
|
*/ |
42 |
|
|
43 |
|
/** |
78 |
|
chi_(0.0), |
79 |
|
integralOfChiDt_(0.0), |
80 |
|
eta_(0.0), id_(-1), hasCOM_(false), |
81 |
< |
hasVolume_(false), volume_(0.0) { |
81 |
> |
volume_(0.0), hasVolume_(false) { |
82 |
|
|
83 |
|
} |
84 |
|
|
89 |
|
currentTime_(0), orthoTolerance_(1e-6), |
90 |
|
orthoRhombic_(0), chi_(0.0), |
91 |
|
integralOfChiDt_(0.0), eta_(0.0), id_(-1), |
92 |
< |
hasCOM_(false), hasVolume_(false), |
93 |
< |
volume_(0.0) { |
92 |
> |
hasCOM_(false), volume_(0.0), |
93 |
> |
hasVolume_(false) { |
94 |
|
} |
95 |
|
|
96 |
|
/** Returns the id of this Snapshot */ |
150 |
|
|
151 |
|
/** Wrapping the vector according to periodic boundary condition*/ |
152 |
|
void wrapVector(Vector3d& v); |
153 |
+ |
/** Scaling a vector to multiples of the periodic box */ |
154 |
+ |
Vector3d scaleVector(Vector3d &v); |
155 |
+ |
|
156 |
+ |
|
157 |
|
Vector3d getCOM(); |
158 |
|
Vector3d getCOMvel(); |
159 |
|
Vector3d getCOMw(); |
201 |
|
eta_ = eta; |
202 |
|
} |
203 |
|
|
204 |
+ |
Mat3x3d getTau() { |
205 |
+ |
return tau_; |
206 |
+ |
} |
207 |
+ |
|
208 |
+ |
void setTau(const Mat3x3d& tau) { |
209 |
+ |
tau_ = tau; |
210 |
+ |
} |
211 |
+ |
|
212 |
|
bool hasCOM() { |
213 |
|
return hasCOM_; |
214 |
|
} |
243 |
|
int id_; /**< identification number of the snapshot */ |
244 |
|
bool hasCOM_; |
245 |
|
bool hasVolume_; |
246 |
< |
|
246 |
> |
Mat3x3d tau_; |
247 |
|
}; |
248 |
|
|
249 |
|
typedef DataStorage (Snapshot::*DataStoragePointer); |