ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/SHAPES/RigidBody.cpp
(Generate patch)

Comparing trunk/SHAPES/RigidBody.cpp (file contents):
Revision 1281 by chrisfen, Mon Jun 21 13:38:55 2004 UTC vs.
Revision 1283 by gezelter, Mon Jun 21 15:54:27 2004 UTC

# Line 48 | Line 48 | void RigidBody::setEuler( double phi, double theta, do
48      A[2][1] = -cos(phi) * sin(theta);
49      A[2][2] = cos(theta);
50  
51 +    printf("A[2][x] = %lf\t%lf\t%lf\n", A[2][0], A[2][1], A[2][2]);
52 +
53   }
54  
55   void RigidBody::getQ( double q[4] ){
# Line 321 | Line 323 | void RigidBody::calcRefCoords( ) {
323    //calculate the proper rotation matrix
324    transposeMat3(pAxisMat, pAxisRotMat);
325    
326 +  
327    for (i=0; i<myAtoms.size(); i++){
328 <    getAtomPos(test, i);
329 <    printf("%d\t%d\t%d\n",test[0],test[1],test[2]);
328 >    apos = refCoords[i];
329 >    printf("%f\t%f\t%f\n",apos[0],apos[1],apos[2]);
330    }
331      
332    //rotate the rigid body to the principle axis frame
# Line 333 | Line 336 | void RigidBody::calcRefCoords( ) {
336    }
337    
338    for (i=0; i<myAtoms.size(); i++){
339 <    getAtomPos(test,i);
340 <    printf("%d\t%d\t%d\n",test[0],test[1],test[2]);
339 >    apos = refCoords[i];
340 >    printf("%f\t%f\t%f\n",apos[0],apos[1],apos[2]);
341    }
342    
343    identityMat3(iMat);

Diff Legend

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