Changeset 5657 for trunk/psLib/test/astro/tst_psEarthOrientation.c

Timestamp:

Dec 1, 2005, 6:40:00 PM (21 years ago)

Author:

drobbin

Message:

Debugging/Test code updates. Aberration, Gravity, and CEOtoTEO mainly

File:

• trunk/psLib/test/astro/tst_psEarthOrientation.c (modified) (7 diffs)

trunk/psLib/test/astro/tst_psEarthOrientation.c

@@ -5 +5 @@
 *  @author d-Rob, MHPCC
 *
-*  @version $Revision: 1.17 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-11-30 23:50:40 $
+*  @version $Revision: 1.18 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-12-02 04:40:00 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -56 +56 @@
 
 #define timesec 1049160600
+#define objR DEG_TO_RAD(122.9153182445501)
+#define objD DEG_TO_RAD(48.562968978679194)
+static psSphere *obj = NULL;
+static void objSetup(void);
+
+static void objSetup(void)
+{
+    if (obj == NULL) {
+        obj = psSphereAlloc();
+        obj->r = objR;
+        obj->d = objD;
+    }
+}
 
 psS32 testAberration(void)
@@ -86 +99 @@
     double speed = sqrt(cubeDir->x*cubeDir->x + cubeDir->y*cubeDir->y + cubeDir->z*cubeDir->z);
     // Speed of light in vacuum (src:NIST)
-    double c;
-    c = 299792458.0; /* m/s */
+    double c = 299792458.0; /* m/s */
     speed = speed / c;
 
@@ -104 +116 @@
 
     apparent = psAberration(apparent, actual, direction, speed);
-    if (apparent == NULL)
-        printf("\nApparent is NULL\n");
-    else {
-        printf("\napparent = r,d  = %.8g, %.8g\n", apparent->r, apparent->d);
-        psFree(cubeDir);
-        cubeDir = psSphereToCube(apparent);
-        printf("  ><><Apparent = x,y,z = %.13g, %.13g, %.13g \n", cubeDir->x, cubeDir->y, cubeDir->z);
-    }
+
+    printf("\nSphere Difference  =  r,d  =  %.13g, %.13g\n",
+           (actual->r - apparent->r), (actual->d - apparent->d));
+    psCube *outCube = psSphereToCube(apparent);
+    printf(" -- resultCube = x,y,z = %.13g, %.13g, %.13g  -- \n",
+           outCube->x, outCube->y, outCube->z);
+
+    double x, y, z;
+    x = -0.35963388069046304;
+    y = 0.5555192509816625;
+    z = 0.7497078321908413;
+
+    printf(" -- expectedCube = x,y,z = %.13g, %.13g, %.13g  -- \n\n", x, y, z);
+
+    if ( fabs(x - outCube->x) > DBL_EPSILON || fabs(y - outCube->y) > DBL_EPSILON ||
+            fabs(z - outCube->z) > DBL_EPSILON ) {
+        psError(PS_ERR_BAD_PARAMETER_VALUE, false,
+                "psAberration returned incorrect values.\n");
+        printf("Cube Difference  =  x,y,z  = %.13g, %.13g, %.13g \n\n",
+               (x - outCube->x), (y - outCube->y), (z - outCube->z) );
+        return 3;
+    }
+
+    psFree(outCube);
+
     psFree(cubeDir);
     psFree(apparent);
@@ -138 +167 @@
     sunCube->y = 2.5880956908748722e10;
     sunCube->z = 1.1220046291457653e10;
-    //    double sum = sqrt(sunCube->x*sunCube->x + sunCube->y*sunCube->y + sunCube->z*sunCube->z);
-    //    sunCube->x = sunCube->x/sum;
-    //    sunCube->y = sunCube->y/sum;
-    //    sunCube->z = sunCube->z/sum;
-    //    printf("\n <<<<SUNvec = x,y,z = %.13g, %.13g %.13g", sunCube->x, sunCube->y, sunCube->z);
     psSphere *sun = psCubeToSphere(sunCube);
-    printf("\n <<<<SUN = r,d = %.13g, %.13g\n", sun->r, sun->d);
 
     empty = psGravityDeflection(apparent, empty, sun);
@@ -161 +184 @@
     apparent = psGravityDeflection(NULL, actual, sun);
     psSphere *result = psSphereSetOffset(actual, apparent, PS_SPHERICAL, PS_RADIAN);
+
     printf("\nActual r,d = %.13g,%.13g    Apparent r,d = %.16g, %.16g \n",
            actual->r, actual->d, result->r, result->d);
+    printf("Sphere Difference  =  r,d  =  %.13g, %.13g\n",
+           (actual->r - result->r), (actual->d - result->d));
     psCube *outCube = psSphereToCube(result);
