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 1290 by cli2, Wed Sep 10 19:51:45 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),
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 < class Bond {
77 <    public:
78 <        Bond(Atom* atom1, Atom* atom2, BondType* bt) : atom1_(atom1), atom2_(atom2), bondType_(bt) {}
79 <
80 <        void calcForce() {
81 <            double len;
82 <            double dvdr;
83 <            Vector3d r12;
84 <            Vector3d force;
85 <
86 <            r12 = atom2_->getPos() - atom1_->getPos();
87 <            len = r12.length();            
88 <            bondType_->calcForce(len,  potential_, dvdr);
89 <
90 <            force = r12 * (-dvdr / len);
91 <            
92 <            atom1_->addFrc(-force);
93 <            atom2_->addFrc(force);
94 <       }
95 <        
96 <        double getPotential() {
97 <            return potential_;
98 <        }
99 <
100 <        Atom* getAtomA() {
101 <            return atom1_;
102 <        }
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 <
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