Changeset 5450 for trunk/psLib/test/astro/tst_psSphereOps.c

Timestamp:

Oct 27, 2005, 4:25:22 PM (21 years ago)

Author:

drobbin

Message:

Implemented and tested psAberration, but results have not been verified.

File:

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

trunk/psLib/test/astro/tst_psSphereOps.c

@@ -6 +6 @@
 *  @author d-Rob, MHPCC
 *
-*  @version $Revision: 1.10 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-10-26 01:20:15 $
+*  @version $Revision: 1.11 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-10-28 02:25:22 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -71 +71 @@
     double q3 = cos(a3)*cos(DELTA_P/2);
 
-    if (FLT_EPSILON < fabs(q0 - myST->q0)) {
+    if (DBL_EPSILON < fabs(q0 - myST->q0)) {
         psError(PS_ERR_UNKNOWN,true,"myST->q0 is %lf, should be %lf\n", myST->q0, q0);
         return 2;
     }
-    if (FLT_EPSILON < fabs(q1 - myST->q1)) {
+    if (DBL_EPSILON < fabs(q1 - myST->q1)) {
         psError(PS_ERR_UNKNOWN,true,"myST->q1 is %f, should be %f\n", myST->q1, q1);
         return 3;
     }
-    if (FLT_EPSILON < fabs(q2 - myST->q2)) {
+    if (DBL_EPSILON < fabs(q2 - myST->q2)) {
         psError(PS_ERR_UNKNOWN,true,"myST->q2 is %f, should be %f\n", myST->q2, q2);
         return 4;
     }
-    if (FLT_EPSILON < fabs(q3 - myST->q3)) {
+    if (DBL_EPSILON < fabs(q3 - myST->q3)) {
         psError(PS_ERR_UNKNOWN,true,"myST->q0 is %f, should be %f\n", myST->q3, q3);
         return 5;
@@ -239 +239 @@
 
         // check ecliptic transforms for correctness
-        if (abs(RAD_TO_DEG(ecliptic->r) - lambda[x]) > TOLERANCE ||
-                abs(RAD_TO_DEG(ecliptic->d) - beta[x]) > TOLERANCE) {
+        if (fabs(RAD_TO_DEG(ecliptic->r) - lambda[x]) > TOLERANCE ||
+                fabs(RAD_TO_DEG(ecliptic->d) - beta[x]) > TOLERANCE) {
             psError(PS_ERR_UNKNOWN, false,
                     "Ecliptic tranformation incorrect.  Result is (%g,%g), expected (%g,%g)",
@@ -247 +247 @@
             return 1;
         }
-        if (abs(RAD_TO_DEG(icrsFromEcliptic->r) - alpha[x]) > TOLERANCE ||
-                abs(RAD_TO_DEG(icrsFromEcliptic->d) - delta[x]) > TOLERANCE) {
+        if (fabs(RAD_TO_DEG(icrsFromEcliptic->r) - alpha[x]) > TOLERANCE ||
+                fabs(RAD_TO_DEG(icrsFromEcliptic->d) - delta[x]) > TOLERANCE) {
             psError(PS_ERR_UNKNOWN, false,
                     "ICRS for Ecliptic tranformation incorrect.  Result is (%g,%g), expected (%g,%g)",
@@ -269 +269 @@
 
         // check ecliptic transforms for correctness
-        if (abs(RAD_TO_DEG(galactic->r) - l[x]) > TOLERANCE ||
-                abs(RAD_TO_DEG(galactic->d) - b[x]) > TOLERANCE) {
+        if (fabs(RAD_TO_DEG(galactic->r) - l[x]) > TOLERANCE ||
+                fabs(RAD_TO_DEG(galactic->d) - b[x]) > TOLERANCE) {
             psError(PS_ERR_UNKNOWN, false,
                     "Galactic tranformation incorrect.  Result is (%g,%g), expected (%g,%g)",
@@ -280 +280 @@
             return 3;
         }
-        if (abs(RAD_TO_DEG(icrsFromGalactic->r) - alpha[x]) > TOLERANCE ||
-                abs(RAD_TO_DEG(icrsFromGalactic->d) - delta[x]) > TOLERANCE) {
+        if (fabs(RAD_TO_DEG(icrsFromGalactic->r) - alpha[x]) > TOLERANCE ||
+                fabs(RAD_TO_DEG(icrsFromGalactic->d) - delta[x]) > TOLERANCE) {
             psError(PS_ERR_UNKNOWN, false,
                     "ICRS for Galactic tranformation incorrect.  Result is (%g,%g), expected (%g,%g)",
@@ -498 +498 @@
     output->d *= -1.0;
     empty = psSphereGetOffset(origin, output, PS_LINEAR, PS_RADIAN);
-    if ( fabs(offset->r - empty->r) > 0.0001 || fabs(offset->d - empty->r) > 0.0001 ) {
+    if ( fabs(offset->r - empty->r) > 0.0001 || abs(offset->d - empty->r) > 0.0001 ) {
         psError(PS_ERR_BAD_PARAMETER_VALUE, false,
                 "psSphereGetOffset failed to return correct linear offset values.\n");