--- trunk/OOPSE-4/src/types/CubicBendType.hpp 2005/01/12 22:41:40 1930 +++ trunk/OOPSE-4/src/types/CubicBendType.hpp 2005/04/15 22:04:00 2204 @@ -1,4 +1,4 @@ - /* +/* * Copyright (c) 2005 The University of Notre Dame. All Rights Reserved. * * The University of Notre Dame grants you ("Licensee") a @@ -52,49 +52,49 @@ namespace oopse { #include "types/BendType.hpp" namespace oopse { -/** - * @class CubicBendType - * @todo document - */ -class CubicBendType : public BendType { + /** + * @class CubicBendType + * @todo document + */ + class CubicBendType : public BendType { - public: + public: - CubicBendType(double theta, double k3, double k2, double k1, double k0) - : BendType(theta),k3_(k3), k2_(k2), k1_(k1), k0_(k0){ - } + CubicBendType(double theta, double k3, double k2, double k1, double k0) + : BendType(theta),k3_(k3), k2_(k2), k1_(k1), k0_(k0){ + } - void setForceConstant(double k3, double k2, double k1, double k0) { - k3_ = k3; - k2_ = k2; - k1_ = k1; - k0_ = k0; - } + void setForceConstant(double k3, double k2, double k1, double k0) { + k3_ = k3; + k2_ = k2; + k1_ = k1; + k0_ = k0; + } - void getForceConstant(double& k3, double& k2, double& k1, double& k0) { - k3 = k3_; - k2 = k2_; - k1 = k1_; - k0 = k0_; - } + void getForceConstant(double& k3, double& k2, double& k1, double& k0) { + k3 = k3_; + k2 = k2_; + k1 = k1_; + k0 = k0_; + } - virtual void calcForce(double theta, double& V, double& dVdTheta) { - double delta = theta- theta0_; - double delta2 = delta * delta; - double delta3 = delta2 * delta; + virtual void calcForce(double theta, double& V, double& dVdTheta) { + double delta = theta- theta0_; + double delta2 = delta * delta; + double delta3 = delta2 * delta; - V =k0_ + k1_ * delta + k2_*delta2 + k3_*delta3; - dVdTheta = k1_ + 2.0*k2_ * delta + 3.0 * k3_*delta2; - } + V =k0_ + k1_ * delta + k2_*delta2 + k3_*delta3; + dVdTheta = k1_ + 2.0*k2_ * delta + 3.0 * k3_*delta2; + } - private: + private: - double k3_; - double k2_; - double k1_; - double k0_; + double k3_; + double k2_; + double k1_; + double k0_; -}; + }; }//end namespace oopse #endif //TYPES_CUBICBENDTYPE_HPP