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

Comparing trunk/OOPSE/libmdtools/ForceFields.cpp (file contents):
Revision 829 by gezelter, Tue Oct 28 16:03:37 2003 UTC vs.
Revision 841 by mmeineke, Wed Oct 29 17:55:28 2003 UTC

# Line 26 | Line 26 | void ForceFields::calcRcut( void ){
26  
27    //calc rCut and rList
28  
29 <  entry_plug->setRcut( 2.5 * bigSigma );  
29 >  entry_plug->setDefaultRcut( 2.5 * bigSigma );  
30      
31   }
32  
# Line 69 | Line 69 | void ForceFields::doForces( int calcPot, int calcStres
69      entry_plug->tau[i] = 0.0;
70    }
71  
72 +  std::cerr << "before\n"
73 +            << "  force[0] = " << frc[0] << "\n"
74 +            << "  pos[0]   = " << pos[0] << "\n"
75 +            << "  trq[0]   = " << trq[0] << "\n"
76 +            << "  A[0]     = " << A[0] << "\n"
77 +            << "  tau      = " << entry_plug->tau[0] << ", " << entry_plug->tau[1] << ", " << entry_plug->tau[2] << "\n"
78 +            << "             " << entry_plug->tau[3] << ", " << entry_plug->tau[4] << ", " << entry_plug->tau[5] << "\n"
79 +            << "             " << entry_plug->tau[6] << ", " << entry_plug->tau[7] << ", " << entry_plug->tau[8] << "\n\n";
80 +
81    fortranForceLoop( pos,
82                      A,
83                      u_l,
# Line 80 | Line 89 | void ForceFields::doForces( int calcPot, int calcStres
89                      &passedCalcStress,
90                      &isError );
91  
92 +
93 +  std::cerr << "after\n"
94 +            << "  force[0] = " << frc[0] << "\n"
95 +            << "  pos[0]   = " << pos[0] << "\n"
96 +            << "  trq[0]   = " << trq[0] << "\n"
97 +            << "  A[0]     = " << A[0] << "\n"
98 +            << "  tau      = " << entry_plug->tau[0] << ", " << entry_plug->tau[1] << ", " << entry_plug->tau[2] << "\n"
99 +            << "             " << entry_plug->tau[3] << ", " << entry_plug->tau[4] << ", " << entry_plug->tau[5] << "\n"
100 +            << "             " << entry_plug->tau[6] << ", " << entry_plug->tau[7] << ", " << entry_plug->tau[8] << "\n\n";
101 +
102    if( isError ){
103      sprintf( painCave.errMsg,
104               "Error returned from the fortran force calculation.\n" );

Diff Legend

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