--- trunk/OOPSE/libmdtools/ForceFields.cpp	2003/10/28 16:03:37	829
+++ trunk/OOPSE/libmdtools/ForceFields.cpp	2003/10/29 17:55:28	841
@@ -26,7 +26,7 @@ void ForceFields::calcRcut( void ){
 
   //calc rCut and rList
 
-  entry_plug->setRcut( 2.5 * bigSigma );  
+  entry_plug->setDefaultRcut( 2.5 * bigSigma );  
     
 }
 
@@ -69,6 +69,15 @@ void ForceFields::doForces( int calcPot, int calcStres
     entry_plug->tau[i] = 0.0;
   }
 
+  std::cerr << "before\n"
+	    << "  force[0] = " << frc[0] << "\n"
+	    << "  pos[0]   = " << pos[0] << "\n"
+	    << "  trq[0]   = " << trq[0] << "\n"
+	    << "  A[0]     = " << A[0] << "\n"
+	    << "  tau      = " << entry_plug->tau[0] << ", " << entry_plug->tau[1] << ", " << entry_plug->tau[2] << "\n"
+	    << "             " << entry_plug->tau[3] << ", " << entry_plug->tau[4] << ", " << entry_plug->tau[5] << "\n"
+	    << "             " << entry_plug->tau[6] << ", " << entry_plug->tau[7] << ", " << entry_plug->tau[8] << "\n\n"; 
+
   fortranForceLoop( pos,
 		    A,
 		    u_l,
@@ -80,6 +89,16 @@ void ForceFields::doForces( int calcPot, int calcStres
 		    &passedCalcStress,
 		    &isError );
 
+
+  std::cerr << "after\n"
+	    << "  force[0] = " << frc[0] << "\n"
+	    << "  pos[0]   = " << pos[0] << "\n"
+	    << "  trq[0]   = " << trq[0] << "\n"
+	    << "  A[0]     = " << A[0] << "\n"
+	    << "  tau      = " << entry_plug->tau[0] << ", " << entry_plug->tau[1] << ", " << entry_plug->tau[2] << "\n"
+	    << "             " << entry_plug->tau[3] << ", " << entry_plug->tau[4] << ", " << entry_plug->tau[5] << "\n"
+	    << "             " << entry_plug->tau[6] << ", " << entry_plug->tau[7] << ", " << entry_plug->tau[8] << "\n\n";
+
   if( isError ){
     sprintf( painCave.errMsg,
 	     "Error returned from the fortran force calculation.\n" );