25 |
|
|
26 |
|
/** |
27 |
|
* @file PolynomialTorsionType.hpp |
28 |
< |
* @author teng lin |
28 |
> |
* @author teng lin |
29 |
|
* @date 11/16/2004 |
30 |
|
* @version 1.0 |
31 |
|
*/ |
45 |
|
|
46 |
|
public: |
47 |
|
PolynomialTorsionType() {} |
48 |
< |
PolynomialTorsionType(double r0, const std::vector<int>& power, |
49 |
< |
const std::vector<double>& coeff) : BondType(r0){ |
48 |
> |
|
49 |
> |
double getCoefficient(int power) { |
50 |
> |
return polynomial_.getCoefficient(power); |
51 |
|
} |
52 |
|
|
53 |
+ |
void addCoefficient(int power, double coefficient) { |
54 |
+ |
polynomial_.addCoefficient(power, coefficient); |
55 |
+ |
} |
56 |
+ |
|
57 |
|
void setCoefficient(int power, double coefficient) { |
58 |
|
polynomial_.setCoefficient(power, coefficient); |
59 |
|
} |
60 |
|
|
61 |
< |
double getCoefficient(int power) { |
62 |
< |
return polynomial_.getCoefficient(power); |
61 |
> |
void setPolynomial(const DoublePolynomial& p) { |
62 |
> |
polynomial_ = p; |
63 |
|
} |
64 |
|
|
65 |
|
virtual void calcForce(double cosPhi, double sinPhi, double& V, double& dVdPhi) { |