Changeset 2199 for trunk/psLib/test/astronomy/tst_psCoord.c
- Timestamp:
- Oct 26, 2004, 11:57:42 AM (22 years ago)
- File:
-
- 1 edited
-
trunk/psLib/test/astronomy/tst_psCoord.c (modified) (7 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/test/astronomy/tst_psCoord.c
r2108 r2199 6 6 * @author GLG, MHPCC 7 7 * 8 * @version $Revision: 1. 3$ $Name: not supported by cvs2svn $9 * @date $Date: 2004-10- 14 01:53:44$8 * @version $Revision: 1.4 $ $Name: not supported by cvs2svn $ 9 * @date $Date: 2004-10-26 21:57:32 $ 10 10 * 11 11 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 33 33 {test3, 0000, "psPlaneDistortApply()", 0, false}, 34 34 {test4, 0000, "psPSphereTransformApply()", 0, false}, 35 {test5, 0000, "psSphereTransformICRS toEcliptic()", 0, false},36 {test6, 0000, "psSphereTransformEcliptic toICRS()", 0, false},37 {test7, 0000, "psSphereTransformICRS toGalatic()", 0, false},38 {test8, 0000, "psSphereTransformGalatic toICRS()", 0, false},35 {test5, 0000, "psSphereTransformICRSToEcliptic()", 0, false}, 36 {test6, 0000, "psSphereTransformEclipticToICRS()", 0, false}, 37 {test7, 0000, "psSphereTransformICRSToGalatic()", 0, false}, 38 {test8, 0000, "psSphereTransformGalaticToICRS()", 0, false}, 39 39 {test20, 0000, "psProject()", 0, false}, 40 40 {test21, 0000, "psDeProject()", 0, false}, … … 76 76 float y) 77 77 { 78 if (y!= 0.0) { 79 return(atan(y/x)); 80 } else { 81 return(0.5 * M_PI); 82 } 83 return (0.0); 78 return((float) atan2((double) x, (double) y)); 84 79 } 85 80 … … 251 246 double T = year / 100.0; 252 247 double phi = -23.452294 + 0.013013 * T + 0.000001639 * T * T - 0.000000503 * T * T * T; 253 psSphereTransform *myST = psSphereTransformICRS toEcliptic(*now);248 psSphereTransform *myST = psSphereTransformICRSToEcliptic(*now); 254 249 255 250 if (FLT_EPSILON < fabs(sin(phi) - myST->sinPhi)) { … … 290 285 double T = year / 100.0; 291 286 double phi = +23.452294 - 0.013013 * T - 0.000001639 * T * T + 0.000000503 * T * T * T; 292 psSphereTransform *myST = psSphereTransformEcliptic toICRS(*now);287 psSphereTransform *myST = psSphereTransformEclipticToICRS(*now); 293 288 294 289 if (FLT_EPSILON < fabs(sin(phi) - myST->sinPhi)) { … … 326 321 double Xo = 282.25; 327 322 double xo = 33.0; 328 psSphereTransform *myST = psSphereTransformICRS toGalatic();323 psSphereTransform *myST = psSphereTransformICRSToGalatic(); 329 324 330 325 if (FLT_EPSILON < fabs(sin(phi) - myST->sinPhi)) { … … 359 354 double Xo = 33.0; 360 355 double xo = 282.25; 361 psSphereTransform *myST = psSphereTransformGalatic toICRS();356 psSphereTransform *myST = psSphereTransformGalaticToICRS(); 362 357 363 358 if (FLT_EPSILON < fabs(sin(phi) - myST->sinPhi)) {
Note:
See TracChangeset
for help on using the changeset viewer.
