ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/trunk/src/UseTheForce/DUFF.cpp
(Generate patch)

Comparing trunk/src/UseTheForce/DUFF.cpp (file contents):
Revision 801 by chuckv, Wed Dec 7 23:10:53 2005 UTC vs.
Revision 1269 by gezelter, Tue Jul 1 13:28:23 2008 UTC

# Line 45 | Line 45
45   #include "UseTheForce/DarkSide/sticky_interface.h"
46   #include "UseTheForce/DarkSide/gb_interface.h"
47   #include "UseTheForce/ForceFieldFactory.hpp"
48 + #include "io/BaseAtomTypesSectionParser.hpp"
49   #include "io/DirectionalAtomTypesSectionParser.hpp"
50   #include "io/AtomTypesSectionParser.hpp"
51   #include "io/LennardJonesAtomTypesSectionParser.hpp"
# Line 67 | Line 68 | namespace oopse {
68      setForceFieldFileName("DUFF2.frc");
69  
70      //The order of adding section parsers is important.
71 +    //OptionSectionParser must come first to set options for other parsers
72      //DirectionalAtomTypesSectionParser should be added before
73      //AtomTypesSectionParser, and these two section parsers will actually
74      //create "real" AtomTypes (AtomTypesSectionParser will create AtomType and
# Line 79 | Line 81 | namespace oopse {
81      //and AtomTypesSectionParser. The order of BondTypesSectionParser,
82      //BendTypesSectionParser and TorsionTypesSectionParser are not important.
83      spMan_.push_back(new OptionSectionParser(forceFieldOptions_));    
84 <    spMan_.push_back(new DirectionalAtomTypesSectionParser());
84 >    spMan_.push_back(new DirectionalAtomTypesSectionParser(forceFieldOptions_));
85 >    spMan_.push_back(new BaseAtomTypesSectionParser());
86      spMan_.push_back(new AtomTypesSectionParser());
87      spMan_.push_back(new LennardJonesAtomTypesSectionParser(forceFieldOptions_));
88 <    spMan_.push_back(new ChargeAtomTypesSectionParser());
89 <    spMan_.push_back(new MultipoleAtomTypesSectionParser());
90 <    spMan_.push_back(new StickyAtomTypesSectionParser());
91 <    spMan_.push_back(new StickyPowerAtomTypesSectionParser());
92 <    spMan_.push_back(new GayBerneAtomTypesSectionParser());
93 <    spMan_.push_back(new BondTypesSectionParser());
94 <    spMan_.push_back(new BendTypesSectionParser());
95 <    spMan_.push_back(new TorsionTypesSectionParser());
88 >    spMan_.push_back(new ChargeAtomTypesSectionParser(forceFieldOptions_));
89 >    spMan_.push_back(new MultipoleAtomTypesSectionParser(forceFieldOptions_));
90 >    spMan_.push_back(new StickyAtomTypesSectionParser(forceFieldOptions_));
91 >    spMan_.push_back(new StickyPowerAtomTypesSectionParser(forceFieldOptions_));
92 >    spMan_.push_back(new GayBerneAtomTypesSectionParser(forceFieldOptions_));
93 >    spMan_.push_back(new BondTypesSectionParser(forceFieldOptions_));
94 >    spMan_.push_back(new BendTypesSectionParser(forceFieldOptions_));
95 >    spMan_.push_back(new TorsionTypesSectionParser(forceFieldOptions_));
96      
97    }
98  
99    void DUFF::parse(const std::string& filename) {
100      ifstrstream* ffStream;
101 +    bool hasGBtypes;
102 +
103      ffStream = openForceFieldFile(filename);
104  
105      spMan_.parse(*ffStream, *this);
# Line 112 | Line 117 | namespace oopse {
117        at->complete();
118      }
119  
120 +    hasGBtypes = false;
121 +    for (at = atomTypeCont_.beginType(i); at != NULL;
122 +         at = atomTypeCont_.nextType(i)) {
123 +      if (at->isGayBerne())
124 +        hasGBtypes = true;
125 +    }
126 +    
127      int isError = 0;
128  
129 +    if (hasGBtypes) {
130 +      completeGBFF(&isError);
131 +    }
132 +
133      delete ffStream;
134      
135    }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines