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

Timestamp:

Nov 29, 2005, 4:17:17 PM (21 years ago)

Author:

drobbin

Message:

Edited psCubeToSphere, tests. EOC - gravity still fails. Others still unchecked.

File:

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

trunk/psLib/test/astro/tst_psEarthOrientation.c

@@ -5 +5 @@
 *  @author d-Rob, MHPCC
 *
-*  @version $Revision: 1.15 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-11-17 04:04:30 $
+*  @version $Revision: 1.16 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-11-30 02:17:17 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -102 +102 @@
     psSphere *actual = psSphereAlloc();
     psSphere *apparent = NULL;
-    psSphere *sun = psSphereAlloc();
+    //    psSphere *sun = psSphereAlloc();
     psSphere *empty = NULL;
 
-    sun->r = 0.2;
-    sun->d = 0.2;
-    actual->r = 0.61001;
-    actual->d = 0.01999;
+    //    sun->r = 0.2;
+    //    sun->d = 0.2;
+    //    actual->r = 0.61001;
+    //    actual->d = 0.01999;
+    actual->r = DEG_TO_RAD(122.9153182445501);
+    actual->d = DEG_TO_RAD(48.562968978679194);
+    psCube *sunCube = psCubeAlloc();
+    sunCube->x = 1.467797790127511e11;
+    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);
@@ -123 +136 @@
     }
 
-    apparent = psGravityDeflection(apparent, actual, sun);
+    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);
+    psCube *outCube = psSphereToCube(result);
+    printf(" -- outCube = x,y,z = %.13g, %.13g, %.13g  -- \n", outCube->x, outCube->y, outCube->z);
+    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);
+
     psFree(result);
     psFree(actual);