ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/trunk/src/primitives/Bond.hpp
(Generate patch)

Comparing trunk/src/primitives/Bond.hpp (file contents):
Revision 246 by gezelter, Wed Jan 12 22:41:40 2005 UTC vs.
Revision 1211 by gezelter, Wed Jan 23 16:38:22 2008 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 54 | Line 54 | namespace oopse {
54   #include "types/BondType.hpp"
55  
56   namespace oopse {
57 <
58 < class Bond {
59 <    public:
60 <        Bond(Atom* atom1, Atom* atom2, BondType* bt) : atom1_(atom1), atom2_(atom2), bondType_(bt) {}
61 <
62 <        void calcForce() {
63 <            double len;
64 <            double dvdr;
65 <            Vector3d r12;
66 <            Vector3d force;
67 <
68 <            r12 = atom2_->getPos() - atom1_->getPos();
69 <            len = r12.length();            
70 <            bondType_->calcForce(len,  potential_, dvdr);
71 <
72 <            force = r12 * (-dvdr / len);
73 <            
74 <            atom1_->addFrc(-force);
75 <            atom2_->addFrc(force);
76 <       }
77 <        
78 <        double getPotential() {
79 <            return potential_;
80 <        }
81 <
82 <        Atom* getAtomA() {
83 <            return atom1_;
84 <        }
85 <
86 <        Atom* getAtomB() {
87 <            return atom2_;
88 <        }
89 <
90 <        BondType* getBondType() {
91 <            return bondType_;
92 <        }
93 <        
94 <    private:
95 <        double potential_;
96 <        Atom* atom1_;
97 <        Atom* atom2_;        
98 <        BondType* bondType_; /**< bond type */
99 <
100 < };    
101 <
102 <
103 <
57 >  
58 >  class Bond {
59 >  public:
60 >    Bond(Atom* atom1, Atom* atom2, BondType* bt) : atom1_(atom1),
61 >                                                   atom2_(atom2),
62 >                                                   bondType_(bt) {}
63 >    virtual ~Bond() {}
64 >    void calcForce() {
65 >      RealType len;
66 >      RealType dvdr;
67 >      Vector3d r12;
68 >      Vector3d force;
69 >      
70 >      r12 = atom2_->getPos() - atom1_->getPos();
71 >      len = r12.length();            
72 >      bondType_->calcForce(len,  potential_, dvdr);
73 >      
74 >      force = r12 * (-dvdr / len);
75 >      
76 >      atom1_->addFrc(-force);
77 >      atom2_->addFrc(force);
78 >    }
79 >    
80 >    RealType getPotential() {
81 >      return potential_;
82 >    }
83 >    
84 >    Atom* getAtomA() {
85 >      return atom1_;
86 >    }
87 >    
88 >    Atom* getAtomB() {
89 >      return atom2_;
90 >    }
91 >    
92 >    BondType* getBondType() {
93 >      return bondType_;
94 >    }
95 >    
96 >  private:
97 >    RealType potential_;
98 >    Atom* atom1_;
99 >    Atom* atom2_;        
100 >    BondType* bondType_; /**< bond type */
101 >    
102 >  };    
103   } //end namespace oopse
104   #endif //PRIMITIVES_BOND_HPP

Diff Legend

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