Index: /trunk/psLib/src/astro/psEarthOrientation.c =================================================================== --- /trunk/psLib/src/astro/psEarthOrientation.c (revision 7117) +++ /trunk/psLib/src/astro/psEarthOrientation.c (revision 7118) @@ -8,6 +8,6 @@ * @author Robert Daniel DeSonia, MHPCC * - * @version $Revision: 1.35 $ $Name: not supported by cvs2svn $ - * @date $Date: 2006-02-14 00:09:26 $ + * @version $Revision: 1.36 $ $Name: not supported by cvs2svn $ + * @date $Date: 2006-05-16 01:08:57 $ * * Copyright 2005 Maui High Performance Computing Center, University of Hawaii @@ -1018,5 +1018,5 @@ // Convert psTime to MJD double MJD = psTimeToMJD(time); - printf("\nMJD check = %.13g\n", MJD); + // printf("\nMJD check = %.13g\n", MJD); if (MJD == NAN) { psError(PS_ERR_BAD_PARAMETER_VALUE, false, Index: /trunk/psLib/test/astro/tst_psEarthOrientation.c =================================================================== --- /trunk/psLib/test/astro/tst_psEarthOrientation.c (revision 7117) +++ /trunk/psLib/test/astro/tst_psEarthOrientation.c (revision 7118) @@ -5,6 +5,6 @@ * @author d-Rob, MHPCC * -* @version $Revision: 1.31 $ $Name: not supported by cvs2svn $ -* @date $Date: 2006-04-26 01:09:13 $ +* @version $Revision: 1.32 $ $Name: not supported by cvs2svn $ +* @date $Date: 2006-05-16 01:08:58 $ * * Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii @@ -64,4 +64,5 @@ #define objR DEG_TO_RAD(122.9153182445501) #define objD DEG_TO_RAD(48.562968978679194) +#define VERBOSE 0 static psSphere *obj = NULL; static void objSetup(void); @@ -149,6 +150,8 @@ apparent = psAberration(apparent, actual, direction, speed); psCube *outCube = psSphereToCube(apparent); - printf("\n -- resultCube = x,y,z = %.13g, %.13g, %.13g -- \n", - outCube->x, outCube->y, outCube->z); + if (VERBOSE) { + printf("\n -- resultCube = x,y,z = %.13g, %.13g, %.13g -- \n", + outCube->x, outCube->y, outCube->z); + } //expected cube values double x, y, z; @@ -157,7 +160,9 @@ z = 0.7497078321908413; - printf(" -- expectedCube = x,y,z = %.13g, %.13g, %.13g -- \n", x, y, z); - printf("Cube Difference = x,y,z = %.13g, %.13g, %.13g \n\n", - (x - outCube->x), (y - outCube->y), (z - outCube->z) ); + if (VERBOSE) { + printf(" -- expectedCube = x,y,z = %.13g, %.13g, %.13g -- \n", x, y, z); + printf("Cube Difference = x,y,z = %.13g, %.13g, %.13g \n\n", + (x - outCube->x), (y - outCube->y), (z - outCube->z) ); + } psFree(actual); actualCube->x = x; @@ -166,7 +171,8 @@ actual = psCubeToSphere(actualCube); double xxx = greatCircle(actual, apparent); - printf(" The great circle angular distance between expected & result = %.13g\n", - xxx); - + if (VERBOSE) { + printf(" The great circle angular distance between expected & result = %.13g\n", + xxx); + } if ( fabs(x - outCube->x) > FLT_EPSILON || fabs(y - outCube->y) > FLT_EPSILON || fabs(z - outCube->z) > FLT_EPSILON ) { @@ -184,6 +190,8 @@ actual = psCubeToSphere(actualCube); x = greatCircle(actual, apparent); - printf(" The great circle angular distance between expected & result = %.13g\n", - x); + if (VERBOSE) { + printf(" The great circle angular distance between expected & result = %.13g\n", + x); + } return 3; } @@ -216,7 +224,12 @@ sunCube->y /= sunLength; sunCube->z /= sunLength; - printf("sunCube = x,y,z = %.13g, %.13g, %.13g\n", sunCube->x, sunCube->y, sunCube->z); + if (VERBOSE) { + printf("sunCube = x,y,z = %.13g, %.13g, %.13g\n", + sunCube->x, sunCube->y, sunCube->z); + } psSphere *sun = psCubeToSphere(sunCube); - printf("sunSphere = r, d = %.13g, %.13g\n", sun->r, sun->d); + if (VERBOSE) { + printf("sunSphere = r, d = %.13g, %.13g\n", sun->r, sun->d); + } psCube *outCube = psCubeAlloc(); @@ -241,12 +254,16 @@ // psSphere *result = psSphereSetOffset(actual, apparent, PS_SPHERICAL, PS_RADIAN); // psSphere *result = psSphereGetOffset(apparent, actual, PS_SPHERICAL, PS_RADIAN); - printf(" -- actualCube = x,y,z = %.13g, %.13g, %.13g -- \n", - actualCube->x, actualCube->y, actualCube->z); + if (VERBOSE) { + printf(" -- actualCube = x,y,z = %.13g, %.13g, %.13g -- \n", + actualCube->x, actualCube->y, actualCube->z); + } // psCube *outCube = psSphereToCube(result); // printf(" -- resultCube = x,y,z = %.13g, %.13g, %.13g -- \n", // outCube->x, outCube->y, outCube->z); psCube *outCube2 = psSphereToCube(apparent); - printf(" -- resultCube2= x,y,z = %.13g, %.13g, %.13g -- \n", - outCube2->x, outCube2->y, outCube2->z); + if (VERBOSE) { + printf(" -- resultCube2= x,y,z = %.13g, %.13g, %.13g -- \n", + outCube2->x, outCube2->y, outCube2->z); + } double x, y, z; x = -0.35961949760293604; @@ -254,17 +271,19 @@ z = 0.7496835020836093; - printf(" -- expectCube = 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("expect-actual= x,y,z = %.13g, %.13g, %.13g \n", - (x - actualCube->x), (y - actualCube->y), (z - actualCube->z) ); - printf("expect-result= x,y,z = %.13g, %.13g, %.13g \n", - (x - outCube2->x), (y - outCube2->y), (z - outCube2->z) ); - printf("result-actual= x,y,z = %.13g, %.13g, %.13g \n", - (outCube2->x - actualCube->x), (outCube2->y - actualCube->y), - (outCube2->z - actualCube->z) ); + if (VERBOSE) { + printf(" -- expectCube = 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("expect-actual= x,y,z = %.13g, %.13g, %.13g \n", + (x - actualCube->x), (y - actualCube->y), (z - actualCube->z) ); + printf("expect-result= x,y,z = %.13g, %.13g, %.13g \n", + (x - outCube2->x), (y - outCube2->y), (z - outCube2->z) ); + printf("result-actual= x,y,z = %.13g, %.13g, %.13g \n", + (outCube2->x - actualCube->x), (outCube2->y - actualCube->y), + (outCube2->z - actualCube->z) ); + } // return 1; // } @@ -278,6 +297,8 @@ psFree(result2); result2 = psSphereGetOffset(actual, apparent, PS_SPHERICAL, PS_RADIAN); - printf("The apparent output sphere = r,d = %.13g, %.13g\n", result2->r, result2->d); - printf("The expected output sphere = r,d = %.13g, %.13g\n\n", result->r, result->d); + if (VERBOSE) { + printf("The apparent output sphere = r,d = %.13g, %.13g\n", result2->r, result2->d); + printf("The expected output sphere = r,d = %.13g, %.13g\n\n", result->r, result->d); + } psFree(result2); @@ -339,9 +360,11 @@ return 4; } - printf("\n PrecessionModel output = x,y,s = %.13g, %.13g, %.13g\n", - pmodel->x, pmodel->y, pmodel->s); - printf(" Expected output = x,y,s = %.13g, %.13g, %.13g\n", x, y, s); - printf(" A difference of: %.13g, %.13g, %.13g\n", - (pmodel->x - x), (pmodel->y - y), (pmodel->s - s) ); + if (VERBOSE) { + printf("\n PrecessionModel output = x,y,s = %.13g, %.13g, %.13g\n", + pmodel->x, pmodel->y, pmodel->s); + printf(" Expected output = x,y,s = %.13g, %.13g, %.13g\n", x, y, s); + printf(" A difference of: %.13g, %.13g, %.13g\n", + (pmodel->x - x), (pmodel->y - y), (pmodel->s - s) ); + } } psFree(pmodel); @@ -393,6 +416,8 @@ return 7; } else { - printf("\nPrecessionCorr output (IERSA) = x,y,s = %.13g, %.13g, %.13g\n", - pcorr->x, pcorr->y, pcorr->s); + if (VERBOSE) { + printf("\nPrecessionCorr output (IERSA) = x,y,s = %.13g, %.13g, %.13g\n", + pcorr->x, pcorr->y, pcorr->s); + } } psFree(pcorr); @@ -415,9 +440,11 @@ // psError(PS_ERR_BAD_PARAMETER_VALUE, false, // " psEOC_PrecessionCorr return incorrect values.\n"); - printf("PrecessionCorr output (IERSB) = x,y,s = %.13g, %.13g, %.13g\n", - pcorr->x, pcorr->y, pcorr->s); - printf("Expected output = x,y,s = %.13g, %.13g, %.13g\n", xx, yy, ss); - printf(" A difference of: %.13g, %.13g, %.13g\n\n", - (pcorr->x - xx), (pcorr->y - yy), (pcorr->s - ss) ); + if (VERBOSE) { + printf("PrecessionCorr output (IERSB) = x,y,s = %.13g, %.13g, %.13g\n", + pcorr->x, pcorr->y, pcorr->s); + printf("Expected output = x,y,s = %.13g, %.13g, %.13g\n", xx, yy, ss); + printf(" A difference of: %.13g, %.13g, %.13g\n\n", + (pcorr->x - xx), (pcorr->y - yy), (pcorr->s - ss) ); + } // return 10; // } @@ -449,9 +476,11 @@ // printf(" Output from CEOtoGCRS only = %.13g,%.13g,%.13g,%.13g\n", // pni->q0, pni->q1, pni->q2, pni->q3); - printf(" Expected sphere rotation = %.13g, %.13g, %.13g\n", q0,q1,q2); - printf(" Result sphere rotation = %.13g, %.13g, %.13g\n", - precessNutInv->q0, precessNutInv->q1, precessNutInv->q2); - printf(" Difference = %.13g, %.13g, %.13g\n\n", - precessNutInv->q0-q0, precessNutInv->q1-q1, precessNutInv->q2-q2); + if (VERBOSE) { + printf(" Expected sphere rotation = %.13g, %.13g, %.13g\n", q0,q1,q2); + printf(" Result sphere rotation = %.13g, %.13g, %.13g\n", + precessNutInv->q0, precessNutInv->q1, precessNutInv->q2); + printf(" Difference = %.13g, %.13g, %.13g\n\n", + precessNutInv->q0-q0, precessNutInv->q1-q1, precessNutInv->q2-q2); + } psCube *objC = psCubeAlloc(); // objC->x = -3.5963388069046304; @@ -479,5 +508,7 @@ y = 0.5555012823608123; z = 0.7496183628158023; - printf("\n<x, objC->y, objC->z); - printf(" Difference = %.13g, %.13g, %.13g\n\n", objC->x-x, objC->y-y, objC->z-z); - double xx = greatCircle(result, expect); - printf("GREAT CIRCLE DIFFERENCE = %.13g \n", xx); + if (VERBOSE) { + printf("<x, objC->y, objC->z); + printf(" Difference = %.13g, %.13g, %.13g\n\n", + objC->x-x, objC->y-y, objC->z-z); + printf("GREAT CIRCLE DIFFERENCE = %.13g \n", xx); + } + psFree(precess); psFree(precessNut); @@ -556,6 +590,8 @@ // psError(PS_ERR_BAD_PARAMETER_VALUE, false, // "psEOC_GetPolarMotion returned incorrect values.\n"); - printf(" <>PolarMotion output (IERSA) = x,y,s = %.13g, %.13g, %.13g\n", - polarMotion->x, polarMotion->y, polarMotion->s); + if (VERBOSE) { + printf(" <>PolarMotion output (IERSA) = x,y,s = %.13g, %.13g, %.13g\n", + polarMotion->x, polarMotion->y, polarMotion->s); + } // printf(" <>PolarMotion expected (IERSA) = x,y,s = %.13g, %.13g, %.13g\n", // x, y, s); @@ -578,6 +614,8 @@ // psError(PS_ERR_BAD_PARAMETER_VALUE, false, // "psEOC_GetPolarMotion returned incorrect values.\n"); - printf(" <>PolarMotion output (IERSB) = x,y,s = %.13g, %.13g, %.13g\n", - polarMotion->x, polarMotion->y, polarMotion->s); + if (VERBOSE) { + printf(" <>PolarMotion output (IERSB) = x,y,s = %.13g, %.13g, %.13g\n", + polarMotion->x, polarMotion->y, polarMotion->s); + } // printf(" <>PolarMotion expected (IERSB) = x,y,s = %.13g, %.13g, %.13g\n", // x, y, s); @@ -603,9 +641,11 @@ // psError(PS_ERR_BAD_PARAMETER_VALUE, false, // " psEOC_GetPolarMotion returned incorrect values.\n"); - printf("\n PolarMotion + NutationCorr out = x,y,s = %.13g, %.13g, %.13g\n", - polarMotion->x, polarMotion->y, polarMotion->s); - printf(" Expected output = x,y,s = %.13g, %.13g, %.13g\n", x, y, s); - printf(" Difference = x,y,s = %.13g, %.13g, %.13g\n", - polarMotion->x - x, polarMotion->y - y, polarMotion->s - s); + if (VERBOSE) { + printf("\n PolarMotion + NutationCorr out = x,y,s = %.13g, %.13g, %.13g\n", + polarMotion->x, polarMotion->y, polarMotion->s); + printf(" Expected output = x,y,s = %.13g, %.13g, %.13g\n", x, y, s); + printf(" Difference = x,y,s = %.13g, %.13g, %.13g\n", + polarMotion->x - x, polarMotion->y - y, polarMotion->s - s); + } // } @@ -645,6 +685,8 @@ return 2; } else { - printf("\nPolarTideCorr output = x,y,s = %.13g, %.13g, %.13g\n", - eop->x, eop->y, eop->s); + if (VERBOSE) { + printf("\nPolarTideCorr output = x,y,s = %.13g, %.13g, %.13g\n", + eop->x, eop->y, eop->s); + } } @@ -689,6 +731,8 @@ return 3; } else { - printf("Nutation Correction output = x,y,s = %.13g, %.13g, %.13g\n\n", - nute->x, nute->y, nute->s); + if (VERBOSE) { + printf("Nutation Correction output = x,y,s = %.13g, %.13g, %.13g\n\n", + nute->x, nute->y, nute->s); + } } psFree(nute); @@ -734,12 +778,16 @@ psError(PS_ERR_BAD_PARAMETER_VALUE, false, "psSphereRot_TEOtoCEO failed to return matching values for different time types.\n"); - printf("\n Output Rotation1 = q0,q1,q2,q3 = %.13g, %.13g, %.13g, %.13g\n", + if (VERBOSE) { + printf("\n Output Rotation1 = q0,q1,q2,q3 = %.13g, %.13g, %.13g, %.13g\n", + teoceo->q0, teoceo->q1, teoceo->q2, teoceo->q3 ); + printf("\n Output Rotation2 = q0,q1,q2,q3 = %.13g, %.13g, %.13g, %.13g\n", + rot->q0, rot->q1, rot->q2, rot->q3 ); + } + return 2; + } + if (VERBOSE) { + printf("\n Output Rotation = q0,q1,q2,q3 = %.13g, %.13g, %.13g, %.13g\n", teoceo->q0, teoceo->q1, teoceo->q2, teoceo->q3 ); - printf("\n Output Rotation2 = q0,q1,q2,q3 = %.13g, %.13g, %.13g, %.13g\n", - rot->q0, rot->q1, rot->q2, rot->q3 ); - return 2; - } - printf("\n Output Rotation = q0,q1,q2,q3 = %.13g, %.13g, %.13g, %.13g\n", - teoceo->q0, teoceo->q1, teoceo->q2, teoceo->q3 ); + } objSetup(); @@ -757,9 +805,11 @@ 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)); + if (VERBOSE) { + 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; } @@ -808,9 +858,11 @@ } - printf("\n Output sphere rotation = %.13g, %.13g, %.13g, %.13g\n", - rot->q0, rot->q1, rot->q2, rot->q3); - printf(" Expected sphere rotation = %.13g, %.13g, %.13g, %.13g\n", q0,q1,q2,q3); - printf(" difference: %.13g, %.13g, %.13g \n", - (rot->q0-q0), (rot->q1-q1), (rot->q2-q2) ); + if (VERBOSE) { + printf("\n Output sphere rotation = %.13g, %.13g, %.13g, %.13g\n", + rot->q0, rot->q1, rot->q2, rot->q3); + printf(" Expected sphere rotation = %.13g, %.13g, %.13g, %.13g\n", q0,q1,q2,q3); + printf(" difference: %.13g, %.13g, %.13g \n", + (rot->q0-q0), (rot->q1-q1), (rot->q2-q2) ); + } if (fabs(rot->q0-q0) > FLT_EPSILON || fabs(rot->q1-q1) > FLT_EPSILON || fabs(rot->q2-q2) > FLT_EPSILON || fabs(rot->q3+q3) > FLT_EPSILON) { @@ -832,8 +884,10 @@ y = 0.5555012823608123; z = 0.7496183628158023; - printf("\n Output 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)); + if (VERBOSE) { + printf("\n Output 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)); + } psCube *expect = psCubeAlloc(); expect->x = x; @@ -842,5 +896,7 @@ psSphere *expected = psCubeToSphere(expect); double d = greatCircle(result, expected); - printf(" The great circle angular distance between expected & result = %.13g\n", d); + if (VERBOSE) { + printf(" The great circle angular distance between expected & result = %.13g\n", d); + } psFree(expect); @@ -895,14 +951,19 @@ psError(PS_ERR_BAD_PARAMETER_VALUE, false, "psSphereRot_ITRStoTEO failed to return expected values.\n"); + if (VERBOSE) { + printf("\n Output sphere rotation = %.13g, %.13g, %.13g, %.13g\n", + rot->q0, rot->q1, rot->q2, rot->q3); + printf(" Expected sphere rotation = %.13g, %.13g, %.13g, %.13g\n", q0,q1,q2,q3); + printf(" a difference: %.13g, %.13g, %.13g, %.13g \n", (rot->q0-q0), (rot->q1-q1), + (rot->q2-q2), (rot->q3-q3) ); + } + printf("Error #3\n"); + // return 3; + } + if (VERBOSE) { printf("\n Output sphere rotation = %.13g, %.13g, %.13g, %.13g\n", rot->q0, rot->q1, rot->q2, rot->q3); printf(" Expected sphere rotation = %.13g, %.13g, %.13g, %.13g\n", q0,q1,q2,q3); - printf(" a difference: %.13g, %.13g, %.13g, %.13g \n", (rot->q0-q0), (rot->q1-q1), - (rot->q2-q2), (rot->q3-q3) ); - // return 3; - } - printf("\n Output sphere rotation = %.13g, %.13g, %.13g, %.13g\n", - rot->q0, rot->q1, rot->q2, rot->q3); - printf(" Expected sphere rotation = %.13g, %.13g, %.13g, %.13g\n", q0,q1,q2,q3); + } psCube *temp = psCubeAlloc(); @@ -920,10 +981,13 @@ test = psSphereRotApply(NULL, newRot, obj); temp = psSphereToCube(test); - printf("\n Cube -test- has x,y,z = %.13g, %.13g, %.13g \n", temp->x, temp->y, temp->z); double x, y, z; x = 0.01698577185310146; y = -0.6616538927902393; z = 0.7496169753347885; - printf("\n Cube -expected- has x,y,z = %.13g, %.13g, %.13g \n", x, y, z ); + if (VERBOSE) { + printf("\n Cube -test- has x,y,z = %.13g, %.13g, %.13g \n", + temp->x, temp->y, temp->z); + printf("\n Cube -expected- has x,y,z = %.13g, %.13g, %.13g \n", x, y, z ); + } temp->x = x; temp->y = y; @@ -931,5 +995,7 @@ psSphere *sphere = psCubeToSphere(temp); double d = greatCircle(sphere, test); - printf("Great circle difference of: %.13g \n", d); + if (VERBOSE) { + printf("Great circle difference of: %.13g \n", d); + } psFree(sphere); Index: /trunk/psLib/test/astro/verified/tst_psEarthOrientation.stdout =================================================================== --- /trunk/psLib/test/astro/verified/tst_psEarthOrientation.stdout (revision 7117) +++ /trunk/psLib/test/astro/verified/tst_psEarthOrientation.stdout (revision 7118) @@ -1,76 +1,5 @@ - - -- resultCube = x,y,z = -0.359633878262, 0.5555192580977, 0.7497078280829 -- - -- expectedCube = x,y,z = -0.3596338806905, 0.5555192509817, 0.7497078321908 -- -Cube Difference = x,y,z = -2.428477863248e-09, -7.116060296575e-09, 4.107927886743e-09 - - The great circle angular distance between expected & result = 0 -sunCube = x,y,z = 0.9820293796219, 0.1731564131522, 0.0750676637684 -sunSphere = r, d = 0.1745310859793, 0.0751383461293 Theta = 1.772849258423 r0 = -7.3028651010049865723e+11 Theta = -1.913704437867e-08 phi = 7.882970755914e-09 - -- actualCube = x,y,z = -0.3596195125758, 0.5555613903456, 0.7496834983725 -- - -- resultCube2= x,y,z = -0.3596194987327, 0.5555613922666, 0.7496835035894 -- - -- expectCube = x,y,z = -0.3596194976029, 0.5555613950298, 0.7496835020836 -- - -expect-actual= x,y,z = 1.497289375818e-08, 4.684221899254e-09, 3.711128404582e-09 -expect-result= x,y,z = 1.129748239403e-09, 2.763238637904e-09, -1.505794267054e-09 -result-actual= x,y,z = 1.384314551878e-08, 1.92098326135e-09, 5.216922671636e-09 -The apparent output sphere = r,d = -1.91370443936e-08, 7.88297072063e-09 -The expected output sphere = r,d = -2.283899158684e-08, 5.607657671547e-09 - - - PrecessionModel output = x,y,s = 0.0002857180163421, 2.394794321783e-05, -1.397070873644e-08 - Expected output = x,y,s = 0.0002857175590089, 2.396873937773e-05, -1.39700664579e-08 - A difference of: 4.57333152468e-10, -2.079615990713e-08, -6.422785401238e-13 - -PrecessionCorr output (IERSA) = x,y,s = 1.947922616393e-10, 1.721410388291e-10, 0 -PrecessionCorr output (IERSB) = x,y,s = 3.054127043514e-10, -1.396712376063e-10, 0 -Expected output = x,y,s = 3.052243007204e-10, -1.394413392358e-10, 0 - A difference of: 1.884036310428e-13, -2.298983704727e-13, 0 - - Expected sphere rotation = -1.198452240676e-05, 0.0001428589335861, 1.219119351891e-10 - Result sphere rotation = -1.197390289426e-05, 0.0001428591622618, 6.98535429644e-09 - Difference = 1.061951249492e-08, 2.286756638811e-10, 6.863442361251e-09 - - -<PolarMotion output (IERSA) = x,y,s = -6.454357026873e-07, 2.112611340557e-06, 0 - <>PolarMotion output (IERSB) = x,y,s = -6.453778846156e-07, 2.112824593826e-06, 0 - -MJD check = 52730.06284722 - - PolarMotion + NutationCorr out = x,y,s = -6.453259754967e-07, 2.113939222748e-06, -7.396173345101e-12 - Expected output = x,y,s = -6.43607313124e-07, 2.113514369736e-06, -7.396175813246e-12 - Difference = x,y,s = -1.718662372636e-09, 4.248530125326e-10, 2.468145872205e-18 - -PolarTideCorr output = x,y,s = 5.576394041032e-11, 1.141315275255e-09, -2.661051633224e-05 - -MJD check = 52730.06284722 -Nutation Correction output = x,y,s = -3.854821468189e-12, -2.668635296536e-11, -7.396173345101e-12 - - - Output Rotation = q0,q1,q2,q3 = 0, 0, -0.9625401013456, 0.2711393614023 - - Output sphere rotation = -1.198437080989e-05, 0.0001428587808888, 6.985033157173e-09, 0.9999999897239 - Expected sphere rotation = -1.198452240676e-05, 0.0001428589335861, 1.219119351891e-10, -0.9999999897238 - difference: 1.515968626689e-10, -1.526973275742e-10, 6.863121221984e-09 - - Output cube = x,y,z = -0.3598480648445, 0.555501287526, 0.7496183627584 -Expected cube = x,y,z = -0.3598480726985, 0.5555012823608, 0.7496183628158 - A difference of: -7.854020045439e-09, -5.165185501532e-09, 5.738021169321e-11 - - The great circle angular distance between expected & result = 1.490116119385e-08 - - Output sphere rotation = -1.056757184868e-06, 3.218036565618e-07, -3.410955728293e-13, 0.9999999999994 - Expected sphere rotation = -1.056757184866e-06, 3.218036562932e-07, -3.35801958072e-12, -0.9999999999994 - - Cube -test- has x,y,z = 0.01698577184911, -0.6616538927903, 0.7496169753348 - - Cube -expected- has x,y,z = 0.0169857718531, -0.6616538927902, 0.7496169753348 -Great circle difference of: 0 +Error #3 Index: /trunk/psLib/test/math/tst_psPolyFit4D.c =================================================================== --- /trunk/psLib/test/math/tst_psPolyFit4D.c (revision 7117) +++ /trunk/psLib/test/math/tst_psPolyFit4D.c (revision 7118) @@ -33,5 +33,5 @@ #define ERROR_TOLERANCE 0.10 #define YERR 10.0 -#define VERBOSE 0 +#define VERBOSE 1 #define NUM_ITERATIONS 5 #define CLIP_SIGMA 4.0