# | Line 1 | Line 1 | |
---|---|---|
1 | #include <cmath> | |
2 | ||
3 | #include "SRI.hpp" | |
4 | + | #include "simError.h" |
5 | ||
6 | QuadraticBend::QuadraticBend( Atom &a, Atom &b, Atom &c ){ | |
7 | ||
# | Line 25 | Line 26 | double QuadraticBend::bend_force( double theta ){ | |
26 | double dt, dt2; | |
27 | double force; | |
28 | ||
29 | + | |
30 | + | |
31 | + | |
32 | dt = ( theta - theta0 ) * M_PI / 180.0; | |
33 | dt2 = dt * dt; | |
34 |
– | Removed lines |
+ | Added lines |
< | Changed lines |
> | Changed lines |