-    printf(" -- outCube = x,y,z = %.13g, %.13g, %.13g  -- \n", outCube->x, outCube->y, outCube->z);
+    printf(" -- resultCube = x,y,z = %.13g, %.13g, %.13g  -- \n",
+           outCube->x, outCube->y, outCube->z);
+    psCube *outCube2 = psSphereToCube(actual);
+    printf(" -- actualCube = x,y,z = %.13g, %.13g, %.13g  -- \n",
+           outCube2->x, outCube2->y, outCube2->z);
+
+    double x, y, z;
+    x = -0.35961949760293604;
+    y = 0.5555613950298085;
+    z = 0.7496835020836093;
+
+    printf(" -- expectedCube = x,y,z = %.13g, %.13g, %.13g  -- \n\n", x, y, z);
+
+    if ( fabs(x - outCube->x) > DBL_EPSILON || fabs(y - outCube->y) > DBL_EPSILON ||
+            fabs(z - outCube->z) > DBL_EPSILON ) {
+        psError(PS_ERR_BAD_PARAMETER_VALUE, false,
+                "psGravityDeflection returned incorrect values.\n");
+        printf("Cube Difference  =  x,y,z  = %.13g, %.13g, %.13g \n",
+               (x - outCube->x), (y - outCube->y), (z - outCube->z) );
+        return 1;
+    }
+
     psFree(outCube);
-    psCube *outCube2 = psSphereToCube(apparent);
-    printf(" -- outCube = x,y,z = %.13g, %.13g, %.13g  -- \n", outCube2->x, outCube2->y, outCube2->z);
     psFree(outCube2);
     psFree(sunCube);
@@ -416 +461 @@
 psS32 testSphereRot_TEOtoCEO(void)
 {
+    psSphereRot *rot = NULL;
+    psTime *empty = NULL;
     //    psTime *now = psTimeAlloc(PS_TIME_UT1);
     //    now->sec = 1128530000;
     //    now->nsec = 931154510;
-    psTime *time = psTimeAlloc(PS_TIME_UTC);
-    time->sec = timesec;
-    time->nsec = 0;
+    psTime *time = psTimeAlloc(PS_TIME_UT1);
+    time->sec = timesec-1;
+    time->nsec = 657017200;
     time->leapsecond = false;
+
     //    psSphereRot *teoceo = psSphereRot_TEOtoCEO(now);
+
+    //return NULL for NULL input time
+    rot = psSphereRot_TEOtoCEO(empty);
+    if (rot != NULL) {
+        psError(PS_ERR_BAD_PARAMETER_VALUE, false,
+                "psSphereRot_TEOtoCEO failed to return NULL for NULL input time.\n");
+        return 1;
+    }
+
     psSphereRot *teoceo = psSphereRot_TEOtoCEO(time);
-    printf("Output sphere rotation = %lf,%lf,%lf,%lf\n",
+    //Make sure values match for other psTime type
+    empty = psTimeAlloc(PS_TIME_UTC);
+    empty->sec = timesec;
+    empty->nsec = 0;
+    empty->leapsecond = false;
+    rot = psSphereRot_TEOtoCEO(empty);
+    if (fabs(rot->q0-teoceo->q0) > DBL_EPSILON || fabs(rot->q1-teoceo->q1) > DBL_EPSILON ||
+            fabs(rot->q2-teoceo->q2) > DBL_EPSILON || fabs(rot->q3-teoceo->q3) > DBL_EPSILON) {
+        psError(PS_ERR_BAD_PARAMETER_VALUE, false,
+                "psSphereRot_TEOtoCEO failed to return matching values for different time types.\n");
+        return 2;
+    }
+    printf("\nOutput sphere rotation = %.13g, %.13g, %.13g, %.13g\n\n",
            teoceo->q0, teoceo->q1, teoceo->q2, teoceo->q3);
 
+    objSetup();
+    psSphereRot *earthRot = psSphereRotConjugate(NULL, teoceo);
+    psSphere *result = psSphereRotApply(NULL, earthRot, obj);
+    //    psCube *cube = psSphereToCube(obj);
+    psCube *cube = psSphereToCube(result);
+
+
+    double x, y, z;
+    x = 0.01698625430807123;
+    y = -0.6616523084626379;
+    z = 0.7496183628158023;
+    if ( fabs(x-cube->x) > DBL_EPSILON  || fabs(y-cube->y) > DBL_EPSILON ||
+            fabs(z-cube->z) > DBL_EPSILON) {
+        psError(PS_ERR_BAD_PARAMETER_VALUE, false,
+                "psSphereRot_TEOtoCEO returned incorrect values.\n");
+        printf("\nOutput cube = x,y,z = %.13g, %.13g, %.13g\n",
+               cube->x, cube->y, cube->z);
+        printf("Expected cube = x,y,z = %.13g, %.13g, %.13g\n", x, y, z);
+        printf("A difference of:   %.13g, %.13g, %.13g\n\n",
+               (x-cube->x), (y-cube->y), (z-cube->z));
+        return 3;
+    }
     //    psFree(now);
+    psFree(result);
+    psFree(obj);
+    psFree(earthRot);
+    psFree(cube);
     psFree(time);
     psFree(teoceo);