src/io/DumpReader.cpp

/* 
 * Copyright (c) 2005 The University of Notre Dame. All Rights Reserved. 
 * 
 * The University of Notre Dame grants you ("Licensee") a 
 * non-exclusive, royalty free, license to use, modify and 
 * redistribute this software in source and binary code form, provided 
 * that the following conditions are met: 
 * 
 * 1. Acknowledgement of the program authors must be made in any 
 *    publication of scientific results based in part on use of the 
 *    program.  An acceptable form of acknowledgement is citation of 
 *    the article in which the program was described (Matthew 
 *    A. Meineke, Charles F. Vardeman II, Teng Lin, Christopher 
 *    J. Fennell and J. Daniel Gezelter, "OOPSE: An Object-Oriented 
 *    Parallel Simulation Engine for Molecular Dynamics," 
 *    J. Comput. Chem. 26, pp. 252-271 (2005)) 
 * 
 * 2. Redistributions of source code must retain the above copyright 
 *    notice, this list of conditions and the following disclaimer. 
 * 
 * 3. Redistributions in binary form must reproduce the above copyright 
 *    notice, this list of conditions and the following disclaimer in the 
 *    documentation and/or other materials provided with the 
 *    distribution. 
 * 
 * This software is provided "AS IS," without a warranty of any 
 * kind. All express or implied conditions, representations and 
 * warranties, including any implied warranty of merchantability, 
 * fitness for a particular purpose or non-infringement, are hereby 
 * excluded.  The University of Notre Dame and its licensors shall not 
 * be liable for any damages suffered by licensee as a result of 
 * using, modifying or distributing the software or its 
 * derivatives. In no event will the University of Notre Dame or its 
 * licensors be liable for any lost revenue, profit or data, or for 
 * direct, indirect, special, consequential, incidental or punitive 
 * damages, however caused and regardless of the theory of liability, 
 * arising out of the use of or inability to use software, even if the 
 * University of Notre Dame has been advised of the possibility of 
 * such damages. 
 */ 
  
#define _LARGEFILE_SOURCE64 
#define _FILE_OFFSET_BITS 64 
 
#include <sys/types.h> 
#include <sys/stat.h> 
 
#include <iostream> 
#include <math.h> 
 
#include <stdio.h> 
#include <stdlib.h> 
#include <string.h> 
 
#include "io/DumpReader.hpp" 
#include "primitives/Molecule.hpp" 
#include "utils/simError.h" 
#include "utils/MemoryUtils.hpp" 
#include "utils/StringTokenizer.hpp" 
 
#ifdef IS_MPI 
 
#include <mpi.h> 
#define TAKE_THIS_TAG_CHAR 0 
#define TAKE_THIS_TAG_INT 1 
 
#endif // is_mpi 
 
 
namespace oopse { 
   
  DumpReader::DumpReader(SimInfo* info, const std::string& filename) 
    : info_(info), filename_(filename), isScanned_(false), nframes_(0), needCOMprops_(false) { 
    
#ifdef IS_MPI 
    
    if (worldRank == 0) { 
#endif 
      
      inFile_ = new std::ifstream(filename_.c_str()); 
      
      if (inFile_->fail()) { 
        sprintf(painCave.errMsg, 
                "DumpReader: Cannot open file: %s\n", 
                filename_.c_str()); 
        painCave.isFatal = 1; 
        simError(); 
      } 
      
#ifdef IS_MPI 
      
    } 
    
    strcpy(checkPointMsg, "Dump file opened for reading successfully."); 
    errorCheckPoint(); 
    
#endif 
    
    return; 
  } 
  
  DumpReader::~DumpReader() { 
    
#ifdef IS_MPI 
    
    if (worldRank == 0) { 
#endif 
      
      delete inFile_; 
      
#ifdef IS_MPI 
      
    } 
    
    strcpy(checkPointMsg, "Dump file closed successfully."); 
    errorCheckPoint(); 
    
#endif 
    
    return; 
  } 
  
  int DumpReader::getNFrames(void) { 
     
    if (!isScanned_) 
      scanFile(); 
     
    return nframes_; 
  } 
   
