ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/branches/development/src/applications/dump2Xyz/Dump2XYZ.cpp
(Generate patch)

Comparing branches/development/src/applications/dump2Xyz/Dump2XYZ.cpp (file contents):
Revision 1465 by chuckv, Fri Jul 9 23:08:25 2010 UTC vs.
Revision 1725 by gezelter, Sat May 26 18:13:43 2012 UTC

# Line 36 | Line 36
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   #include <iostream>
# Line 63 | Line 64 | using namespace OpenMD;
64  
65   using namespace OpenMD;
66  
67 + using namespace std;
68   int main(int argc, char* argv[]){
69    
68  //register force fields
69  registerForceFields();
70  
70    gengetopt_args_info args_info;
71 <  std::string dumpFileName;
72 <  std::string xyzFileName;
73 <  bool printVel;
74 <  bool printFrc;
75 <  bool printVec;
76 <  bool printChrg;
71 >  string dumpFileName;
72 >  string xyzFileName;
73 >
74 >  bool printVel(false);
75 >  bool printFrc(false);
76 >  bool printVec(false);
77 >  bool printChrg(false);
78    
79    //parse the command line option
80    if (cmdline_parser (argc, argv, &args_info) != 0) {
# Line 85 | Line 85 | int main(int argc, char* argv[]){
85    if (args_info.input_given){
86      dumpFileName = args_info.input_arg;
87    } else {
88 <    std::cerr << "Does not have input file name" << std::endl;
88 >    cerr << "Does not have input file name" << endl;
89      exit(1);
90    }
91    
# Line 114 | Line 114 | int main(int argc, char* argv[]){
114    SSDAtomVisitor* ssdVisitor = new SSDAtomVisitor(info);
115    compositeVisitor->addVisitor(ssdVisitor, 800);
116    
117 +  //create GBtail atom visitor
118 +  GBtailVisitor* gbtVisitor = new GBtailVisitor(info);
119 +  compositeVisitor->addVisitor(gbtVisitor, 790);
120 +  
121 +  //create GBhead atom visitor
122 +  GBheadVisitor* gbhVisitor = new GBheadVisitor(info);
123 +  compositeVisitor->addVisitor(gbhVisitor, 789);
124 +  
125    //create default atom visitor
126    DefaultAtomVisitor* defaultAtomVisitor = new DefaultAtomVisitor(info);
127    compositeVisitor->addVisitor(defaultAtomVisitor, 700);
# Line 130 | Line 138 | int main(int argc, char* argv[]){
138    if (args_info.basetype_flag) {
139      AtomNameVisitor* atomNameVisitor = new AtomNameVisitor(info);
140      compositeVisitor->addVisitor(atomNameVisitor, 550);    
141 <    std::cout << compositeVisitor->toString();
141 >    cout << compositeVisitor->toString();
142    }
143    
144    //create ZconsVisitor
# Line 172 | Line 180 | int main(int argc, char* argv[]){
180                                                         args_info.refsele_arg),
181                                   250);
182    } else if (args_info.refsele_given || args_info.originsele_given) {
183 <    std::cerr << "Both of --refsele and --originsele should appear by pair"
184 <              << std::endl;
183 >    cerr << "Both of --refsele and --originsele should appear by pair"
184 >              << endl;
185      exit(1);
186    }
187    
# Line 212 | Line 220 | int main(int argc, char* argv[]){
220    DumpReader* dumpReader = new DumpReader(info, dumpFileName);
221    int nframes = dumpReader->getNFrames();
222    
223 <  std::ofstream xyzStream(xyzFileName.c_str());
223 >  ofstream xyzStream(xyzFileName.c_str());
224    
225    SimInfo::MoleculeIterator miter;
226    Molecule::IntegrableObjectIterator  iiter;

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines