61 |
|
|
62 |
|
public: |
63 |
|
|
64 |
< |
UreyBradleyBendType(double theta0, double ktheta, double s0, double kub) |
64 |
> |
UreyBradleyBendType(RealType theta0, RealType ktheta, RealType s0, RealType kub) |
65 |
|
: HarmonicBendType(theta0, ktheta), hbt_(HarmonicBondType(s0, kub)){ |
66 |
|
} |
67 |
|
|
68 |
< |
double getUBForceConstant() { |
68 |
> |
RealType getUBForceConstant() { |
69 |
|
return hbt_.getForceConstant(); |
70 |
|
} |
71 |
|
|
72 |
< |
void setUBForceConstant(double kub) { |
72 |
> |
void setUBForceConstant(RealType kub) { |
73 |
|
hbt_.setForceConstant(kub); |
74 |
|
} |
75 |
|
|
76 |
< |
double getUBLength() { |
76 |
> |
RealType getUBLength() { |
77 |
|
return hbt_.getEquilibriumBondLength(); |
78 |
|
} |
79 |
|
|
80 |
< |
void setUBLength(double s) { |
80 |
> |
void setUBLength(RealType s) { |
81 |
|
hbt_.setEquilibriumBondLength(s); |
82 |
|
} |
83 |
|
|
87 |
|
|
88 |
|
private: |
89 |
|
|
90 |
< |
double k_; |
90 |
> |
RealType k_; |
91 |
|
HarmonicBondType hbt_; |
92 |
|
}; |
93 |
|
|