| 35 |  | * | 
| 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). | 
| 38 | > | * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). | 
| 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 |  | #ifndef SELECTION_SELECTIONEVALUATOR_HPP | 
| 55 |  | #include "selection/DistanceFinder.hpp" | 
| 56 |  | #include "selection/HullFinder.hpp" | 
| 57 |  | #include "selection/IndexFinder.hpp" | 
| 58 | < | #include "utils/OpenMDBitSet.hpp" | 
| 58 | > | #include "selection/SelectionSet.hpp" | 
| 59 |  | #include "primitives/StuntDouble.hpp" | 
| 60 |  | #include "utils/StringUtils.hpp" | 
| 61 |  | namespace OpenMD { | 
| 73 |  | bool loadScriptString(const std::string& script); | 
| 74 |  | bool loadScriptFile(const std::string& filename); | 
| 75 |  |  | 
| 76 | < | OpenMDBitSet evaluate(); | 
| 76 | > | SelectionSet evaluate(); | 
| 77 | > | SelectionSet evaluate(int frame); | 
| 78 |  |  | 
| 79 |  | /** | 
| 80 |  | * Tests if the result from evaluation of script is dynamic. | 
| 105 |  | } | 
| 106 |  | return script.substr(ichBegin, ichEnd); | 
| 107 |  | } | 
| 108 | + | bool hasSurfaceArea() { return hasSurfaceArea_; } | 
| 109 | + | RealType getSurfaceArea() { | 
| 110 | + | if (hasSurfaceArea_) { | 
| 111 | + | return surfaceArea_; | 
| 112 | + | } else { | 
| 113 | + | sprintf( painCave.errMsg, | 
| 114 | + | "SelectionEvaluator Error: %s\n", "No Surface Area For You!"); | 
| 115 | + | painCave.severity = OPENMD_ERROR; | 
| 116 | + | painCave.isFatal = 1; | 
| 117 | + | simError(); | 
| 118 | + | return 0.0; | 
| 119 | + | } | 
| 120 | + | } | 
| 121 | + |  | 
| 122 |  |  | 
| 123 |  | private: | 
| 124 |  |  | 
| 128 |  |  | 
| 129 |  | bool loadScriptFileInternal(const std::string& filename); | 
| 130 |  |  | 
| 131 | < |  | 
| 131 | > | SelectionSet createSelectionSets(); | 
| 132 |  | void clearDefinitionsAndLoadPredefined(); | 
| 133 |  |  | 
| 134 |  | void define(); | 
| 135 | < | void select(OpenMDBitSet& bs); | 
| 135 | > | void select(SelectionSet& bs); | 
| 136 | > | void select(SelectionSet& bs, int frame); | 
| 137 |  | void predefine(const std::string& script); | 
| 138 |  |  | 
| 139 | < | void instructionDispatchLoop(OpenMDBitSet& bs); | 
| 139 | > | void instructionDispatchLoop(SelectionSet& bs); | 
| 140 | > | void instructionDispatchLoop(SelectionSet& bs, int frame); | 
| 141 |  |  | 
| 142 | < | void withinInstruction(const Token& instruction, OpenMDBitSet& bs); | 
| 142 | > | void withinInstruction(const Token& instruction, SelectionSet& bs); | 
| 143 | > | void withinInstruction(const Token& instruction, SelectionSet& bs, int frame); | 
| 144 | > | SelectionSet allInstruction(); | 
| 145 |  |  | 
| 146 | < | OpenMDBitSet comparatorInstruction(const Token& instruction); | 
| 147 | < | void compareProperty(StuntDouble* sd, OpenMDBitSet& bs, int property, int comparator, float comparisonValue); | 
| 148 | < | OpenMDBitSet nameInstruction(const std::string& name); | 
| 149 | < | OpenMDBitSet indexInstruction(const boost::any& value); | 
| 150 | < | OpenMDBitSet expression(const std::vector<Token>& tokens, int pc); | 
| 146 | > | SelectionSet comparatorInstruction(const Token& instruction); | 
| 147 | > | SelectionSet comparatorInstruction(const Token& instruction, int frame); | 
| 148 | > | void compareProperty(StuntDouble* sd, SelectionSet& bs, int property, int comparator, float comparisonValue); | 
| 149 | > | void compareProperty(StuntDouble* sd, SelectionSet& bs, int property, int comparator, float comparisonValue, int frame); | 
| 150 | > | SelectionSet nameInstruction(const std::string& name); | 
| 151 | > | SelectionSet indexInstruction(const boost::any& value); | 
| 152 | > | SelectionSet expression(const std::vector<Token>& tokens, int pc); | 
| 153 | > | SelectionSet expression(const std::vector<Token>& tokens, int pc, int frame); | 
| 154 |  |  | 
| 155 | < | OpenMDBitSet lookupValue(const std::string& variable); | 
| 155 | > | SelectionSet lookupValue(const std::string& variable); | 
| 156 |  |  | 
| 157 | < | OpenMDBitSet hull(); | 
| 157 | > | SelectionSet hull(); | 
| 158 | > | SelectionSet hull(int frame); | 
| 159 |  |  | 
| 160 |  | void evalError(const std::string& message) { | 
| 161 |  | sprintf( painCave.errMsg, | 
| 193 |  | bool containDynamicToken(const std::vector<Token>& tokens); | 
| 194 |  |  | 
| 195 |  | RealType getCharge(Atom* atom); | 
| 196 | + | RealType getCharge(Atom* atom, int frame); | 
| 197 |  |  | 
| 198 |  | SelectionCompiler compiler; | 
| 199 |  |  | 
| 207 |  | std::vector<int> linenumbers; | 
| 208 |  | std::vector<int> lineIndices; | 
| 209 |  | std::vector<std::vector<Token> > aatoken; | 
| 210 | < | int pc; // program counter | 
| 210 | > | unsigned int pc; // program counter | 
| 211 |  |  | 
| 212 |  | bool error; | 
| 213 |  | std::string errorMessage; | 
| 220 |  | DistanceFinder distanceFinder; | 
| 221 |  | HullFinder hullFinder; | 
| 222 |  | IndexFinder indexFinder; | 
| 223 | < | int nStuntDouble;   //nGLOBALatoms + nGLOBALrigidbodies | 
| 223 | > | vector<int> nObjects; | 
| 224 |  |  | 
| 225 |  | typedef std::map<std::string, boost::any > VariablesType; | 
| 226 |  | VariablesType variables; | 
| 227 |  |  | 
| 228 |  | bool isDynamic_; | 
| 229 |  | bool isLoaded_; | 
| 230 | + | bool hasSurfaceArea_; | 
| 231 | + | RealType surfaceArea_; | 
| 232 |  |  | 
| 233 |  | }; | 
| 234 |  |  |