--- trunk/OOPSE/libmdtools/BendExtensions.cpp 2003/03/28 22:34:02 437 +++ trunk/OOPSE/libmdtools/BendExtensions.cpp 2003/03/31 21:50:59 438 @@ -1,6 +1,7 @@ #include #include "SRI.hpp" +#include "simError.h" QuadraticBend::QuadraticBend( Atom &a, Atom &b, Atom &c ){ @@ -25,6 +26,9 @@ double QuadraticBend::bend_force( double theta ){ double dt, dt2; double force; + + + dt = ( theta - theta0 ) * M_PI / 180.0; dt2 = dt * dt;