ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE-2.0/src/types/QuarticBondType.hpp
(Generate patch)

Comparing trunk/OOPSE-2.0/src/types/QuarticBondType.hpp (file contents):
Revision 1930 by gezelter, Wed Jan 12 22:41:40 2005 UTC vs.
Revision 2204 by gezelter, Fri Apr 15 22:04:00 2005 UTC

# Line 1 | Line 1
1 < /*
1 > /*
2   * Copyright (c) 2005 The University of Notre Dame. All Rights Reserved.
3   *
4   * The University of Notre Dame grants you ("Licensee") a
# Line 52 | Line 52 | namespace oopse {
52   #include "types/BondType.hpp"
53  
54   namespace oopse {
55 < /**
56 < * @class QuarticBondType
57 < * @todo document
58 < */
59 < class QuarticBondType : public BondType {
55 >  /**
56 >   * @class QuarticBondType
57 >   * @todo document
58 >   */
59 >  class QuarticBondType : public BondType {
60  
61 <    public:
61 >  public:
62          
63          
64 <        QuarticBondType(double r0, double k4, double k3, double k2, double k1, double k0)
65 <            : BondType(r0), k4_(k4), k3_(k3), k2_(k2),  k1_(k1), k0_(k0){
66 <        }
64 >    QuarticBondType(double r0, double k4, double k3, double k2, double k1, double k0)
65 >      : BondType(r0), k4_(k4), k3_(k3), k2_(k2),  k1_(k1), k0_(k0){
66 >      }
67  
68 <        void setForceConstant(double k4, double k3, double k2, double k1, double k0) {
69 <            k4_ = k4;
70 <            k3_ = k3;
71 <            k2_ = k2;
72 <            k1_ = k1;
73 <            k0_ = k0;
68 >    void setForceConstant(double k4, double k3, double k2, double k1, double k0) {
69 >      k4_ = k4;
70 >      k3_ = k3;
71 >      k2_ = k2;
72 >      k1_ = k1;
73 >      k0_ = k0;
74  
75 <        }
75 >    }
76  
77 <        void getForceConstant(double& k4, double& k3, double& k2, double& k1, double& k0) {
78 <            k4 = k4_;
79 <            k3 = k3_;
80 <            k2  = k2_;
81 <            k1 = k1_;
82 <            k0 = k0_;
83 <        }
77 >    void getForceConstant(double& k4, double& k3, double& k2, double& k1, double& k0) {
78 >      k4 = k4_;
79 >      k3 = k3_;
80 >      k2  = k2_;
81 >      k1 = k1_;
82 >      k0 = k0_;
83 >    }
84  
85 <        virtual void calcForce(double r, double& V, double& dVdr) {
86 <            double dr =  r- r0;
87 <            double dr2 = dr * dr;
88 <            double dr3 = dr2 * dr;
89 <            double dr4 = dr3 * dr;
85 >    virtual void calcForce(double r, double& V, double& dVdr) {
86 >      double dr =  r- r0;
87 >      double dr2 = dr * dr;
88 >      double dr3 = dr2 * dr;
89 >      double dr4 = dr3 * dr;
90              
91 <            V =k0_ + k1_ * dr + k2_*dr2 + k3_*dr3 + k4_*dr4;
92 <            dVdr = k1_ + 2.0*k2_ * dr + 3.0 * k3_*dr2 + 4.0*k4_*dr3;            
93 <        }        
91 >      V =k0_ + k1_ * dr + k2_*dr2 + k3_*dr3 + k4_*dr4;
92 >      dVdr = k1_ + 2.0*k2_ * dr + 3.0 * k3_*dr2 + 4.0*k4_*dr3;            
93 >    }        
94          
95 <    private:
96 <        double k4_;
97 <        double k3_;
98 <        double k2_;
99 <        double k1_;
100 <        double k0_;
95 >  private:
96 >    double k4_;
97 >    double k3_;
98 >    double k2_;
99 >    double k1_;
100 >    double k0_;
101  
102 < };
102 >  };
103  
104   }//end namespace oopse
105   #endif //TYPES_QUADRATICBONDTYPE_HPP

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines