ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/branches/development/src/io/RestWriter.cpp
(Generate patch)

Comparing trunk/src/io/RestWriter.cpp (file contents):
Revision 990 by chrisfen, Mon Jun 19 01:36:06 2006 UTC vs.
Revision 1364 by cli2, Wed Oct 7 20:49:50 2009 UTC

# Line 1 | Line 1
1   /*
2 < * Copyright (c) 2005 The University of Notre Dame. All Rights Reserved.
2 > * Copyright (c) 2009 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
# Line 38 | Line 38
38   * University of Notre Dame has been advised of the possibility of
39   * such damages.
40   */
41 <
42 < #include <algorithm>
41 >
42 >
43   #include <iostream>
44 – #include <map>
44  
46 – #include "primitives/Molecule.hpp"
45   #include "io/RestWriter.hpp"
46   #include "utils/simError.h"
47 + #include "brains/SnapshotManager.hpp"
48 + #ifdef IS_MPI
49 + #include <mpi.h>
50 + #endif
51  
50 –
52   namespace oopse {
53 <  RestWriter::RestWriter(SimInfo* info) :
54 <    info_(info) {
53 >  RestWriter::RestWriter(SimInfo* info, const std::string& filename,
54 >                         std::vector<Restraint*> restraints ) :
55 >    info_(info){
56  
57 <    // only the master node writes to the disk
57 >    //use master - slave mode, only master node writes to disk
58   #ifdef IS_MPI
59 <    if (worldRank == 0) {
60 < #endif // is_mpi
59 >    if(worldRank == 0){
60 > #endif
61 >      
62 >      output_.open(filename.c_str());
63 >      
64 >      if(!output_){
65 >        sprintf( painCave.errMsg,
66 >                 "Could not open %s for restraint output.\n",
67 >                 filename.c_str());
68 >        painCave.isFatal = 1;
69 >        simError();
70 >      }
71  
72 <      outName = info_->getRestFileName();
72 >      output_ << "#time\t";
73 >      
74 >      // TODO:  get Restraint info from slave nodes:
75 >      std::vector<Restraint*>::const_iterator resti;
76 >      for(resti=restraints.begin(); resti != restraints.end(); ++resti){
77  
78 +        if ((*resti)->getPrintRestraint()) {
79 +          std::string myName = (*resti)->getRestraintName();
80 +          int myType = (*resti)->getRestraintType();
81 +          
82 +          output_ << myName << ":";
83 +          
84 +          if (myType & Restraint::rtDisplacement)
85 +            output_ << "\tPosition(angstroms)\tEnergy(kcal/mol)";
86 +          
87 +          if (myType & Restraint::rtTwist)
88 +            output_ << "\tTwistAngle(radians)\tEnergy(kcal/mol)";
89 +          
90 +          if (myType & Restraint::rtSwingX)
91 +            output_ << "\tSwingXAngle(radians)\tEnergy(kcal/mol)";
92 +          
93 +          if (myType & Restraint::rtSwingY)
94 +            output_ << "\tSwingYAngle(radians)\tEnergy(kcal/mol)";
95 +          
96 +        }
97 +      }
98 +      output_ << "\n";
99   #ifdef IS_MPI
100      }
101 < #endif // is_mpi
101 > #endif      
102    }
103    
104 <  RestWriter::~RestWriter() {}
68 <  
69 <  void RestWriter::writeZangle(){
70 <    const int BUFFERSIZE = 2000;
71 <    char tempBuffer[BUFFERSIZE];
72 <    char writeLine[BUFFERSIZE];
73 <    
74 <    std::ofstream finalOut;
75 <    
76 <    Molecule* mol;
77 <    StuntDouble* integrableObject;
78 <    SimInfo::MoleculeIterator mi;
79 <    Molecule::IntegrableObjectIterator ii;
80 <
104 >  RestWriter::~RestWriter() {
105   #ifdef IS_MPI
106      if(worldRank == 0 ){
107 < #endif    
108 <      finalOut.open( outName.c_str(), std::ios::out | std::ios::trunc );
109 <      if( !finalOut ){
86 <        sprintf( painCave.errMsg,
87 <                 "Could not open \"%s\" for zAngle output.\n",
88 <                 outName.c_str() );
89 <        painCave.isFatal = 1;
90 <        simError();
91 <      }
92 < #ifdef IS_MPI
107 > #endif  
108 >      output_.close();  
109 > #ifdef IS_MPI  
110      }
111 < #endif // is_mpi
112 <    
113 < #ifndef IS_MPI
114 <    // first we do output for the single processor version
98 <    finalOut
99 <      << info_->getSnapshotManager()->getCurrentSnapshot()->getTime()
100 <      << " : omega values at this time\n";
101 <    
102 <    for (mol = info_->beginMolecule(mi); mol != NULL;
103 <         mol = info_->nextMolecule(mi)) {
111 > #endif
112 >  }
113 >  
114 >  void RestWriter::writeRest(std::vector<std::map<int, Restraint::RealPair> > restInfo){
115        
105 –      for (integrableObject = mol->beginIntegrableObject(ii);
106 –           integrableObject != NULL;
107 –           integrableObject = mol->nextIntegrableObject(ii)) {    
108 –        
109 –        sprintf( tempBuffer,
110 –                 "%14.10lf\n",
111 –                 integrableObject->getZangle());
112 –        strcpy( writeLine, tempBuffer );
113 –        
114 –        finalOut << writeLine;      
115 –      }
116 –      
117 –    }
116      
117 < #else
120 <    int nproc;
121 <    MPI_Comm_size(MPI_COMM_WORLD, &nproc);
122 <    const int masterNode = 0;
123 <    int myNode = worldRank;
124 <    std::vector<int> tmpNIntObjects(nproc, 0);
125 <    std::vector<int> nIntObjectsInProc(nproc, 0);
126 <    tmpNIntObjects[myNode] = info_->getNGlobalIntegrableObjects();
127 <    
128 <    //do MPI_ALLREDUCE to exchange the total number of atoms, rigidbodies and cutoff groups
129 <    MPI_Allreduce(&tmpNIntObjects[0], &nIntObjectsInProc[0], nproc, MPI_INT,
130 <                  MPI_SUM, MPI_COMM_WORLD);
131 <    
132 <    MPI_Status ierr;
133 <    int intObIndex;
134 <    RealType zAngle;
135 <  
136 <    if (worldRank == masterNode) {
137 <      std::map<int, RealType> zAngData;
138 <      for(int i = 0 ; i < nproc; ++i) {
139 <        if (i == masterNode) {
140 <          for (mol = info_->beginMolecule(mi); mol != NULL;
141 <               mol = info_->nextMolecule(mi)) {
142 <            
143 <            for (integrableObject = mol->beginIntegrableObject(ii);
144 <                 integrableObject != NULL;
145 <                 integrableObject = mol->nextIntegrableObject(ii)) {
146 <              
147 <              intObIndex = integrableObject->getGlobalIndex();
148 <
149 <              zAngle = integrableObject->getZangle();
150 <              zAngData.insert(std::pair<int, RealType>(intObIndex, zAngle));
151 <            }      
152 <          }
153 <        } else {
154 <          for (mol = info_->beginMolecule(mi); mol != NULL;
155 <               mol = info_->nextMolecule(mi)) {
156 <            
157 <            for (integrableObject = mol->beginIntegrableObject(ii);
158 <                 integrableObject != NULL;
159 <                 integrableObject = mol->nextIntegrableObject(ii)) {
160 <              
161 <              MPI_Recv(&intObIndex, 1, MPI_INT, i, 0, MPI_COMM_WORLD,&ierr);
162 <              MPI_Recv(&zAngle, 1, MPI_REALTYPE, i, 0, MPI_COMM_WORLD,&ierr);
163 <              zAngData.insert(std::pair<int, RealType>(intObIndex, zAngle));
164 <            }
165 <          }
166 <        }
167 <      }
168 <      
169 <      finalOut
170 <        << info_->getSnapshotManager()->getCurrentSnapshot()->getTime()
171 <        << " : omega values at this time\n";
172 <      
173 <      std::map<int, RealType>::iterator l;
174 <      for (l = zAngData.begin(); l != zAngData.end(); ++l) {
175 <        sprintf( tempBuffer,
176 <                 "%14.10lf\n",
177 <                 l->second);
178 <        strcpy( writeLine, tempBuffer );
117 >    output_ << info_->getSnapshotManager()->getCurrentSnapshot()->getTime();
118          
119 <        finalOut << writeLine;      
119 >    // output some information about the molecules
120 >    std::vector<std::map<int, Restraint::RealPair> >::const_iterator i;
121 >    std::map<int, Restraint::RealPair>::const_iterator j;
122 >    for( i = restInfo.begin(); i != restInfo.end(); ++i){
123 >      for(j = (*i).begin(); j != (*i).end(); ++j){                
124 >        output_ << "\t" << (j->second).first << "\t" << (j->second).second;
125        }
126 <
183 <      finalOut.close();
184 <      
185 <    } else {
186 <      for(int j = 1; j < nproc; ++j) {
187 <        for (mol = info_->beginMolecule(mi); mol != NULL;
188 <             mol = info_->nextMolecule(mi)) {
189 <          
190 <          for (integrableObject = mol->beginIntegrableObject(ii);
191 <               integrableObject != NULL;
192 <               integrableObject = mol->nextIntegrableObject(ii)) {
193 <            intObIndex = integrableObject->getGlobalIndex();            
194 <            zAngle = integrableObject->getZangle();
195 <
196 <            MPI_Send(&intObIndex, 1, MPI_INT, masterNode, 0, MPI_COMM_WORLD);
197 <            MPI_Send(&zAngle, 1, MPI_REALTYPE, masterNode, 0, MPI_COMM_WORLD);
198 <          }
199 <        }
200 <      }
126 >      output_ << std::endl;
127      }
202 –
203 – #endif
128    }
129    
130 < }
130 > }// end oopse
131 >  

Diff Legend

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