Changeset 5306 for trunk/psLib/src/astro/psSphereOps.c

Timestamp:

Oct 13, 2005, 10:23:57 AM (21 years ago)

Author:

drobbin

Message:

Fixed psSphereRotAlloc and test. Correctly verified calculations.

File:

• trunk/psLib/src/astro/psSphereOps.c (modified) (2 diffs)

trunk/psLib/src/astro/psSphereOps.c

@@ -7 +7 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-08-30 01:14:10 $
+ *  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-10-13 20:23:57 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -153 +153 @@
 
     // following came from ADD
-    out->q0 = b3*a0 + b2*a1 - b1*a2 + b0*a3;
-    out->q1 = b3*a1 - b2*a0 + b1*a3 + b0*a2;
-    out->q2 = b3*a2 + b2*a3 + b1*a0 - b0*a1;
-    out->q3 = b3*a3 - b3*a2 - b1*a1 - b0*a0;
+    //    out->q0 = b3*a0 + b2*a1 - b1*a2 + b0*a3;
+    //    out->q1 = b3*a1 - b2*a0 + b1*a3 + b0*a2;
+    //    out->q2 = b3*a2 + b2*a3 + b1*a0 - b0*a1;
+    out->q0 = a3*b0 + a0*b3 + a1*b2 - a2*b1;
+    out->q1 = a3*b1 - a0*b2 + a1*b3 + a2*b0;
+    out->q2 = a3*b2 + a0*b1 - a1*b0 + a2*b3;
+
+    out->q3 = b3*a3 - b2*a2 - b1*a1 - b0*a0;
 
     return out;