  void DumpReader::scanFile(void) { 
    int lineNo = 0; 
    std::streampos prevPos;
    std::streampos  currPos; 
    
#ifdef IS_MPI 
    
    if (worldRank == 0) { 
#endif // is_mpi 
      
      currPos = inFile_->tellg();
      prevPos = currPos;
      bool foundOpenSnapshotTag = false;
      bool foundClosedSnapshotTag = false;
      while(inFile_->getline(buffer, bufferSize)) {
        ++lineNo;
        
        std::string line = buffer;
        currPos = inFile_->tellg(); 
        if (line.find("<Snapshot>")!= std::string::npos) {
          if (foundOpenSnapshotTag) {
            sprintf(painCave.errMsg, 
                    "DumpReader:<Snapshot> is multiply nested at line %d in %s \n", lineNo, 
                    filename_.c_str()); 
            painCave.isFatal = 1; 
            simError();           
          }
          foundOpenSnapshotTag = true;
          foundClosedSnapshotTag = false;
          framePos_.push_back(prevPos);
          
        } else if (line.find("</Snapshot>") != std::string::npos){
          if (!foundOpenSnapshotTag) {
            sprintf(painCave.errMsg, 
                    "DumpReader:</Snapshot> appears before <Snapshot> at line %d in %s \n", lineNo, 
                    filename_.c_str()); 
            painCave.isFatal = 1; 
            simError(); 
          }
          
          if (foundClosedSnapshotTag) {
            sprintf(painCave.errMsg, 
                    "DumpReader:</Snapshot> appears multiply nested at line %d in %s \n", lineNo, 
                    filename_.c_str()); 
            painCave.isFatal = 1; 
            simError(); 
          }
          foundClosedSnapshotTag = true;
          foundOpenSnapshotTag = false;
        }
        prevPos = currPos;
      }
      
      // only found <Snapshot> for the last frame means the file is corrupted, we should discard
      // it and give a warning message
      if (foundOpenSnapshotTag) {
        sprintf(painCave.errMsg, 
                "DumpReader: last frame in %s is invalid\n", filename_.c_str()); 
        painCave.isFatal = 0; 
        simError();       
        framePos_.pop_back();
      }
      
      nframes_ = framePos_.size(); 
      
      if (nframes_ == 0) {
        sprintf(painCave.errMsg, 
                "DumpReader: %s does not contain a valid frame\n", filename_.c_str()); 
        painCave.isFatal = 1; 
        simError();      
      }
#ifdef IS_MPI 
    } 
     
    MPI_Bcast(&nframes_, 1, MPI_INT, 0, MPI_COMM_WORLD); 
    
#endif // is_mpi 
    
    isScanned_ = true; 
  } 
   
  void DumpReader::readFrame(int whichFrame) { 
    if (!isScanned_) 
      scanFile(); 
         
    int storageLayout = info_->getSnapshotManager()->getStorageLayout(); 
     
    if (storageLayout & DataStorage::dslPosition) { 
      needPos_ = true; 
    } else { 
      needPos_ = false; 
    } 
     
    if (storageLayout & DataStorage::dslVelocity) { 
      needVel_ = true; 
    } else { 
      needVel_ = false; 
    } 
     
    if (storageLayout & DataStorage::dslAmat || storageLayout & DataStorage::dslElectroFrame) { 
      needQuaternion_ = true; 
    } else { 
      needQuaternion_ = false; 
    } 
     
    if (storageLayout & DataStorage::dslAngularMomentum) { 
      needAngMom_ = true; 
    } else { 
      needAngMom_ = false;     
    } 
     
    readSet(whichFrame); 

    if (needCOMprops_) {
      Snapshot* s = info_->getSnapshotManager()->getCurrentSnapshot();
      Vector3d com;
      Vector3d comvel;
      Vector3d comw;
      if (needPos_ && needVel_){
        info_->getComAll(com, comvel);
        comw = info_->getAngularMomentum();
      }else{
        com = info_->getCom();
        comvel = 0.0;
        comw   = 0.0;
      }
      s->setCOMprops(com, comvel, comw);      
    }

  } 
   
  void DumpReader::readSet(int whichFrame) {     
    std::string line;

#ifndef IS_MPI 

    inFile_->clear();  
    inFile_->seekg(framePos_[whichFrame]); 

    std::istream& inputStream = *inFile_;     

#else 
    int masterNode = 0;
    std::stringstream sstream;
    if (worldRank == masterNode) {
      std::string sendBuffer;

      inFile_->clear();  
      inFile_->seekg(framePos_[whichFrame]); 
      
      while (inFile_->getline(buffer, bufferSize)) {

        line = buffer;
        sendBuffer += line;
        sendBuffer += '\n';
        if (line.find("</Snapshot>") != std::string::npos) {
          break;
        }        
      }

      int sendBufferSize = sendBuffer.size();
      MPI_Bcast(&sendBufferSize, 1, MPI_INT, masterNode, MPI_COMM_WORLD);     
      MPI_Bcast((void *)sendBuffer.c_str(), sendBufferSize, MPI_CHAR, masterNode, MPI_COMM_WORLD);     
      
      sstream.str(sendBuffer);
    } else {
      int sendBufferSize;
      MPI_Bcast(&sendBufferSize, 1, MPI_INT, masterNode, MPI_COMM_WORLD);     
      char * recvBuffer = new char[sendBufferSize+1];
      MPI_Bcast(recvBuffer, sendBufferSize, MPI_CHAR, masterNode, MPI_COMM_WORLD);     
      sstream.str(recvBuffer);
    }      

    std::istream& inputStream = sstream;  
#endif

    inputStream.getline(buffer, bufferSize);

    line = buffer;
    if (line.find("<Snapshot>") == std::string::npos) {
      sprintf(painCave.errMsg, 
              "DumpReader Error: can not find <Snapshot>\n"); 
      painCave.isFatal = 1; 
      simError(); 
    } 
    
    //read frameData
    readFrameProperties(inputStream);

    //read StuntDoubles
    readStuntDoubles(inputStream);     

    inputStream.getline(buffer, bufferSize);
    line = buffer;
    if (line.find("</Snapshot>") == std::string::npos) {
      sprintf(painCave.errMsg, 
              "DumpReader Error: can not find </Snapshot>\n"); 
      painCave.isFatal = 1; 
      simError(); 
    }        
   
  } 
   
