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/BaseAtomTypesSectionParser.hpp" |
51 |
|
#include "io/AtomTypesSectionParser.hpp" |
52 |
+ |
#include "io/BaseAtomTypesSectionParser.hpp" |
53 |
|
#include "io/LennardJonesAtomTypesSectionParser.hpp" |
54 |
|
#include "io/ChargeAtomTypesSectionParser.hpp" |
55 |
|
#include "io/MultipoleAtomTypesSectionParser.hpp" |
84 |
|
//BendTypesSectionParser and TorsionTypesSectionParser are not important. |
85 |
|
spMan_.push_back(new OptionSectionParser(forceFieldOptions_)); |
86 |
|
spMan_.push_back(new DirectionalAtomTypesSectionParser(forceFieldOptions_)); |
87 |
+ |
spMan_.push_back(new BaseAtomTypesSectionParser()); |
88 |
|
spMan_.push_back(new AtomTypesSectionParser()); |
89 |
|
spMan_.push_back(new LennardJonesAtomTypesSectionParser(forceFieldOptions_)); |
90 |
|
spMan_.push_back(new ChargeAtomTypesSectionParser(forceFieldOptions_)); |
100 |
|
|
101 |
|
void DUFF::parse(const std::string& filename) { |
102 |
|
ifstrstream* ffStream; |
103 |
+ |
bool hasGBtypes; |
104 |
+ |
|
105 |
|
ffStream = openForceFieldFile(filename); |
106 |
|
|
107 |
|
spMan_.parse(*ffStream, *this); |
119 |
|
at->complete(); |
120 |
|
} |
121 |
|
|
122 |
+ |
hasGBtypes = false; |
123 |
+ |
for (at = atomTypeCont_.beginType(i); at != NULL; |
124 |
+ |
at = atomTypeCont_.nextType(i)) { |
125 |
+ |
if (at->isGayBerne()) |
126 |
+ |
hasGBtypes = true; |
127 |
+ |
} |
128 |
+ |
|
129 |
|
int isError = 0; |
130 |
|
|
131 |
+ |
if (hasGBtypes) { |
132 |
+ |
completeGBFF(&isError); |
133 |
+ |
} |
134 |
+ |
|
135 |
|
delete ffStream; |
136 |
|
|
137 |
|
} |