[ViewVC] Diff of: group/trunk/OOPSE/utils/sysbuilder/MoLocator.cpp

Comparing trunk/OOPSE/utils/sysbuilder/MoLocator.cpp (file contents):
Revision 700 by chuckv, Mon Aug 18 20:59:47 2003 UTC vs.
Revision 821 by mmeineke, Fri Oct 24 22:17:45 2003 UTC

-+
-+
+void getRandomRot( double rot[3][3] ){
-+
-+
+  double theta, phi, psi;
-+
+  double cosTheta;
-+
-+
+  // select random phi, psi, and cosTheta
-+
-+
+  phi = 2.0 * M_PI * drand48();
-+
+  psi = 2.0 * M_PI * drand48();
-+
+  cosTheta = (2.0 * drand48()) - 1.0; // sample cos -1 to 1
-+
-+
+  theta = acos( cosTheta );
-+
-+
+  getEulerRot( theta, phi, psi, rot );
-+
+}
-+
-+
-+
+void getEulerRot( double theta, double phi, double psi, double rot[3][3] ){
-+
-+
+  rot[0][0] = (cos(phi) * cos(psi)) - (sin(phi) * cos(theta) * sin(psi));
-+
+  rot[0][1] = (sin(phi) * cos(psi)) + (cos(phi) * cos(theta) * sin(psi));
-+
+  rot[0][2] = sin(theta) * sin(psi);
-+
-+
+  rot[1][0] = -(cos(phi) * sin(psi)) - (sin(phi) * cos(theta) * cos(psi));
-+
+  rot[1][1] = -(sin(phi) * sin(psi)) + (cos(phi) * cos(theta) * cos(psi));
-+
+  rot[1][2] = sin(theta) * cos(psi);
-+
-+
+  rot[2][0] = sin(phi) * sin(theta);
-+
+  rot[2][1] = -cos(phi) * sin(theta);
-+
+  rot[2][2] = cos(theta);
-+
+}
-+
-+
-+
+void getUnitRot( double u[3], double rot[3][3] ){
-+
-+
+  double theta, phi, psi;
-+
-+
+  theta = acos(u[2]);
-+
+  phi = atan(u[1] / u[0]);
-+
+  psi = 0.0;
-+
-+
+  getEulerRot( theta, phi, psi, rot );
-+
+}
-+
-+

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