  void DumpReader::parseDumpLine(const std::string& line) { 

       
    StringTokenizer tokenizer(line); 
    int nTokens; 
     
    nTokens = tokenizer.countTokens(); 
     
    if (nTokens < 2) { 
      sprintf(painCave.errMsg, 
              "DumpReader Error: Not enough Tokens.\n%s\n", line.c_str()); 
      painCave.isFatal = 1; 
      simError(); 
    } 

    int index = tokenizer.nextTokenAsInt();
 
    StuntDouble* integrableObject = info_->getIOIndexToIntegrableObject(index);

    if (integrableObject == NULL) {
      return;
    }
    std::string type = tokenizer.nextToken(); 
    int size = type.size();

    for(int i = 0; i < size; ++i) {
      switch(type[i]) {
        
        case 'p': {
            Vector3d pos;
            pos[0] = tokenizer.nextTokenAsDouble(); 
            pos[1] = tokenizer.nextTokenAsDouble(); 
            pos[2] = tokenizer.nextTokenAsDouble(); 
            if (needPos_) { 
              integrableObject->setPos(pos); 
            }             
            break;
        }
        case 'v' : {
            Vector3d vel;
            vel[0] = tokenizer.nextTokenAsDouble(); 
            vel[1] = tokenizer.nextTokenAsDouble(); 
            vel[2] = tokenizer.nextTokenAsDouble(); 
            if (needVel_) { 
              integrableObject->setVel(vel); 
            } 
            break;
        }

        case 'q' : {
           Quat4d q;
           if (integrableObject->isDirectional()) { 
              
             q[0] = tokenizer.nextTokenAsDouble(); 
             q[1] = tokenizer.nextTokenAsDouble(); 
             q[2] = tokenizer.nextTokenAsDouble(); 
             q[3] = tokenizer.nextTokenAsDouble(); 
              
             RealType qlen = q.length(); 
             if (qlen < oopse::epsilon) { //check quaternion is not equal to 0 
                
               sprintf(painCave.errMsg, 
                       "DumpReader Error: initial quaternion error (q0^2 + q1^2 + q2^2 + q3^2) ~ 0\n"); 
               painCave.isFatal = 1; 
               simError(); 
                
             }  
              
             q.normalize(); 
             if (needQuaternion_) {            
               integrableObject->setQ(q); 
             }               
           }            
           break;
        }  
        case 'j' : {
          Vector3d ji;
          if (integrableObject->isDirectional()) {
             ji[0] = tokenizer.nextTokenAsDouble(); 
             ji[1] = tokenizer.nextTokenAsDouble(); 
             ji[2] = tokenizer.nextTokenAsDouble(); 
             if (needAngMom_) { 
               integrableObject->setJ(ji); 
             } 
          }
          break;
        }  
        case 'f': {

          Vector3d force;
          force[0] = tokenizer.nextTokenAsDouble(); 
          force[1] = tokenizer.nextTokenAsDouble(); 
          force[2] = tokenizer.nextTokenAsDouble();           
          integrableObject->setFrc(force); 
          break;
        }
        case 't' : {

           Vector3d torque;
           torque[0] = tokenizer.nextTokenAsDouble(); 
           torque[1] = tokenizer.nextTokenAsDouble(); 
           torque[2] = tokenizer.nextTokenAsDouble();           
           integrableObject->setTrq(torque);          
           break;
        }
        default: {
               sprintf(painCave.errMsg, 
                       "DumpReader Error: %s is an unrecognized type\n", type.c_str()); 
               painCave.isFatal = 1; 
               simError(); 
          break;   
        }

      }
    }
     
  } 
   

