| 46 |
|
* @time 13:51am |
| 47 |
|
* @version 1.0 |
| 48 |
|
*/ |
| 49 |
< |
|
| 49 |
> |
#include <exception> |
| 50 |
|
#include <iostream> |
| 51 |
|
#include <sstream> |
| 52 |
|
#include <string> |
| 63 |
|
#include "mdParser/MDParser.hpp" |
| 64 |
|
#include "mdParser/MDTreeParser.hpp" |
| 65 |
|
#include "mdParser/SimplePreprocessor.hpp" |
| 66 |
+ |
#include "antlr/ANTLRException.hpp" |
| 67 |
+ |
#include "antlr/TokenStreamRecognitionException.hpp" |
| 68 |
+ |
#include "antlr/TokenStreamIOException.hpp" |
| 69 |
+ |
#include "antlr/TokenStreamException.hpp" |
| 70 |
+ |
#include "antlr/RecognitionException.hpp" |
| 71 |
+ |
#include "antlr/CharStreamException.hpp" |
| 72 |
|
|
| 73 |
+ |
#include "antlr/MismatchedCharException.hpp" |
| 74 |
+ |
#include "antlr/MismatchedTokenException.hpp" |
| 75 |
+ |
#include "antlr/NoViableAltForCharException.hpp" |
| 76 |
+ |
#include "antlr/NoViableAltException.hpp" |
| 77 |
|
|
| 78 |
|
#ifdef IS_MPI |
| 79 |
|
#include "math/ParallelRandNumGen.hpp" |
| 102 |
|
streamSize = ppStream.str().size() +1; |
| 103 |
|
commStatus = MPI_Bcast(&streamSize, 1, MPI_LONG, masterNode, MPI_COMM_WORLD); |
| 104 |
|
|
| 105 |
< |
commStatus = MPI_Bcast(ppStream.str().c_str(), streamSize, MPI_CHAR, masterNode, MPI_COMM_WORLD); |
| 105 |
> |
commStatus = MPI_Bcast(static_cast<void*>(const_cast<char*>(ppStream.str().c_str())), streamSize, MPI_CHAR, masterNode, MPI_COMM_WORLD); |
| 106 |
|
|
| 107 |
|
|
| 108 |
|
} else { |
| 147 |
|
simParams = treeParser.walkTree(parser.getAST()); |
| 148 |
|
|
| 149 |
|
} |
| 140 |
– |
catch (exception& e) { |
| 141 |
– |
cerr << "parser exception: " << e.what() << endl; |
| 142 |
– |
} |
| 150 |
|
|
| 151 |
+ |
|
| 152 |
+ |
catch(antlr::MismatchedCharException& e) { |
| 153 |
+ |
sprintf(painCave.errMsg, |
| 154 |
+ |
"parser exception: %s %s:%d:%d\n", |
| 155 |
+ |
e.getMessage().c_str(),e.getFilename().c_str(), e.getLine(), e.getColumn()); |
| 156 |
+ |
painCave.isFatal = 1; |
| 157 |
+ |
simError(); |
| 158 |
+ |
} |
| 159 |
+ |
catch(antlr::MismatchedTokenException &e) { |
| 160 |
+ |
sprintf(painCave.errMsg, |
| 161 |
+ |
"parser exception: %s %s:%d:%d\n", |
| 162 |
+ |
e.getMessage().c_str(),e.getFilename().c_str(), e.getLine(), e.getColumn()); |
| 163 |
+ |
painCave.isFatal = 1; |
| 164 |
+ |
simError(); |
| 165 |
+ |
} |
| 166 |
+ |
catch(antlr::NoViableAltForCharException &e) { |
| 167 |
+ |
sprintf(painCave.errMsg, |
| 168 |
+ |
"parser exception: %s %s:%d:%d\n", |
| 169 |
+ |
e.getMessage().c_str(),e.getFilename().c_str(), e.getLine(), e.getColumn()); |
| 170 |
+ |
painCave.isFatal = 1; |
| 171 |
+ |
simError(); |
| 172 |
+ |
} |
| 173 |
+ |
catch(antlr::NoViableAltException &e) { |
| 174 |
+ |
sprintf(painCave.errMsg, |
| 175 |
+ |
"parser exception: %s %s:%d:%d\n", |
| 176 |
+ |
e.getMessage().c_str(),e.getFilename().c_str(), e.getLine(), e.getColumn()); |
| 177 |
+ |
painCave.isFatal = 1; |
| 178 |
+ |
simError(); |
| 179 |
+ |
} |
| 180 |
+ |
|
| 181 |
+ |
catch(antlr::TokenStreamRecognitionException& e) { |
| 182 |
+ |
sprintf(painCave.errMsg, |
| 183 |
+ |
"parser exception: %s %s:%d:%d\n", |
| 184 |
+ |
e.getMessage().c_str(),e.getFilename().c_str(), e.getLine(), e.getColumn()); |
| 185 |
+ |
painCave.isFatal = 1; |
| 186 |
+ |
simError(); |
| 187 |
+ |
} |
| 188 |
+ |
|
| 189 |
+ |
catch(antlr::TokenStreamIOException& e) { |
| 190 |
+ |
sprintf(painCave.errMsg, |
| 191 |
+ |
"parser exception: %s\n", |
| 192 |
+ |
e.getMessage().c_str()); |
| 193 |
+ |
painCave.isFatal = 1; |
| 194 |
+ |
simError(); |
| 195 |
+ |
} |
| 196 |
+ |
|
| 197 |
+ |
catch(antlr::TokenStreamException& e) { |
| 198 |
+ |
sprintf(painCave.errMsg, |
| 199 |
+ |
"parser exception: %s\n", |
| 200 |
+ |
e.getMessage().c_str()); |
| 201 |
+ |
painCave.isFatal = 1; |
| 202 |
+ |
simError(); |
| 203 |
+ |
} |
| 204 |
+ |
catch (antlr::RecognitionException& e) { |
| 205 |
+ |
sprintf(painCave.errMsg, |
| 206 |
+ |
"parser exception: %s %s:%d:%d\n", |
| 207 |
+ |
e.getMessage().c_str(),e.getFilename().c_str(), e.getLine(), e.getColumn()); |
| 208 |
+ |
painCave.isFatal = 1; |
| 209 |
+ |
simError(); |
| 210 |
+ |
} |
| 211 |
+ |
catch (antlr::CharStreamException& e) { |
| 212 |
+ |
sprintf(painCave.errMsg, |
| 213 |
+ |
"parser exception: %s\n", |
| 214 |
+ |
e.getMessage().c_str()); |
| 215 |
+ |
painCave.isFatal = 1; |
| 216 |
+ |
simError(); |
| 217 |
+ |
} |
| 218 |
+ |
catch (OOPSEException& e) { |
| 219 |
+ |
sprintf(painCave.errMsg, |
| 220 |
+ |
"%s\n", |
| 221 |
+ |
e.getMessage().c_str()); |
| 222 |
+ |
painCave.isFatal = 1; |
| 223 |
+ |
simError(); |
| 224 |
+ |
} |
| 225 |
+ |
catch (std::exception& e) { |
| 226 |
+ |
sprintf(painCave.errMsg, |
| 227 |
+ |
"parser exception: %s\n", |
| 228 |
+ |
e.what()); |
| 229 |
+ |
painCave.isFatal = 1; |
| 230 |
+ |
simError(); |
| 231 |
+ |
} |
| 232 |
+ |
|
| 233 |
|
return simParams; |
| 234 |
|
} |
| 235 |
|
|
| 275 |
|
} |
| 276 |
|
|
| 277 |
|
ff->parse(forcefieldFileName); |
| 278 |
< |
|
| 278 |
> |
ff->setFortranForceOptions(); |
| 279 |
|
//create SimInfo |
| 280 |
|
SimInfo * info = new SimInfo(ff, simParams); |
| 281 |
|
|
| 282 |
< |
//gather parameters (SimCreator only retrieves part of the parameters) |
| 282 |
> |
//gather parameters (SimCreator only retrieves part of the |
| 283 |
> |
//parameters) |
| 284 |
|
gatherParameters(info, mdFileName); |
| 285 |
|
|
| 286 |
|
//divide the molecules and determine the global index of molecules |
| 292 |
|
createMolecules(info); |
| 293 |
|
|
| 294 |
|
|
| 295 |
< |
//allocate memory for DataStorage(circular reference, need to break it) |
| 295 |
> |
//allocate memory for DataStorage(circular reference, need to |
| 296 |
> |
//break it) |
| 297 |
|
info->setSnapshotManager(new SimSnapshotManager(info)); |
| 298 |
|
|
| 299 |
< |
//set the global index of atoms, rigidbodies and cutoffgroups (only need to be set once, the |
| 300 |
< |
//global index will never change again). Local indices of atoms and rigidbodies are already set by |
| 301 |
< |
//MoleculeCreator class which actually delegates the responsibility to LocalIndexManager. |
| 299 |
> |
//set the global index of atoms, rigidbodies and cutoffgroups |
| 300 |
> |
//(only need to be set once, the global index will never change |
| 301 |
> |
//again). Local indices of atoms and rigidbodies are already set |
| 302 |
> |
//by MoleculeCreator class which actually delegates the |
| 303 |
> |
//responsibility to LocalIndexManager. |
| 304 |
|
setGlobalIndex(info); |
| 305 |
|
|
| 306 |
< |
//Alought addExculdePairs is called inside SimInfo's addMolecule method, at that point |
| 307 |
< |
//atoms don't have the global index yet (their global index are all initialized to -1). |
| 308 |
< |
//Therefore we have to call addExcludePairs explicitly here. A way to work around is that |
| 309 |
< |
//we can determine the beginning global indices of atoms before they get created. |
| 306 |
> |
//Although addExcludePairs is called inside SimInfo's addMolecule |
| 307 |
> |
//method, at that point atoms don't have the global index yet |
| 308 |
> |
//(their global index are all initialized to -1). Therefore we |
| 309 |
> |
//have to call addExcludePairs explicitly here. A way to work |
| 310 |
> |
//around is that we can determine the beginning global indices of |
| 311 |
> |
//atoms before they get created. |
| 312 |
|
SimInfo::MoleculeIterator mi; |
| 313 |
|
Molecule* mol; |
| 314 |
|
for (mol= info->beginMolecule(mi); mol != NULL; mol = info->nextMolecule(mi)) { |