リビジョン | 871cbe79b83fb500803bc276850f8cbe8417a41f (tree) |
---|---|
日時 | 2011-09-28 00:30:57 |
作者 | Mikiya Fujii <mikiya.fujii@gmai...> |
コミッター | Mikiya Fujii |
A bug in Molecule::GetCoreRepulsionFirstDerivative is fixed.
git-svn-id: https://svn.sourceforge.jp/svnroot/molds/MolDS/trunk@143 1136aad2-a195-0410-b898-f5ea1d11b9d8
@@ -209,7 +209,7 @@ double Molecule::GetCoreRepulsionFirstDerivative(int indexAtomA, int indexAtomB, | ||
209 | 209 | CartesianType axisA){ |
210 | 210 | double value=0.0; |
211 | 211 | Atom* atomA = (*this->atomVect)[indexAtomA]; |
212 | - Atom* atomB = (*this->atomVect)[indexAtomA]; | |
212 | + Atom* atomB = (*this->atomVect)[indexAtomB]; | |
213 | 213 | double distance = this->GetDistanceAtoms(indexAtomA, indexAtomB); |
214 | 214 | value = atomA->GetCoreCharge()*atomB->GetCoreCharge(); |
215 | 215 | value *= (atomA->GetXyz()[axisA] - atomB->GetXyz()[axisA])/distance; |