34 |
|
* |
35 |
|
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
36 |
|
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
37 |
< |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). |
38 |
< |
* [4] Vardeman & Gezelter, in progress (2009). |
37 |
> |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). |
38 |
> |
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
39 |
> |
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
40 |
|
*/ |
41 |
|
|
42 |
|
#include "applications/staticProps/NanoLength.hpp" |
54 |
|
NanoLength::NanoLength(SimInfo* info, |
55 |
|
const std::string& filename, |
56 |
|
const std::string& sele) |
57 |
< |
: StaticAnalyser(info, filename), selectionScript_(sele), evaluator_(info), seleMan_(info) { |
57 |
> |
: StaticAnalyser(info, filename), selectionScript_(sele), evaluator_(info), |
58 |
> |
seleMan_(info) { |
59 |
|
setOutputName(getPrefix(filename) + ".length"); |
60 |
|
|
61 |
|
osq.open(getOutputFileName().c_str()); |
69 |
|
|
70 |
|
void NanoLength::process() { |
71 |
|
Molecule* mol; |
70 |
– |
Atom* atom; |
72 |
|
RigidBody* rb; |
72 |
– |
int myIndex; |
73 |
|
SimInfo::MoleculeIterator mi; |
74 |
|
Molecule::RigidBodyIterator rbIter; |
75 |
– |
Molecule::AtomIterator ai; |
75 |
|
StuntDouble* sd; |
76 |
|
Vector3d vec; |
77 |
< |
int i,j; |
77 |
> |
int i; |
78 |
|
|
79 |
|
DumpReader reader(info_, dumpFilename_); |
80 |
|
int nFrames = reader.getNFrames(); |
119 |
|
osq << time << "\t" << rodLength << std::endl; |
120 |
|
} |
121 |
|
} |
122 |
+ |
osq.close(); |
123 |
|
} |
124 |
|
|
125 |
|
RealType NanoLength::getLength(std::vector<StuntDouble*> atoms) { |