Changeset 5908 for trunk/doc/pslib/eoc_testing.txt

Timestamp:

Jan 5, 2006, 2:16:50 PM (21 years ago)

Author:

Paul Price

Message:

Fixed up inverses to clarify.

File:

• trunk/doc/pslib/eoc_testing.txt (modified) (5 diffs)

trunk/doc/pslib/eoc_testing.txt

@@ -100 +100 @@
 Check that:
 
-    X = 3.05224300720406e-07
-    Y = -1.39441339235822e-07
+    X = 3.05224300720406e-10
+    Y = -1.39441339235822e-10
     S = 0
 
@@ -134 +134 @@
 (a) Generate the transformation:
 
-    psSphereRot *earthRotInv = psSphereRot_TEOtoCEO(time); // This is TEO->CEO
-    psSphereRot *earthRot = psSphereRotInvert(earthRotInv); // This is CEO->TEO
-
-Check that the rotation quaternion of earthRotInv is:
+    psSphereRot *earthRot = psSphereRot_TEOtoCEO(time); // This is TEO->CEO
+
+Check that the rotation quaternion of earthRot is:
 
     0.0
     0.0
     0.9625401009002903
-    0.2711393629830588
+    -0.2711393629830588
 
 Within psSphereRot_TEOtoCEO, the Earth rotation angle should be 428251.4641536639 degrees
@@ -148 +147 @@
 (b) Apply the rotation to the position:
 
-    obj = psSphereRotApply(obj, earthRot, obj);
+    psSphereRot *earthRotInv = psSphereRotInvert(earthRot); // This is CEO->TEO
+    obj = psSphereRotApply(obj, earthRotInv, obj);
 
 Convert obj to a psCube and check that:
@@ -172 +172 @@
     x=-6.43607313124045e-07
     y=2.11351436973568e-06
-    s=-2.05449328145735e-15
+    s=-7.39617581324646e-12
 
 (c) Generate the transformation:
 
-   psSphereRot *polarMotionInv = psSphereRot_ITRStoTEO(earthPM); // This is ITRS->TEO
-   psSphereRot *polarMotion = psSphereRotInvert(polarMotionInv); // This is TEO->ITRS
-
-Check that the rotation quaternion of polarMotionInv is:
+   psSphereRot *polarMotion = psSphereRot_ITRStoTEO(earthPM); // This is ITRS->TEO
+
+Check that the rotation quaternion of polarMotion is:
 
     -1.0567571848664005E-6
@@ -188 +187 @@
 (d) Apply these to the position:
 
-   obj = psSphereRotApply(obj, polarMotion, obj);
+   psSphereRot *polarMotionInv = psSphereRotInvert(polarMotion); // This is TEO->ITRS
+   obj = psSphereRotApply(obj, polarMotionInv, obj);
 
 Convert obj to a psCube, should get: