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

Comparing trunk/src/brains/SimCreator.cpp (file contents):
Revision 770 by tim, Fri Dec 2 15:38:03 2005 UTC vs.
Revision 1051 by gezelter, Mon Sep 25 22:08:33 2006 UTC

# Line 46 | Line 46
46   * @time 13:51am
47   * @version 1.0
48   */
49 <
49 > #include <exception>
50   #include <iostream>
51   #include <sstream>
52   #include <string>
# Line 63 | Line 63
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"
# Line 71 | Line 81 | namespace oopse {
81  
82   namespace oopse {
83    
84 < Globals* SimCreator::parseFile(const std::string mdFileName){
85 <        Globals* simParams = NULL;
86 <        try {
84 >  Globals* SimCreator::parseFile(std::istream& rawMetaDataStream, const std::string& filename, int startOfMetaDataBlock ){
85 >    Globals* simParams = NULL;
86 >    try {
87  
88 <            // Create a preprocessor that preprocesses md file into an ostringstream
89 <            std::stringstream ppStream;
88 >      // Create a preprocessor that preprocesses md file into an ostringstream
89 >      std::stringstream ppStream;
90   #ifdef IS_MPI            
91 <            int streamSize;
92 <            const int masterNode = 0;
93 <            int commStatus;
94 <            if (worldRank == masterNode) {
91 >      int streamSize;
92 >      const int masterNode = 0;
93 >      int commStatus;
94 >      if (worldRank == masterNode) {
95   #endif
96                  
97 <                SimplePreprocessor preprocessor;
98 <                preprocessor.preprocess(mdFileName, ppStream);
97 >        SimplePreprocessor preprocessor;
98 >        preprocessor.preprocess(rawMetaDataStream, filename, startOfMetaDataBlock, ppStream);
99                  
100   #ifdef IS_MPI            
101 <                //brocasting the stream size
102 <                streamSize = ppStream.str().size() +1;
103 <                commStatus = MPI_Bcast(&streamSize, 1, MPI_LONG, masterNode, MPI_COMM_WORLD);                  
101 >        //brocasting the stream size
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 {
109 <                //get stream size
110 <                commStatus = MPI_Bcast(&streamSize, 1, MPI_LONG, masterNode, MPI_COMM_WORLD);  
108 >      } else {
109 >        //get stream size
110 >        commStatus = MPI_Bcast(&streamSize, 1, MPI_LONG, masterNode, MPI_COMM_WORLD);  
111                  
112 <                  char* buf = new char[streamSize];
113 <                  assert(buf);
112 >        char* buf = new char[streamSize];
113 >        assert(buf);
114                  
115 <                  //receive file content
116 <                  commStatus = MPI_Bcast(buf, streamSize, MPI_CHAR, masterNode, MPI_COMM_WORLD);
115 >        //receive file content
116 >        commStatus = MPI_Bcast(buf, streamSize, MPI_CHAR, masterNode, MPI_COMM_WORLD);
117                  
118 <                  ppStream.str(buf);
119 <                  delete buf;
118 >        ppStream.str(buf);
119 >        delete buf;
120  
121 <            }
121 >      }
122   #endif            
123 <            // Create a scanner that reads from the input stream
124 <            MDLexer lexer(ppStream);
125 <            lexer.setFilename(mdFileName);
126 <            lexer.initDeferredLineCount();
123 >      // Create a scanner that reads from the input stream
124 >      MDLexer lexer(ppStream);
125 >      lexer.setFilename(filename);
126 >      lexer.initDeferredLineCount();
127      
128 <            // Create a parser that reads from the scanner
129 <            MDParser parser(lexer);
130 <            parser.setFilename(mdFileName);
128 >      // Create a parser that reads from the scanner
129 >      MDParser parser(lexer);
130 >      parser.setFilename(filename);
131  
132 <            // Create an observer that synchorizes file name change
133 <            FilenameObserver observer;
134 <            observer.setLexer(&lexer);
135 <            observer.setParser(&parser);
136 <            lexer.setObserver(&observer);
132 >      // Create an observer that synchorizes file name change
133 >      FilenameObserver observer;
134 >      observer.setLexer(&lexer);
135 >      observer.setParser(&parser);
136 >      lexer.setObserver(&observer);
137      
138 <            antlr::ASTFactory factory;
139 <            parser.initializeASTFactory(factory);
140 <            parser.setASTFactory(&factory);
141 <            parser.mdfile();
138 >      antlr::ASTFactory factory;
139 >      parser.initializeASTFactory(factory);
140 >      parser.setASTFactory(&factory);
141 >      parser.mdfile();
142  
143 <            // Create a tree parser that reads information into Globals
144 <            MDTreeParser treeParser;
145 <            treeParser.initializeASTFactory(factory);
146 <            treeParser.setASTFactory(&factory);
147 <             simParams = treeParser.walkTree(parser.getAST());
143 >      // Create a tree parser that reads information into Globals
144 >      MDTreeParser treeParser;
145 >      treeParser.initializeASTFactory(factory);
146 >      treeParser.setASTFactory(&factory);
147 >      simParams = treeParser.walkTree(parser.getAST());
148  
149 <        }
140 <        catch (exception& e) {
141 <            cerr << "parser exception: " << e.what() << endl;
142 <        }
149 >    }
150  
151 <        return simParams;
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    
236    SimInfo*  SimCreator::createSim(const std::string & mdFileName,
237                                    bool loadInitCoords) {
238  
239 +    const int bufferSize = 65535;
240 +    char buffer[bufferSize];
241 +    int lineNo = 0;
242 +    std::string mdRawData;
243 +    int metaDataBlockStart = -1;
244 +    int metaDataBlockEnd = -1;
245 +    int i;
246 +    int mdOffset;
247 +
248 + #ifdef IS_MPI            
249 +    const int masterNode = 0;
250 +    if (worldRank == masterNode) {
251 + #endif
252 +
253 +      std::ifstream mdFile_(mdFileName.c_str());
254 +      
255 +      if (mdFile_.fail()) {
256 +        sprintf(painCave.errMsg,
257 +                "SimCreator: Cannot open file: %s\n",
258 +                mdFileName.c_str());
259 +        painCave.isFatal = 1;
260 +        simError();
261 +      }
262 +
263 +      mdFile_.getline(buffer, bufferSize);
264 +      ++lineNo;
265 +      std::string line = trimLeftCopy(buffer);
266 +      i = CaseInsensitiveFind(line, "<OOPSE");
267 +      if (i == string::npos) {
268 +        sprintf(painCave.errMsg,
269 +                "SimCreator: File: %s is not an OOPSE file!\n",
270 +                mdFileName.c_str());
271 +        painCave.isFatal = 1;
272 +        simError();
273 +      }
274 +
275 +      //scan through the input stream and find MetaData tag        
276 +      while(mdFile_.getline(buffer, bufferSize)) {
277 +        ++lineNo;
278 +        
279 +        std::string line = trimLeftCopy(buffer);
280 +        if (metaDataBlockStart == -1) {
281 +          i = CaseInsensitiveFind(line, "<MetaData>");
282 +          if (i != string::npos) {
283 +            metaDataBlockStart = lineNo;
284 +            mdOffset = mdFile_.tellg();
285 +          }
286 +        } else {
287 +          i = CaseInsensitiveFind(line, "</MetaData>");
288 +          if (i != string::npos) {
289 +            metaDataBlockEnd = lineNo;
290 +          }
291 +        }
292 +      }
293 +
294 +      if (metaDataBlockStart == -1) {
295 +        sprintf(painCave.errMsg,
296 +                "SimCreator: File: %s did not contain a <MetaData> tag!\n",
297 +                mdFileName.c_str());
298 +        painCave.isFatal = 1;
299 +        simError();
300 +      }
301 +      if (metaDataBlockEnd == -1) {
302 +        sprintf(painCave.errMsg,
303 +                "SimCreator: File: %s did not contain a closed MetaData block!\n",
304 +                mdFileName.c_str());
305 +        painCave.isFatal = 1;
306 +        simError();
307 +      }
308 +        
309 +      mdFile_.clear();
310 +      mdFile_.seekg(0);
311 +      mdFile_.seekg(mdOffset);
312 +
313 +      mdRawData.clear();
314 +
315 +      for (int i = 0; i < metaDataBlockEnd - metaDataBlockStart - 1; ++i) {
316 +        mdFile_.getline(buffer, bufferSize);
317 +        mdRawData += buffer;
318 +        mdRawData += "\n";
319 +      }
320 +
321 +      mdFile_.close();
322 +
323 + #ifdef IS_MPI
324 +    }
325 + #endif
326 +
327 +    std::stringstream rawMetaDataStream(mdRawData);
328 +
329      //parse meta-data file
330 <    Globals* simParams = parseFile(mdFileName);
330 >    Globals* simParams = parseFile(rawMetaDataStream, mdFileName, metaDataBlockStart+1);
331      
332      //create the force field
333      ForceField * ff = ForceFieldFactory::getInstance()
# Line 163 | Line 342 | Globals* SimCreator::parseFile(const std::string mdFil
342      }
343      
344      if (simParams->haveForceFieldFileName()) {
345 +      std::cout<< simParams->getForceFieldFileName() << "\n";
346        ff->setForceFieldFileName(simParams->getForceFieldFileName());
347      }
348      
# Line 186 | Line 366 | Globals* SimCreator::parseFile(const std::string mdFil
366      }
367      
368      ff->parse(forcefieldFileName);
369 <        
369 >    ff->setFortranForceOptions();
370      //create SimInfo
371      SimInfo * info = new SimInfo(ff, simParams);
372 +
373 +    info->setRawMetaData(mdRawData);
374      
375 <    //gather parameters (SimCreator only retrieves part of the parameters)
375 >    //gather parameters (SimCreator only retrieves part of the
376 >    //parameters)
377      gatherParameters(info, mdFileName);
378      
379      //divide the molecules and determine the global index of molecules
# Line 202 | Line 385 | Globals* SimCreator::parseFile(const std::string mdFil
385      createMolecules(info);
386      
387      
388 <    //allocate memory for DataStorage(circular reference, need to break it)
388 >    //allocate memory for DataStorage(circular reference, need to
389 >    //break it)
390      info->setSnapshotManager(new SimSnapshotManager(info));
391      
392 <    //set the global index of atoms, rigidbodies and cutoffgroups (only need to be set once, the
393 <    //global index will never change again). Local indices of atoms and rigidbodies are already set by
394 <    //MoleculeCreator class which actually delegates the responsibility to LocalIndexManager.
392 >    //set the global index of atoms, rigidbodies and cutoffgroups
393 >    //(only need to be set once, the global index will never change
394 >    //again). Local indices of atoms and rigidbodies are already set
395 >    //by MoleculeCreator class which actually delegates the
396 >    //responsibility to LocalIndexManager.
397      setGlobalIndex(info);
398      
399 <    //Alought addExculdePairs is called inside SimInfo's addMolecule method, at that point
400 <    //atoms don't have the global index yet  (their global index are all initialized to -1).
401 <    //Therefore we have to call addExcludePairs explicitly here. A way to work around is that
402 <    //we can determine the beginning global indices of atoms before they get created.
399 >    //Although addExcludePairs is called inside SimInfo's addMolecule
400 >    //method, at that point atoms don't have the global index yet
401 >    //(their global index are all initialized to -1).  Therefore we
402 >    //have to call addExcludePairs explicitly here. A way to work
403 >    //around is that we can determine the beginning global indices of
404 >    //atoms before they get created.
405      SimInfo::MoleculeIterator mi;
406      Molecule* mol;
407      for (mol= info->beginMolecule(mi); mol != NULL; mol = info->nextMolecule(mi)) {
# Line 221 | Line 409 | Globals* SimCreator::parseFile(const std::string mdFil
409      }
410      
411      if (loadInitCoords)
412 <      loadCoordinates(info);    
412 >      loadCoordinates(info, mdFileName);    
413      
414      return info;
415    }
# Line 257 | Line 445 | Globals* SimCreator::parseFile(const std::string mdFil
445    
446   #ifdef IS_MPI
447    void SimCreator::divideMolecules(SimInfo *info) {
448 <    double numerator;
449 <    double denominator;
450 <    double precast;
451 <    double x;
452 <    double y;
453 <    double a;
448 >    RealType numerator;
449 >    RealType denominator;
450 >    RealType precast;
451 >    RealType x;
452 >    RealType y;
453 >    RealType a;
454      int old_atoms;
455      int add_atoms;
456      int new_atoms;
# Line 374 | Line 562 | Globals* SimCreator::parseFile(const std::string mdFil
562            //           Pacc(x) = exp(- a * x)
563            // where a = penalty / (average atoms per molecule)
564            
565 <          x = (double)(new_atoms - nTarget);
565 >          x = (RealType)(new_atoms - nTarget);
566            y = myRandom->rand();
567            
568            if (y < exp(- a * x)) {
# Line 440 | Line 628 | Globals* SimCreator::parseFile(const std::string mdFil
628      Molecule::AtomIterator ai;
629      Molecule::RigidBodyIterator ri;
630      Molecule::CutoffGroupIterator ci;
631 +    Molecule::IntegrableObjectIterator  ioi;
632      Molecule * mol;
633      Atom * atom;
634      RigidBody * rb;
# Line 448 | Line 637 | Globals* SimCreator::parseFile(const std::string mdFil
637      int beginRigidBodyIndex;
638      int beginCutoffGroupIndex;
639      int nGlobalAtoms = info->getNGlobalAtoms();
640 <    
640 >
641 >    /**@todo fixme */
642   #ifndef IS_MPI
643      
644      beginAtomIndex = 0;
# Line 563 | Line 753 | Globals* SimCreator::parseFile(const std::string mdFil
753   #else
754      info->setGlobalMolMembership(globalMolMembership);
755   #endif
756 +
757 +    // nIOPerMol holds the number of integrable objects per molecule
758 +    // here the molecules are listed by their global indices.
759 +
760 +    std::vector<int> nIOPerMol(info->getNGlobalMolecules(), 0);
761 +    for (mol = info->beginMolecule(mi); mol != NULL; mol = info->nextMolecule(mi)) {
762 +      nIOPerMol[mol->getGlobalIndex()] = mol->getNIntegrableObjects();      
763 +    }
764      
765 + #ifdef IS_MPI
766 +    std::vector<int> numIntegrableObjectsPerMol(info->getNGlobalMolecules(), 0);
767 +    MPI_Allreduce(&nIOPerMol[0], &numIntegrableObjectsPerMol[0],
768 +                  info->getNGlobalMolecules(), MPI_INT, MPI_SUM, MPI_COMM_WORLD);
769 + #else
770 +    std::vector<int> numIntegrableObjectsPerMol = nIOPerMol;
771 + #endif    
772 +
773 + std::vector<int> startingIOIndexForMol(info->getNGlobalMolecules());
774 +
775 + int startingIndex = 0;
776 + for (int i = 0; i < info->getNGlobalMolecules(); i++) {
777 +  startingIOIndexForMol[i] = startingIndex;
778 +  startingIndex += numIntegrableObjectsPerMol[i];
779 + }
780 +
781 + std::vector<StuntDouble*> IOIndexToIntegrableObject(info->getNGlobalIntegrableObjects(), (StuntDouble*)NULL);
782 + for (mol = info->beginMolecule(mi); mol != NULL; mol = info->nextMolecule(mi)) {
783 +      int myGlobalIndex = mol->getGlobalIndex();
784 +      int globalIO = startingIOIndexForMol[myGlobalIndex];
785 +      for (StuntDouble* integrableObject = mol->beginIntegrableObject(ioi); integrableObject != NULL;
786 +           integrableObject = mol->nextIntegrableObject(ioi)) {
787 +            integrableObject->setGlobalIntegrableObjectIndex(globalIO);
788 +            IOIndexToIntegrableObject[globalIO] = integrableObject;
789 +            globalIO++;
790 +      }
791 +    }
792 +
793 +  info->setIOIndexToIntegrableObject(IOIndexToIntegrableObject);
794 +  
795    }
796    
797 <  void SimCreator::loadCoordinates(SimInfo* info) {
797 >  void SimCreator::loadCoordinates(SimInfo* info, const std::string& mdFileName) {
798      Globals* simParams;
799      simParams = info->getSimParams();
800      
573 –    if (!simParams->haveInitialConfig()) {
574 –      sprintf(painCave.errMsg,
575 –              "Cannot intialize a simulation without an initial configuration file.\n");
576 –      painCave.isFatal = 1;;
577 –      simError();
578 –    }
801      
802 <    DumpReader reader(info, simParams->getInitialConfig());
802 >    DumpReader reader(info, mdFileName);
803      int nframes = reader.getNFrames();
804      
805      if (nframes > 0) {
# Line 586 | Line 808 | Globals* SimCreator::parseFile(const std::string mdFil
808        //invalid initial coordinate file
809        sprintf(painCave.errMsg,
810                "Initial configuration file %s should at least contain one frame\n",
811 <              simParams->getInitialConfig().c_str());
811 >              mdFileName.c_str());
812        painCave.isFatal = 1;
813        simError();
814      }

Diff Legend

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