322 |
|
transposeMat3(pAxisMat, pAxisRotMat); |
323 |
|
|
324 |
|
|
325 |
– |
for (i=0; i<myAtoms.size(); i++){ |
326 |
– |
apos = refCoords[i]; |
327 |
– |
printf("%f\t%f\t%f\n",apos[0],apos[1],apos[2]); |
328 |
– |
} |
329 |
– |
|
325 |
|
//rotate the rigid body to the principle axis frame |
326 |
|
for (i = 0; i < myAtoms.size(); i++) { |
327 |
|
matVecMul3(pAxisRotMat, refCoords[i].vec, refCoords[i].vec); |
328 |
|
myAtoms[i]->setPos(refCoords[i].vec); |
334 |
– |
} |
335 |
– |
|
336 |
– |
for (i=0; i<myAtoms.size(); i++){ |
337 |
– |
apos = refCoords[i]; |
338 |
– |
printf("%f\t%f\t%f\n",apos[0],apos[1],apos[2]); |
329 |
|
} |
330 |
|
|
331 |
|
identityMat3(iMat); |
536 |
|
return myAtoms[index]->getEps(); |
537 |
|
|
538 |
|
} |
539 |
+ |
|
540 |
+ |
char *RigidBody::getAtomType(int index){ |
541 |
+ |
|
542 |
+ |
return myAtoms[index]->getType(); |
543 |
+ |
|
544 |
+ |
} |