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

Comparing trunk/src/visitors/ReplacementVisitor.cpp (file contents):
Revision 1456 by gezelter, Fri Jun 25 17:40:24 2010 UTC vs.
Revision 1457 by gezelter, Fri Jun 25 18:46:47 2010 UTC

# Line 120 | Line 120 | namespace OpenMD {
120          
121      pos = datom->getPos();
122      vel = datom->getVel();
123 <    frc = datom->getFrc();
124 <    trq = datom->getTrq();
123 >
124      j   = datom->getJ();
125      I   = datom->getI();
126      A   = datom->getA();
127 +
128      skewMat(0, 0) =  0;
129      skewMat(0, 1) =  j[2] / I(2, 2);
130      skewMat(0, 2) = -j[1] / I(1, 1);    
# Line 134 | Line 134 | namespace OpenMD {
134      skewMat(2, 0) =  j[1] / I(1, 1);
135      skewMat(2, 1) = -j[0] / I(0, 0);
136      skewMat(2, 2) =  0;
137
137      Mat3x3d mat = (A * skewMat).transpose();
138      
139      // We need A^T to convert from body-fixed to space-fixed:

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines