--- trunk/OOPSE/libmdtools/Integrator.cpp	2003/07/14 22:38:13	600
+++ trunk/OOPSE/libmdtools/Integrator.cpp	2003/07/15 17:57:04	614
@@ -27,8 +27,6 @@ Integrator::Integrator( SimInfo *theInfo, ForceFields*
 
   nAtoms = info->n_atoms;
 
-  std::cerr << "integ nAtoms = "  << nAtoms << "\n";
-
   // check for constraints
   
   constrainedA    = NULL;
@@ -75,12 +73,7 @@ void Integrator::checkConstraints( void ){
       
       constrained = theArray[j]->is_constrained();
 
-      std::cerr << "Is the folowing bond constrained \n";
-      theArray[j]->printMe();
-      
       if(constrained){
-	
-	std::cerr << "Yes\n";
 
 	dummy_plug = theArray[j]->get_constraint();
 	temp_con[nConstrained].set_a( dummy_plug->get_a() );
@@ -90,7 +83,6 @@ void Integrator::checkConstraints( void ){
 	nConstrained++;
 	constrained = 0;
       } 
-      else std::cerr << "No.\n";
     }
 
     theArray = (SRI**) molecules[i].getMyBends();
@@ -225,8 +217,6 @@ void Integrator::integrate( void ){
       calcStress = 1;
     }
 
-    std::cerr << currTime << "\n";
-
     integrateStep( calcPot, calcStress );
       
     currTime += dt;
@@ -274,15 +264,34 @@ void Integrator::integrateStep( int calcPot, int calcS
   moveA();
   if( nConstrained ) constrainA();
 
+  
+#ifdef IS_MPI
+  strcpy( checkPointMsg, "Succesful moveA\n" );
+  MPIcheckPoint();
+#endif // is_mpi
+  
+
   // calc forces
 
   myFF->doForces(calcPot,calcStress);
 
+#ifdef IS_MPI
+  strcpy( checkPointMsg, "Succesful doForces\n" );
+  MPIcheckPoint();
+#endif // is_mpi
+  
+
   // finish the velocity  half step
   
   moveB();
   if( nConstrained ) constrainB();
-   
+  
+#ifdef IS_MPI
+  strcpy( checkPointMsg, "Succesful moveB\n" );
+  MPIcheckPoint();
+#endif // is_mpi
+  
+ 
 }