82 |
|
k0 = k0_; |
83 |
|
} |
84 |
|
|
85 |
< |
virtual void calcForce(double cosPhi, double sinPhi, double& V, double& dVdPhi){ |
85 |
> |
virtual void calcForce(double cosPhi, double& V, double& dVdcosPhi){ |
86 |
|
double cosPhi2 = cosPhi * cosPhi; |
87 |
|
double cosPhi3 = cosPhi2 * cosPhi; |
88 |
|
double cosPhi4 = cosPhi3 * cosPhi; |
89 |
|
|
90 |
|
V =k0_ + k1_ * cosPhi + k2_*cosPhi2 + k3_*cosPhi3 + k4_*cosPhi4; |
91 |
< |
dVdPhi = k1_ + 2.0*k2_ * cosPhi + 3.0 * k3_*cosPhi2 + 4.0*k4_*cosPhi3; |
91 |
> |
dVdcosPhi = k1_ + 2.0*k2_ * cosPhi + 3.0 * k3_*cosPhi2 + 4.0*k4_*cosPhi3; |
92 |
|
} |
93 |
|
|
94 |
|
private: |