  void  DumpReader::readStuntDoubles(std::istream& inputStream) {

    inputStream.getline(buffer, bufferSize);
    std::string line(buffer);
    
    if (line.find("<StuntDoubles>") == std::string::npos) {
      sprintf(painCave.errMsg, 
              "DumpReader Error: Missing <StuntDoubles>\n"); 
      painCave.isFatal = 1; 
      simError(); 
    }

    while(inputStream.getline(buffer, bufferSize)) {
      line = buffer;
      
      if(line.find("</StuntDoubles>") != std::string::npos) {
        break;
      }

      parseDumpLine(line);
    }
  
  }

  void DumpReader::readFrameProperties(std::istream& inputStream) {

    Snapshot* s = info_->getSnapshotManager()->getCurrentSnapshot();
    inputStream.getline(buffer, bufferSize);
    std::string line(buffer);

    if (line.find("<FrameData>") == std::string::npos) {
      sprintf(painCave.errMsg, 
              "DumpReader Error: Missing <FrameData>\n"); 
      painCave.isFatal = 1; 
      simError(); 
    }

    while(inputStream.getline(buffer, bufferSize)) {
      line = buffer;
      
      if(line.find("</FrameData>") != std::string::npos) {
        break;
      }
      
      StringTokenizer tokenizer(line, " ;\t\n\r{}:,");
      if (!tokenizer.hasMoreTokens()) {
        sprintf(painCave.errMsg, 
                "DumpReader Error: Not enough Tokens.\n%s\n", line.c_str()); 
        painCave.isFatal = 1; 
        simError();      
      }

      std::string propertyName = tokenizer.nextToken();
      if (propertyName == "Time") {
        RealType currTime = tokenizer.nextTokenAsDouble(); 
        s->setTime(currTime); 
      } else if (propertyName == "Hmat"){
        Mat3x3d hmat;
        hmat(0, 0) = tokenizer.nextTokenAsDouble(); 
        hmat(0, 1) = tokenizer.nextTokenAsDouble(); 
        hmat(0, 2) = tokenizer.nextTokenAsDouble(); 
        hmat(1, 0) = tokenizer.nextTokenAsDouble(); 
        hmat(1, 1) = tokenizer.nextTokenAsDouble(); 
        hmat(1, 2) = tokenizer.nextTokenAsDouble(); 
        hmat(2, 0) = tokenizer.nextTokenAsDouble(); 
        hmat(2, 1) = tokenizer.nextTokenAsDouble(); 
        hmat(2, 2) = tokenizer.nextTokenAsDouble(); 
        s->setHmat(hmat);      
      } else if (propertyName == "Thermostat") {
        RealType chi = tokenizer.nextTokenAsDouble();
        RealType integralOfChiDt = tokenizer.nextTokenAsDouble();
        s->setChi(chi); 
        s->setIntegralOfChiDt(integralOfChiDt);        
     } else if (propertyName == "Barostat") {
        Mat3x3d eta;
        eta(0, 0) = tokenizer.nextTokenAsDouble(); 
        eta(0, 1) = tokenizer.nextTokenAsDouble(); 
        eta(0, 2) = tokenizer.nextTokenAsDouble(); 
        eta(1, 0) = tokenizer.nextTokenAsDouble(); 
        eta(1, 1) = tokenizer.nextTokenAsDouble(); 
        eta(1, 2) = tokenizer.nextTokenAsDouble(); 
        eta(2, 0) = tokenizer.nextTokenAsDouble(); 
        eta(2, 1) = tokenizer.nextTokenAsDouble(); 
        eta(2, 2) = tokenizer.nextTokenAsDouble(); 
        s->setEta(eta); 
      } else {
        sprintf(painCave.errMsg, 
                "DumpReader Error: %s is an invalid property in <FrameData>\n", propertyName.c_str()); 
        painCave.isFatal = 0; 
        simError();         
      }
      
    }

  }

   
}//end namespace oopse 
Revision:	1241
Committed:	Fri Apr 25 15:14:47 2008 UTC (17 years ago) by gezelter
File size:	16223 byte(s)
Log Message:	A bunch of minor changes to make MPI compilation faster than the double compilation we do now...
#	User	Rev	Content
1	chrisfen	721	/*
2			* Copyright (c) 2005 The University of Notre Dame. All Rights Reserved.
3			*
4			* The University of Notre Dame grants you ("Licensee") a
5			* non-exclusive, royalty free, license to use, modify and
6			* redistribute this software in source and binary code form, provided
7			* that the following conditions are met:
8			*
9			* 1. Acknowledgement of the program authors must be made in any
10			* publication of scientific results based in part on use of the
11			* program. An acceptable form of acknowledgement is citation of
12			* the article in which the program was described (Matthew
13			* A. Meineke, Charles F. Vardeman II, Teng Lin, Christopher
14			* J. Fennell and J. Daniel Gezelter, "OOPSE: An Object-Oriented
15			* Parallel Simulation Engine for Molecular Dynamics,"
16			* J. Comput. Chem. 26, pp. 252-271 (2005))
17			*
18			* 2. Redistributions of source code must retain the above copyright
19			* notice, this list of conditions and the following disclaimer.
20			*
21			* 3. Redistributions in binary form must reproduce the above copyright
22			* notice, this list of conditions and the following disclaimer in the
23			* documentation and/or other materials provided with the
24			* distribution.
25			*
26			* This software is provided "AS IS," without a warranty of any
27			* kind. All express or implied conditions, representations and
28			* warranties, including any implied warranty of merchantability,
29			* fitness for a particular purpose or non-infringement, are hereby
30			* excluded. The University of Notre Dame and its licensors shall not
31			* be liable for any damages suffered by licensee as a result of
32			* using, modifying or distributing the software or its
33			* derivatives. In no event will the University of Notre Dame or its
34			* licensors be liable for any lost revenue, profit or data, or for
35			* direct, indirect, special, consequential, incidental or punitive
36			* damages, however caused and regardless of the theory of liability,
37			* arising out of the use of or inability to use software, even if the
38			* University of Notre Dame has been advised of the possibility of
39			* such damages.
40			*/
41
42			#define _LARGEFILE_SOURCE64
43			#define _FILE_OFFSET_BITS 64
44
45			#include <sys/types.h>
46			#include <sys/stat.h>
47
48			#include <iostream>
49			#include <math.h>
50
51			#include <stdio.h>
52			#include <stdlib.h>
53			#include <string.h>
54
55			#include "io/DumpReader.hpp"
56			#include "primitives/Molecule.hpp"
57			#include "utils/simError.h"
58			#include "utils/MemoryUtils.hpp"
59			#include "utils/StringTokenizer.hpp"
60
61			#ifdef IS_MPI
62
63			#include <mpi.h>
64			#define TAKE_THIS_TAG_CHAR 0
65			#define TAKE_THIS_TAG_INT 1
66
67			#endif // is_mpi
68
69
70			namespace oopse {
71
72			DumpReader::DumpReader(SimInfo* info, const std::string& filename)
73	gezelter	1104	: info_(info), filename_(filename), isScanned_(false), nframes_(0), needCOMprops_(false) {
74	chrisfen	996
75	chrisfen	721	#ifdef IS_MPI
76	chrisfen	996
77			if (worldRank == 0) {
78	chrisfen	721	#endif
79	chrisfen	996
80	chrisfen	721	inFile_ = new std::ifstream(filename_.c_str());
81	chrisfen	996
82			if (inFile_->fail()) {
83			sprintf(painCave.errMsg,
84			"DumpReader: Cannot open file: %s\n",
85			filename_.c_str());
86			painCave.isFatal = 1;
87			simError();
88			}
89
90	chrisfen	721	#ifdef IS_MPI
91	chrisfen	996
92			}
93
94			strcpy(checkPointMsg, "Dump file opened for reading successfully.");
95	gezelter	1241	errorCheckPoint();
96	chrisfen	996
97	chrisfen	721	#endif
98	chrisfen	996
99			return;
100			}
101
102	chrisfen	721	DumpReader::~DumpReader() {
103	chrisfen	996
104	chrisfen	721	#ifdef IS_MPI
105	chrisfen	996
106	chrisfen	721	if (worldRank == 0) {
107			#endif
108	chrisfen	996
109	chrisfen	721	delete inFile_;
110	chrisfen	996
111	chrisfen	721	#ifdef IS_MPI
112	chrisfen	996
113	chrisfen	721	}
114	chrisfen	996
115	chrisfen	721	strcpy(checkPointMsg, "Dump file closed successfully.");
116	gezelter	1241	errorCheckPoint();
117	chrisfen	996
118	chrisfen	721	#endif
119	chrisfen	996
120	chrisfen	721	return;
121			}
122	chrisfen	996
123	chrisfen	721	int DumpReader::getNFrames(void) {
124
125			if (!isScanned_)
126			scanFile();
127
128			return nframes_;
129			}
130
131			void DumpReader::scanFile(void) {
132	tim	1024	int lineNo = 0;
133			std::streampos prevPos;
134	chrisfen	721	std::streampos currPos;
135	tim	1024
136	chrisfen	721	#ifdef IS_MPI
137	tim	1024
138	chrisfen	721	if (worldRank == 0) {
139			#endif // is_mpi
140	tim	1024
141			currPos = inFile_->tellg();
142			prevPos = currPos;
143			bool foundOpenSnapshotTag = false;
144			bool foundClosedSnapshotTag = false;
145			while(inFile_->getline(buffer, bufferSize)) {
146			++lineNo;
147
148			std::string line = buffer;
149			currPos = inFile_->tellg();
150			if (line.find("<Snapshot>")!= std::string::npos) {
151			if (foundOpenSnapshotTag) {
152	chrisfen	721	sprintf(painCave.errMsg,
153	tim	1024	"DumpReader:<Snapshot> is multiply nested at line %d in %s \n", lineNo,
154			filename_.c_str());
155	chrisfen	721	painCave.isFatal = 1;
156	tim	1024	simError();
157			}
158			foundOpenSnapshotTag = true;
159			foundClosedSnapshotTag = false;
160			framePos_.push_back(prevPos);
161
162			} else if (line.find("</Snapshot>") != std::string::npos){
163			if (!foundOpenSnapshotTag) {
164			sprintf(painCave.errMsg,
165			"DumpReader:</Snapshot> appears before <Snapshot> at line %d in %s \n", lineNo,
166			filename_.c_str());
167			painCave.isFatal = 1;
168	chrisfen	721	simError();
169	tim	1024	}
170
171			if (foundClosedSnapshotTag) {
172			sprintf(painCave.errMsg,
173			"DumpReader:</Snapshot> appears multiply nested at line %d in %s \n", lineNo,
174			filename_.c_str());
175			painCave.isFatal = 1;
176			simError();
177			}
178			foundClosedSnapshotTag = true;
179			foundOpenSnapshotTag = false;
180			}
181			prevPos = currPos;
182			}
183
184			// only found <Snapshot> for the last frame means the file is corrupted, we should discard
185			// it and give a warning message
186			if (foundOpenSnapshotTag) {
187			sprintf(painCave.errMsg,
188			"DumpReader: last frame in %s is invalid\n", filename_.c_str());
189			painCave.isFatal = 0;
190			simError();
191			framePos_.pop_back();
192			}
193
194	chrisfen	721	nframes_ = framePos_.size();
195	tim	1024
196			if (nframes_ == 0) {
197			sprintf(painCave.errMsg,
198			"DumpReader: %s does not contain a valid frame\n", filename_.c_str());
199			painCave.isFatal = 1;
200			simError();
201			}
202	chrisfen	721	#ifdef IS_MPI
203			}
204
205			MPI_Bcast(&nframes_, 1, MPI_INT, 0, MPI_COMM_WORLD);
206	tim	1024
207	chrisfen	721	#endif // is_mpi
208	tim	1024
209	chrisfen	721	isScanned_ = true;
210			}
211
212			void DumpReader::readFrame(int whichFrame) {
213			if (!isScanned_)
214			scanFile();
215
216			int storageLayout = info_->getSnapshotManager()->getStorageLayout();
217
218			if (storageLayout & DataStorage::dslPosition) {
219			needPos_ = true;
220			} else {
221			needPos_ = false;
222			}
223
224			if (storageLayout & DataStorage::dslVelocity) {
225			needVel_ = true;
226			} else {
227			needVel_ = false;
228			}
229
230			if (storageLayout & DataStorage::dslAmat \|\| storageLayout & DataStorage::dslElectroFrame) {
231			needQuaternion_ = true;
232			} else {
233			needQuaternion_ = false;
234			}
235
236			if (storageLayout & DataStorage::dslAngularMomentum) {
237			needAngMom_ = true;
238			} else {
239			needAngMom_ = false;
240			}
241
242			readSet(whichFrame);
243	gezelter	1104
244			if (needCOMprops_) {
245			Snapshot* s = info_->getSnapshotManager()->getCurrentSnapshot();
246			Vector3d com;
247			Vector3d comvel;
248			Vector3d comw;
249	chuckv	1112	if (needPos_ && needVel_){
250			info_->getComAll(com, comvel);
251			comw = info_->getAngularMomentum();
252			}else{
253			com = info_->getCom();
254			comvel = 0.0;
255			comw = 0.0;
256			}
257	gezelter	1104	s->setCOMprops(com, comvel, comw);
258			}
259
260	chrisfen	721	}
261
262	tim	1024	void DumpReader::readSet(int whichFrame) {
263			std::string line;
264
265	chrisfen	721	#ifndef IS_MPI
266	tim	1024
267	chrisfen	721	inFile_->clear();
268			inFile_->seekg(framePos_[whichFrame]);
269	tim	1024
270			std::istream& inputStream = *inFile_;
271
272			#else
273			int masterNode = 0;
274			std::stringstream sstream;
275			if (worldRank == masterNode) {
276			std::string sendBuffer;
277
278			inFile_->clear();
279			inFile_->seekg(framePos_[whichFrame]);
280
281			while (inFile_->getline(buffer, bufferSize)) {
282
283			line = buffer;
284			sendBuffer += line;
285			sendBuffer += '\n';
286			if (line.find("</Snapshot>") != std::string::npos) {
287			break;
288			}
289			}
290
291			int sendBufferSize = sendBuffer.size();
292			MPI_Bcast(&sendBufferSize, 1, MPI_INT, masterNode, MPI_COMM_WORLD);
293			MPI_Bcast((void *)sendBuffer.c_str(), sendBufferSize, MPI_CHAR, masterNode, MPI_COMM_WORLD);
294
295			sstream.str(sendBuffer);
296			} else {
297			int sendBufferSize;
298			MPI_Bcast(&sendBufferSize, 1, MPI_INT, masterNode, MPI_COMM_WORLD);
299			char * recvBuffer = new char[sendBufferSize+1];
300			MPI_Bcast(recvBuffer, sendBufferSize, MPI_CHAR, masterNode, MPI_COMM_WORLD);
301			sstream.str(recvBuffer);
302			}
303
304			std::istream& inputStream = sstream;
305			#endif
306
307			inputStream.getline(buffer, bufferSize);
308
309			line = buffer;
310			if (line.find("<Snapshot>") == std::string::npos) {
311	chrisfen	721	sprintf(painCave.errMsg,
312	tim	1024	"DumpReader Error: can not find <Snapshot>\n");
313	chrisfen	721	painCave.isFatal = 1;
314			simError();
315			}
316	tim	1024
317			//read frameData
318			readFrameProperties(inputStream);
319
320			//read StuntDoubles
321			readStuntDoubles(inputStream);
322
323			inputStream.getline(buffer, bufferSize);
324			line = buffer;
325			if (line.find("</Snapshot>") == std::string::npos) {
326	chrisfen	721	sprintf(painCave.errMsg,
327	tim	1024	"DumpReader Error: can not find </Snapshot>\n");
328	chrisfen	721	painCave.isFatal = 1;
329			simError();
330	tim	1024	}
331
332	chrisfen	721	}
333
334	tim	1024	void DumpReader::parseDumpLine(const std::string& line) {
335
336
337	chrisfen	721	StringTokenizer tokenizer(line);
338			int nTokens;
339
340			nTokens = tokenizer.countTokens();
341
342	tim	1024	if (nTokens < 2) {
343	chrisfen	721	sprintf(painCave.errMsg,
344	tim	1024	"DumpReader Error: Not enough Tokens.\n%s\n", line.c_str());
345	chrisfen	721	painCave.isFatal = 1;
346			simError();
347			}
348	tim	1024
349			int index = tokenizer.nextTokenAsInt();
350
351			StuntDouble* integrableObject = info_->getIOIndexToIntegrableObject(index);
352
353			if (integrableObject == NULL) {
354			return;
355			}
356			std::string type = tokenizer.nextToken();
357			int size = type.size();
358	gezelter	1037
359	tim	1024	for(int i = 0; i < size; ++i) {
360			switch(type[i]) {
361
362			case 'p': {
363			Vector3d pos;
364			pos[0] = tokenizer.nextTokenAsDouble();
365			pos[1] = tokenizer.nextTokenAsDouble();
366			pos[2] = tokenizer.nextTokenAsDouble();
367			if (needPos_) {
368			integrableObject->setPos(pos);
369			}
370			break;
371			}
372			case 'v' : {
373			Vector3d vel;
374			vel[0] = tokenizer.nextTokenAsDouble();
375			vel[1] = tokenizer.nextTokenAsDouble();
376			vel[2] = tokenizer.nextTokenAsDouble();
377			if (needVel_) {
378			integrableObject->setVel(vel);
379			}
380			break;
381			}
382
383			case 'q' : {
384			Quat4d q;
385			if (integrableObject->isDirectional()) {
386
387			q[0] = tokenizer.nextTokenAsDouble();
388			q[1] = tokenizer.nextTokenAsDouble();
389			q[2] = tokenizer.nextTokenAsDouble();
390			q[3] = tokenizer.nextTokenAsDouble();
391
392			RealType qlen = q.length();
393			if (qlen < oopse::epsilon) { //check quaternion is not equal to 0
394
395			sprintf(painCave.errMsg,
396			"DumpReader Error: initial quaternion error (q0^2 + q1^2 + q2^2 + q3^2) ~ 0\n");
397			painCave.isFatal = 1;
398			simError();
399
400			}
401
402			q.normalize();
403			if (needQuaternion_) {
404			integrableObject->setQ(q);
405			}
406			}
407			break;
408			}
409			case 'j' : {
410			Vector3d ji;
411			if (integrableObject->isDirectional()) {
412			ji[0] = tokenizer.nextTokenAsDouble();
413			ji[1] = tokenizer.nextTokenAsDouble();
414			ji[2] = tokenizer.nextTokenAsDouble();
415			if (needAngMom_) {
416			integrableObject->setJ(ji);
417			}
418			}
419	gezelter	1037	break;
420	tim	1024	}
421			case 'f': {
422
423			Vector3d force;
424			force[0] = tokenizer.nextTokenAsDouble();
425			force[1] = tokenizer.nextTokenAsDouble();
426			force[2] = tokenizer.nextTokenAsDouble();
427			integrableObject->setFrc(force);
428			break;
429			}
430			case 't' : {
431
432			Vector3d torque;
433			torque[0] = tokenizer.nextTokenAsDouble();
434			torque[1] = tokenizer.nextTokenAsDouble();
435			torque[2] = tokenizer.nextTokenAsDouble();
436			integrableObject->setTrq(torque);
437			break;
438			}
439			default: {
440			sprintf(painCave.errMsg,
441			"DumpReader Error: %s is an unrecognized type\n", type.c_str());
442			painCave.isFatal = 1;
443			simError();
444			break;
445			}
446
447			}
448			}
449	chrisfen	721
450	tim	1024	}
451
452
453			void DumpReader::readStuntDoubles(std::istream& inputStream) {
454
455			inputStream.getline(buffer, bufferSize);
456			std::string line(buffer);
457
458			if (line.find("<StuntDoubles>") == std::string::npos) {
459	chrisfen	721	sprintf(painCave.errMsg,
460	tim	1024	"DumpReader Error: Missing <StuntDoubles>\n");
461	chrisfen	721	painCave.isFatal = 1;
462	tim	1024	simError();
463			}
464
465			while(inputStream.getline(buffer, bufferSize)) {
466			line = buffer;
467
468			if(line.find("</StuntDoubles>") != std::string::npos) {
469			break;
470			}
471
472			parseDumpLine(line);
473			}
474
475			}
476
477			void DumpReader::readFrameProperties(std::istream& inputStream) {
478
479			Snapshot* s = info_->getSnapshotManager()->getCurrentSnapshot();
480			inputStream.getline(buffer, bufferSize);
481			std::string line(buffer);
482
483			if (line.find("<FrameData>") == std::string::npos) {
484			sprintf(painCave.errMsg,
485			"DumpReader Error: Missing <FrameData>\n");
486			painCave.isFatal = 1;
487			simError();
488			}
489
490			while(inputStream.getline(buffer, bufferSize)) {
491			line = buffer;
492
493			if(line.find("</FrameData>") != std::string::npos) {
494			break;
495			}
496
497			StringTokenizer tokenizer(line, " ;\t\n\r{}:,");
498			if (!tokenizer.hasMoreTokens()) {
499	chrisfen	721	sprintf(painCave.errMsg,
500	tim	1024	"DumpReader Error: Not enough Tokens.\n%s\n", line.c_str());
501	chrisfen	721	painCave.isFatal = 1;
502	tim	1024	simError();
503			}
504
505			std::string propertyName = tokenizer.nextToken();
506			if (propertyName == "Time") {
507			RealType currTime = tokenizer.nextTokenAsDouble();
508			s->setTime(currTime);
509			} else if (propertyName == "Hmat"){
510			Mat3x3d hmat;
511			hmat(0, 0) = tokenizer.nextTokenAsDouble();
512			hmat(0, 1) = tokenizer.nextTokenAsDouble();
513			hmat(0, 2) = tokenizer.nextTokenAsDouble();
514			hmat(1, 0) = tokenizer.nextTokenAsDouble();
515			hmat(1, 1) = tokenizer.nextTokenAsDouble();
516			hmat(1, 2) = tokenizer.nextTokenAsDouble();
517			hmat(2, 0) = tokenizer.nextTokenAsDouble();
518			hmat(2, 1) = tokenizer.nextTokenAsDouble();
519			hmat(2, 2) = tokenizer.nextTokenAsDouble();
520			s->setHmat(hmat);
521			} else if (propertyName == "Thermostat") {
522			RealType chi = tokenizer.nextTokenAsDouble();
523			RealType integralOfChiDt = tokenizer.nextTokenAsDouble();
524			s->setChi(chi);
525			s->setIntegralOfChiDt(integralOfChiDt);
526			} else if (propertyName == "Barostat") {
527			Mat3x3d eta;
528			eta(0, 0) = tokenizer.nextTokenAsDouble();
529			eta(0, 1) = tokenizer.nextTokenAsDouble();
530			eta(0, 2) = tokenizer.nextTokenAsDouble();
531			eta(1, 0) = tokenizer.nextTokenAsDouble();
532			eta(1, 1) = tokenizer.nextTokenAsDouble();
533			eta(1, 2) = tokenizer.nextTokenAsDouble();
534			eta(2, 0) = tokenizer.nextTokenAsDouble();
535			eta(2, 1) = tokenizer.nextTokenAsDouble();
536			eta(2, 2) = tokenizer.nextTokenAsDouble();
537			s->setEta(eta);
538			} else {
539			sprintf(painCave.errMsg,
540			"DumpReader Error: %s is an invalid property in <FrameData>\n", propertyName.c_str());
541			painCave.isFatal = 0;
542			simError();
543			}
544
545			}
546
547			}
548
549	chrisfen	721
550			}//end namespace oopse