ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/mdtools/interface_implementation/SimSetup.cpp
(Generate patch)

Comparing trunk/mdtools/interface_implementation/SimSetup.cpp (file contents):
Revision 139 by chuckv, Wed Oct 16 21:07:16 2002 UTC vs.
Revision 198 by mmeineke, Fri Dec 6 21:20:42 2002 UTC

# Line 6 | Line 6
6   #include "parse_me.h"
7   #include "LRI.hpp"
8   #include "Integrator.hpp"
9 + #include "simError.h"
10 +
11   #ifdef IS_MPI
12   #include "mpiBASS.h"
13   #include "bassDiag.hpp"
# Line 14 | Line 16 | SimSetup::SimSetup(){
16   SimSetup::SimSetup(){
17    stamps = new MakeStamps();
18    globals = new Globals();
19 +  
20 + #ifdef IS_MPI
21 +  strcpy( checkPointMsg, "SimSetup creation successful" );
22 +  MPIcheckPoint();
23 + #endif // IS_MPI
24   }
25  
26   SimSetup::~SimSetup(){
# Line 23 | Line 30 | void SimSetup::parseFile( char* fileName ){
30  
31   void SimSetup::parseFile( char* fileName ){
32  
26 –  inFileName = fileName;
27 –  set_interface_stamps( stamps, globals );
33   #ifdef IS_MPI
34 <  mpiEventInit();
34 >  if( worldRank == 0 ){
35 > #endif // is_mpi
36 >    
37 >    inFileName = fileName;
38 >    set_interface_stamps( stamps, globals );
39 >    
40 > #ifdef IS_MPI
41 >    mpiEventInit();
42   #endif
43 <  yacc_BASS( fileName );
43 >
44 >    yacc_BASS( fileName );
45 >
46   #ifdef IS_MPI
47 <  throwMPIEvent(NULL);
47 >    throwMPIEvent(NULL);
48 >  }
49 >  else receiveParse();
50   #endif
51  
52   }
# Line 40 | Line 56 | void SimSetup::receiveParse(void){
56  
57      set_interface_stamps( stamps, globals );
58      mpiEventInit();
59 +    MPIcheckPoint();
60      mpiEventLoop();
61  
62   }
# Line 51 | Line 68 | void SimSetup::testMe(void){
68    delete dumpMe;
69   }
70   #endif
71 +
72   void SimSetup::createSim( void ){
73  
74    MakeStamps *the_stamps;
# Line 75 | Line 93 | void SimSetup::createSim( void ){
93    else if( !strcmp( force_field, "DipoleTest" ) ) the_ff = new DipoleTestFF();
94    else if( !strcmp( force_field, "TraPPE_Ex" ) ) the_ff = new TraPPE_ExFF();
95    else{
96 <    std::cerr<< "SimSetup Error. Unrecognized force field -> "
97 <             << force_field << "\n";
98 <    exit(8);
96 >    sprintf( painCave.errMsg,
97 >             "SimSetup Error. Unrecognized force field -> %s\n",
98 >             force_field );
99 >    painCave.isFatal = 1;
100 >    simError();
101    }
102  
103 + #ifdef IS_MPI
104 +  strcpy( checkPointMsg, "ForceField creation successful" );
105 +  MPIcheckPoint();
106 + #endif // is_mpi
107 +
108    // get the components and calculate the tot_nMol and indvidual n_mol
109    the_components = the_globals->getComponents();
110    components_nmol = new int[n_components];
# Line 94 | Line 119 | void SimSetup::createSim( void ){
119  
120        if( !the_components[i]->haveNMol() ){
121          // we have a problem
122 <        std::cerr << "SimSetup Error. No global NMol or component NMol"
123 <                  << " given. Cannot calculate the number of atoms.\n";
124 <        exit( 8 );
122 >        sprintf( painCave.errMsg,
123 >                 "SimSetup Error. No global NMol or component NMol"
124 >                 " given. Cannot calculate the number of atoms.\n" );
125 >        painCave.isFatal = 1;
126 >        simError();
127        }
128  
129        tot_nmol += the_components[i]->getNMol();
# Line 104 | Line 131 | void SimSetup::createSim( void ){
131      }
132    }
133    else{
134 <    std::cerr << "NOT A SUPPORTED FEATURE\n";
135 <
136 < //     tot_nmol = the_globals->getNMol();
137 <
138 < //     //we have the total number of molecules, now we check for molfractions
139 < //     for( i=0; i<n_components; i++ ){
140 <
141 < //       if( !the_components[i]->haveMolFraction() ){
142 <
143 < //      if( !the_components[i]->haveNMol() ){
144 < //        //we have a problem
145 < //        std::cerr << "SimSetup error. Neither molFraction nor "
146 < //                  << " nMol was given in component
147 <
134 >    sprintf( painCave.errMsg,
135 >             "SimSetup error.\n"
136 >             "\tSorry, the ability to specify total"
137 >             " nMols and then give molfractions in the components\n"
138 >             "\tis not currently supported."
139 >             " Please give nMol in the components.\n" );
140 >    painCave.isFatal = 1;
141 >    simError();
142 >    
143 >    
144 >    //     tot_nmol = the_globals->getNMol();
145 >    
146 >    //   //we have the total number of molecules, now we check for molfractions
147 >    //     for( i=0; i<n_components; i++ ){
148 >    
149 >    //       if( !the_components[i]->haveMolFraction() ){
150 >    
151 >    //  if( !the_components[i]->haveNMol() ){
152 >    //    //we have a problem
153 >    //    std::cerr << "SimSetup error. Neither molFraction nor "
154 >    //              << " nMol was given in component
155 >    
156    }
157  
158 + #ifdef IS_MPI
159 +  strcpy( checkPointMsg, "Have the number of components" );
160 +  MPIcheckPoint();
161 + #endif // is_mpi
162 +
163    // make an array of molecule stamps that match the components used.
164 +  // also extract the used stamps out into a separate linked list
165  
166 +  simnfo->nComponents = n_components;
167 +  simnfo->componentsNmol = components_nmol;
168 +  simnfo->compStamps = comp_stamps;
169 +  simnfo->headStamp = new LinkedMolStamp();
170 +  
171 +  char* id;
172 +  LinkedMolStamp* headStamp = simnfo->headStamp;
173 +  LinkedMolStamp* currentStamp = NULL;
174    for( i=0; i<n_components; i++ ){
175  
176 <    comp_stamps[i] =
177 <      the_stamps->getMolecule( the_components[i]->getType() );
176 >    id = the_components[i]->getType();
177 >    comp_stamps[i] = NULL;
178 >    
179 >    // check to make sure the component isn't already in the list
180 >
181 >    comp_stamps[i] = headStamp->match( id );
182 >    if( comp_stamps[i] == NULL ){
183 >      
184 >      // extract the component from the list;
185 >      
186 >      currentStamp = the_stamps->extractMolStamp( id );
187 >      if( currentStamp == NULL ){
188 >        sprintf( painCave.errMsg,
189 >                 "SimSetup error: Component \"%s\" was not found in the "
190 >                 "list of declared molecules\n"
191 >                 id );
192 >        painCave.isFatal = 1;
193 >        simError();
194 >      }
195 >      
196 >      headStamp->add( currentStamp );
197 >      comp_stamps[i] = headStamp->match( id );
198 >    }
199    }
200  
201 + #ifdef IS_MPI
202 +  strcpy( checkPointMsg, "Component stamps successfully extracted\n" );
203 +  MPIcheckPoint();
204 + #endif // is_mpi
205 +  
206  
207  
208 +
209    // caclulate the number of atoms, bonds, bends and torsions
210  
211    tot_atoms = 0;
# Line 151 | Line 227 | void SimSetup::createSim( void ){
227    simnfo->n_bends = tot_bends;
228    simnfo->n_torsions = tot_torsions;
229    simnfo->n_SRI = tot_SRI;
230 +  simnfo->n_mol = tot_nmol;
231  
232 +  
233 + #ifdef IS_MPI
234 +
235 +  // divide the molecules among processors here.
236 +  
237 +  new mpiSimulation( simnfo );
238 +  
239 +  simnfo->mpiSim->divideLabor( n_components, comp_stamps, components_nmol );
240 +
241 + #endif // is_mpi
242 +  
243 +
244    // create the atom and short range interaction arrays
245  
246 <  the_atoms = new Atom*[tot_atoms];
247 <  the_molecules = new Molecule[tot_nmol];
246 >  Atom::createArrays(simnfo->n_atoms);
247 >  the_atoms = new Atom*[simnfo->n_atoms];
248 >  the_molecules = new Molecule[simnfo->n_mol];
249  
250  
251 <  if( tot_SRI ){
252 <    the_sris = new SRI*[tot_SRI];
253 <    the_excludes = new ex_pair[tot_SRI];
251 >  if( simnfo->n_SRI ){
252 >    the_sris = new SRI*[simnfo->n_SRI];
253 >    the_excludes = new ex_pair[simnfo->n_SRI];
254    }
255  
256    // set the arrays into the SimInfo object
# Line 189 | Line 279 | void SimSetup::createSim( void ){
279      makeTorsions();
280    }
281  
192 –  //  makeMolecules();
282  
283    // get some of the tricky things that may still be in the globals
284  
285    if( simnfo->n_dipoles ){
286  
287      if( !the_globals->haveRRF() ){
288 <      std::cerr << "SimSetup Error, system has dipoles, but no rRF was set.\n";
289 <      exit(8);
288 >      sprintf( painCave.errMsg,
289 >               "SimSetup Error, system has dipoles, but no rRF was set.\n");
290 >      painCave.isFatal = 1;
291 >      simError();
292      }
293      if( !the_globals->haveDielectric() ){
294 <      std::cerr << "SimSetup Error, system has dipoles, but no"
295 <                << " dielectric was set.\n";
296 <      exit(8);
294 >      sprintf( painCave.errMsg,
295 >               "SimSetup Error, system has dipoles, but no"
296 >               " dielectric was set.\n" );
297 >      painCave.isFatal = 1;
298 >      simError();
299      }
300  
301      simnfo->rRF        = the_globals->getRRF();
302      simnfo->dielectric = the_globals->getDielectric();
303    }
304  
305 + #ifdef IS_MPI
306 +  strcpy( checkPointMsg, "rRf and dielectric check out" );
307 +  MPIcheckPoint();
308 + #endif // is_mpi
309 +  
310    if( the_globals->haveBox() ){
311      simnfo->box_x = the_globals->getBox();
312      simnfo->box_y = the_globals->getBox();
# Line 224 | Line 322 | void SimSetup::createSim( void ){
322    }
323    else{
324      if( !the_globals->haveBoxX() ){
325 <      std::cerr << "SimSetup error, no periodic BoxX size given.\n";
326 <      exit(8);
325 >      sprintf( painCave.errMsg,
326 >               "SimSetup error, no periodic BoxX size given.\n" );
327 >      painCave.isFatal = 1;
328 >      simError();
329      }
330      simnfo->box_x = the_globals->getBoxX();
331  
332      if( !the_globals->haveBoxY() ){
333 <      std::cerr << "SimSetup error, no periodic BoxY size given.\n";
334 <      exit(8);
333 >      sprintf( painCave.errMsg,
334 >               "SimSetup error, no periodic BoxY size given.\n" );
335 >      painCave.isFatal = 1;
336 >      simError();
337      }
338      simnfo->box_y = the_globals->getBoxY();
339  
340      if( !the_globals->haveBoxZ() ){
341 <      std::cerr << "SimSetup error, no periodic BoxZ size given.\n";
342 <      exit(8);
341 >      sprintf( painCave.errMsg,
342 >               "SimSetup error, no periodic BoxZ size given.\n" );
343 >      painCave.isFatal = 1;
344 >      simError();
345      }
346      simnfo->box_z = the_globals->getBoxZ();
347    }
348  
349 + #ifdef IS_MPI
350 +  strcpy( checkPointMsg, "Box size set up" );
351 +  MPIcheckPoint();
352 + #endif // is_mpi
353  
354 +
355 +
356   //   if( the_globals->haveInitialConfig() ){
357   //        InitializeFromFile* fileInit;
358   //     fileInit = new InitializeFromFile( the_globals->getInitialConfig() );
# Line 253 | Line 363 | void SimSetup::createSim( void ){
363   //   }
364   //   else{
365  
366 <    initFromBass();
366 > #ifdef IS_MPI
367  
368 <
369 < //   }
368 >  // no init from bass
369 >  
370 >  sprintf( painCave.errMsg,
371 >           "Cannot intialize a parallel simulation without an initial configuration file.\n" );
372 >  painCave.isFatal;
373 >  simError();
374 >  
375 > #else
376  
377 < //   if( the_globals->haveFinalConfig() ){
262 < //     strcpy( simnfo->finalName, the_globals->getFinalConfig() );
263 < //   }
264 < //   else{
265 < //     strcpy( simnfo->finalName, inFileName );
266 < //     char* endTest;
267 < //     int nameLength = strlen( simnfo->finalName );
268 < //     endTest = &(simnfo->finalName[nameLength - 5]);
269 < //     if( !strcmp( endTest, ".bass" ) ){
270 < //       strcpy( endTest, ".eor" );
271 < //     }
272 < //     else if( !strcmp( endTest, ".BASS" ) ){
273 < //       strcpy( endTest, ".eor" );
274 < //     }
275 < //     else{
276 < //       endTest = &(simnfo->finalName[nameLength - 4]);
277 < //       if( !strcmp( endTest, ".bss" ) ){
278 < //      strcpy( endTest, ".eor" );
279 < //       }
280 < //       else if( !strcmp( endTest, ".mdl" ) ){
281 < //      strcpy( endTest, ".eor" );
282 < //       }
283 < //       else{
284 < //      strcat( simnfo->finalName, ".eor" );
285 < //       }
286 < //     }
287 < //   }
377 >  initFromBass();
378  
379 < //   // make the sample and status out names
379 > #endif // is_mpi
380 >
381 > #ifdef IS_MPI
382 >  strcpy( checkPointMsg, "Successfully read in the initial configuration" );
383 >  MPIcheckPoint();
384 > #endif // is_mpi
385  
291 – //   strcpy( simnfo->sampleName, inFileName );
292 – //   char* endTest;
293 – //   int nameLength = strlen( simnfo->sampleName );
294 – //   endTest = &(simnfo->sampleName[nameLength - 5]);
295 – //   if( !strcmp( endTest, ".bass" ) ){
296 – //     strcpy( endTest, ".dump" );
297 – //   }
298 – //   else if( !strcmp( endTest, ".BASS" ) ){
299 – //     strcpy( endTest, ".dump" );
300 – //   }
301 – //   else{
302 – //     endTest = &(simnfo->sampleName[nameLength - 4]);
303 – //     if( !strcmp( endTest, ".bss" ) ){
304 – //       strcpy( endTest, ".dump" );
305 – //     }
306 – //     else if( !strcmp( endTest, ".mdl" ) ){
307 – //       strcpy( endTest, ".dump" );
308 – //     }
309 – //     else{
310 – //       strcat( simnfo->sampleName, ".dump" );
311 – //     }
312 – //   }
386  
387 < //   strcpy( simnfo->statusName, inFileName );
388 < //   nameLength = strlen( simnfo->statusName );
389 < //   endTest = &(simnfo->statusName[nameLength - 5]);
390 < //   if( !strcmp( endTest, ".bass" ) ){
391 < //     strcpy( endTest, ".stat" );
392 < //   }
393 < //   else if( !strcmp( endTest, ".BASS" ) ){
394 < //     strcpy( endTest, ".stat" );
395 < //   }
396 < //   else{
397 < //     endTest = &(simnfo->statusName[nameLength - 4]);
398 < //     if( !strcmp( endTest, ".bss" ) ){
399 < //       strcpy( endTest, ".stat" );
400 < //     }
401 < //     else if( !strcmp( endTest, ".mdl" ) ){
402 < //       strcpy( endTest, ".stat" );
403 < //     }
404 < //     else{
405 < //       strcat( simnfo->statusName, ".stat" );
406 < //     }
407 < //   }
408 <
409 <
387 >  
388 >
389 >  
390 >  //   }
391 >  
392 > #ifdef IS_MPI
393 >  if( worldRank == 0 ){
394 > #endif // is_mpi
395 >    
396 >    if( the_globals->haveFinalConfig() ){
397 >      strcpy( simnfo->finalName, the_globals->getFinalConfig() );
398 >    }
399 >    else{
400 >      strcpy( simnfo->finalName, inFileName );
401 >      char* endTest;
402 >      int nameLength = strlen( simnfo->finalName );
403 >      endTest = &(simnfo->finalName[nameLength - 5]);
404 >      if( !strcmp( endTest, ".bass" ) ){
405 >        strcpy( endTest, ".eor" );
406 >      }
407 >      else if( !strcmp( endTest, ".BASS" ) ){
408 >        strcpy( endTest, ".eor" );
409 >      }
410 >      else{
411 >        endTest = &(simnfo->finalName[nameLength - 4]);
412 >        if( !strcmp( endTest, ".bss" ) ){
413 >          strcpy( endTest, ".eor" );
414 >        }
415 >        else if( !strcmp( endTest, ".mdl" ) ){
416 >          strcpy( endTest, ".eor" );
417 >        }
418 >        else{
419 >          strcat( simnfo->finalName, ".eor" );
420 >        }
421 >      }
422 >    }
423 >    
424 >    // make the sample and status out names
425 >    
426 >    strcpy( simnfo->sampleName, inFileName );
427 >    char* endTest;
428 >    int nameLength = strlen( simnfo->sampleName );
429 >    endTest = &(simnfo->sampleName[nameLength - 5]);
430 >    if( !strcmp( endTest, ".bass" ) ){
431 >      strcpy( endTest, ".dump" );
432 >    }
433 >    else if( !strcmp( endTest, ".BASS" ) ){
434 >      strcpy( endTest, ".dump" );
435 >    }
436 >    else{
437 >      endTest = &(simnfo->sampleName[nameLength - 4]);
438 >      if( !strcmp( endTest, ".bss" ) ){
439 >        strcpy( endTest, ".dump" );
440 >      }
441 >      else if( !strcmp( endTest, ".mdl" ) ){
442 >        strcpy( endTest, ".dump" );
443 >      }
444 >      else{
445 >        strcat( simnfo->sampleName, ".dump" );
446 >      }
447 >    }
448 >    
449 >    strcpy( simnfo->statusName, inFileName );
450 >    nameLength = strlen( simnfo->statusName );
451 >    endTest = &(simnfo->statusName[nameLength - 5]);
452 >    if( !strcmp( endTest, ".bass" ) ){
453 >      strcpy( endTest, ".stat" );
454 >    }
455 >    else if( !strcmp( endTest, ".BASS" ) ){
456 >      strcpy( endTest, ".stat" );
457 >    }
458 >    else{
459 >      endTest = &(simnfo->statusName[nameLength - 4]);
460 >      if( !strcmp( endTest, ".bss" ) ){
461 >        strcpy( endTest, ".stat" );
462 >      }
463 >      else if( !strcmp( endTest, ".mdl" ) ){
464 >        strcpy( endTest, ".stat" );
465 >      }
466 >      else{
467 >        strcat( simnfo->statusName, ".stat" );
468 >      }
469 >    }
470 >    
471 > #ifdef IS_MPI
472 >  }
473 > #endif // is_mpi
474 >  
475    // set the status, sample, and themal kick times
476 <
476 >  
477    if( the_globals->haveSampleTime() ){
478      simnfo->sampleTime = the_globals->getSampleTime();
479      simnfo->statusTime = simnfo->sampleTime;
# Line 391 | Line 529 | void SimSetup::makeAtoms( void ){
529          current_atom = comp_stamps[i]->getAtom( k );
530          if( current_atom->haveOrientation() ){
531  
532 <          dAtom = new DirectionalAtom;
532 >          dAtom = new DirectionalAtom(index);
533            simnfo->n_oriented++;
534            the_atoms[index] = dAtom;
535  
# Line 411 | Line 549 | void SimSetup::makeAtoms( void ){
549            dAtom->setSUz( uz );
550          }
551          else{
552 <          the_atoms[index] = new GeneralAtom;
552 >          the_atoms[index] = new GeneralAtom(index);
553          }
554          the_atoms[index]->setType( current_atom->getType() );
555          the_atoms[index]->setIndex( index );
# Line 558 | Line 696 | void SimSetup::initFromBass( void ){
696      n_per_extra = (int)ceil( temp1 );
697  
698      if( n_per_extra > 4){
699 <      std::cerr << "THere has been an error in constructing the non-complete lattice.\n";
700 <      exit(8);
699 >      sprintf( painCave.errMsg,
700 >               "SimSetup error. There has been an error in constructing"
701 >               " the non-complete lattice.\n" );
702 >      painCave.isFatal = 1;
703 >      simError();
704      }
705    }
706    else{
# Line 664 | Line 805 | void SimSetup::makeElement( double x, double y, double
805  
806      current_atom = comp_stamps[current_comp]->getAtom( k );
807      if( !current_atom->havePosition() ){
808 <      std::cerr << "Component " << comp_stamps[current_comp]->getID()
809 <                << ", atom " << current_atom->getType()
810 <                << " does not have a position specified.\n"
811 <                << "The initialization routine is unable to give a start"
812 <                << " position.\n";
813 <      exit(8);
808 >      sprintf( painCave.errMsg,
809 >               "SimSetup:initFromBass error.\n"
810 >               "\tComponent %s, atom %s does not have a position specified.\n"
811 >               "\tThe initialization routine is unable to give a start"
812 >               " position.\n",
813 >               comp_stamps[current_comp]->getID(),
814 >               current_atom->getType() );
815 >      painCave.isFatal = 1;
816 >      simError();
817      }
818  
819      the_atoms[current_atom_ndx]->setX( x + current_atom->getPosX() );

Diff Legend

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