ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/branches/development/src/primitives/DirectionalAtom.cpp
(Generate patch)

Comparing branches/development/src/primitives/DirectionalAtom.cpp (file contents):
Revision 1786 by gezelter, Fri May 18 21:44:02 2012 UTC vs.
Revision 1787 by gezelter, Wed Aug 29 18:13:11 2012 UTC

# Line 54 | Line 54 | namespace OpenMD {
54      I_ = da.getI();
55  
56      MultipoleAdapter ma = MultipoleAdapter(dAtomType);
57 <    if (ma.isMultipole()) {
58 <      electroBodyFrame_ = ma.getElectroBodyFrame();
57 >    if (ma.isDipole()) {
58 >      dipole_ = ma.getDipole();
59      }
60 +    if (ma.isQuadrupole()) {
61 +      quadrupole_ = ma.getQuadrupole();
62 +    }
63  
64      // Check if one of the diagonal inertia tensor of this directional
65      // atom is zero:
# Line 86 | Line 89 | namespace OpenMD {
89    
90    void DirectionalAtom::setPrevA(const RotMat3x3d& a) {
91      ((snapshotMan_->getPrevSnapshot())->*storage_).aMat[localIndex_] = a;
92 +
93      if (atomType_->isMultipole()) {
94 <      ((snapshotMan_->getPrevSnapshot())->*storage_).electroFrame[localIndex_] = a.transpose() * electroBodyFrame_;
94 >      RotMat3x3d atrans = a.transpose();
95 >      
96 >      if (atomType_->isDipole()) {
97 >        ((snapshotMan_->getPrevSnapshot())->*storage_).dipole[localIndex_] = atrans * dipole_;
98 >      }
99 >
100 >      if (atomType_->isQuadrupole()) {
101 >        ((snapshotMan_->getPrevSnapshot())->*storage_).quadrupole[localIndex_] = atrans * quadrupole_ * a;
102 >      }
103      }
104    }
105    
106    
107    void DirectionalAtom::setA(const RotMat3x3d& a) {
108      ((snapshotMan_->getCurrentSnapshot())->*storage_).aMat[localIndex_] = a;
109 <    
109 >
110      if (atomType_->isMultipole()) {
111 <      ((snapshotMan_->getCurrentSnapshot())->*storage_).electroFrame[localIndex_] = a.transpose() * electroBodyFrame_;
111 >      RotMat3x3d atrans = a.transpose();
112 >      
113 >      if (atomType_->isDipole()) {
114 >        ((snapshotMan_->getCurrentSnapshot())->*storage_).dipole[localIndex_] = atrans * dipole_;
115 >      }
116 >
117 >      if (atomType_->isQuadrupole()) {
118 >        ((snapshotMan_->getCurrentSnapshot())->*storage_).quadrupole[localIndex_] = atrans * quadrupole_ * a;
119 >      }
120      }
121 +  
122    }    
123    
124    void DirectionalAtom::setA(const RotMat3x3d& a, int snapshotNo) {
125      ((snapshotMan_->getSnapshot(snapshotNo))->*storage_).aMat[localIndex_] = a;
126 <    
126 >
127      if (atomType_->isMultipole()) {
128 <      ((snapshotMan_->getSnapshot(snapshotNo))->*storage_).electroFrame[localIndex_] = a.transpose() * electroBodyFrame_;    
128 >      RotMat3x3d atrans = a.transpose();
129 >      
130 >      if (atomType_->isDipole()) {
131 >        ((snapshotMan_->getSnapshot(snapshotNo))->*storage_).dipole[localIndex_] = atrans * dipole_;
132 >      }
133 >
134 >      if (atomType_->isQuadrupole()) {
135 >        ((snapshotMan_->getSnapshot(snapshotNo))->*storage_).quadrupole[localIndex_] = atrans * quadrupole_ * a;
136 >      }
137      }
138 +
139    }    
140    
141    void DirectionalAtom::rotateBy(const RotMat3x3d& m) {

Diff Legend

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