ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE/libmdtools/Integrator.cpp
(Generate patch)

Comparing trunk/OOPSE/libmdtools/Integrator.cpp (file contents):
Revision 600 by gezelter, Mon Jul 14 22:38:13 2003 UTC vs.
Revision 637 by gezelter, Thu Jul 17 21:50:01 2003 UTC

# Line 27 | Line 27 | Integrator::Integrator( SimInfo *theInfo, ForceFields*
27  
28    nAtoms = info->n_atoms;
29  
30 –  std::cerr << "integ nAtoms = "  << nAtoms << "\n";
31 –
30    // check for constraints
31    
32    constrainedA    = NULL;
# Line 75 | Line 73 | void Integrator::checkConstraints( void ){
73        
74        constrained = theArray[j]->is_constrained();
75  
78 –      std::cerr << "Is the folowing bond constrained \n";
79 –      theArray[j]->printMe();
80 –      
76        if(constrained){
82 –        
83 –        std::cerr << "Yes\n";
77  
78          dummy_plug = theArray[j]->get_constraint();
79          temp_con[nConstrained].set_a( dummy_plug->get_a() );
# Line 90 | Line 83 | void Integrator::checkConstraints( void ){
83          nConstrained++;
84          constrained = 0;
85        }
93 –      else std::cerr << "No.\n";
86      }
87  
88      theArray = (SRI**) molecules[i].getMyBends();
# Line 225 | Line 217 | void Integrator::integrate( void ){
217        calcStress = 1;
218      }
219  
228 –    std::cerr << currTime << "\n";
229 –
220      integrateStep( calcPot, calcStress );
221        
222      currTime += dt;
223 +    info->setTime(currTime);
224  
225      if( info->setTemp ){
226        if( currTime >= currThermal ){
# Line 273 | Line 264 | void Integrator::integrateStep( int calcPot, int calcS
264    preMove();
265    moveA();
266    if( nConstrained ) constrainA();
267 +
268 +  
269 + #ifdef IS_MPI
270 +  strcpy( checkPointMsg, "Succesful moveA\n" );
271 +  MPIcheckPoint();
272 + #endif // is_mpi
273 +  
274  
275    // calc forces
276  
277    myFF->doForces(calcPot,calcStress);
278  
279 + #ifdef IS_MPI
280 +  strcpy( checkPointMsg, "Succesful doForces\n" );
281 +  MPIcheckPoint();
282 + #endif // is_mpi
283 +  
284 +
285    // finish the velocity  half step
286    
287    moveB();
288    if( nConstrained ) constrainB();
289 <  
289 >  
290 > #ifdef IS_MPI
291 >  strcpy( checkPointMsg, "Succesful moveB\n" );
292 >  MPIcheckPoint();
293 > #endif // is_mpi
294 >  
295 >
296   }
297  
298  

Diff Legend

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