24 |
|
*/ |
25 |
|
|
26 |
|
/** |
27 |
< |
* @file BendType.hpp |
27 |
> |
* @file QuadraticBendType.hpp |
28 |
|
* @author tlin |
29 |
|
* @date 11/01/2004 |
30 |
|
* @version 1.0 |
60 |
|
k4 = k4_; |
61 |
|
} |
62 |
|
|
63 |
< |
virtual void calcForce(double theta, double& V, double& dVdelta) { |
63 |
> |
virtual void calcForce(double theta, double& V, double& dVdTheta) { |
64 |
|
double delta = theta- theta0_; |
65 |
|
double delta2 = delta * delta; |
66 |
|
double delta3 = delta2 * delta; |
67 |
|
double delta4 = delta3 * delta; |
68 |
|
|
69 |
|
V = 0.5 * (k2_*delta2 + k3_*delta3 + k4_*delta4); |
70 |
< |
dVdelta = k2_ * delta + 1.5 * k3_*delta2 + 2.0 * k4_ *delta3; |
70 |
> |
dVdTheta = k2_ * delta + 1.5 * k3_*delta2 + 2.0 * k4_ *delta3; |
71 |
|
} |
72 |
|
|
73 |
|
private: |