00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 #ifdef HAVE_CONFIG_H
00032 #include <config.h>
00033 #endif
00034
00035 #include <cpl_test.h>
00036
00037 #include "irplib_wcs.h"
00038
00039
00040
00041
00042 static void irplib_wcs_all_test(void);
00043
00044
00045
00046
00047 int main (void)
00048 {
00049
00050 cpl_test_init(PACKAGE_BUGREPORT, CPL_MSG_WARNING);
00051
00052 irplib_wcs_all_test();
00053
00054 return cpl_test_end(0);
00055 }
00056
00057 static void irplib_wcs_all_test(void)
00058 {
00059 cpl_propertylist * prop_wcs;
00060 cpl_wcs * wcs = NULL;
00061 double xorig,yorig;
00062 double xnew,ynew;
00063 double ra, dec;
00064
00065
00066 prop_wcs = cpl_propertylist_new();
00067 cpl_test_nonnull(prop_wcs);
00068 cpl_propertylist_append_double(prop_wcs, "CRVAL1", 10.);
00069 cpl_propertylist_append_double(prop_wcs, "CRVAL2", 20.);
00070 cpl_propertylist_append_int(prop_wcs, "CRPIX1", 1);
00071 cpl_propertylist_append_int(prop_wcs, "CRPIX2", 2);
00072 cpl_propertylist_append_double(prop_wcs, "CD1_1", 10.);
00073 cpl_propertylist_append_double(prop_wcs, "CD1_2", 11.);
00074 cpl_propertylist_append_double(prop_wcs, "CD2_1", 13.);
00075 cpl_propertylist_append_double(prop_wcs, "CD2_2", 14.);
00076 wcs = cpl_wcs_new_from_propertylist(prop_wcs);
00077 if(cpl_error_get_code() == CPL_ERROR_NO_WCS)
00078 {
00079 cpl_msg_warning(__func__,"No WCS present. Tests disabled");
00080 cpl_test_error(CPL_ERROR_NO_WCS);
00081 cpl_test_null(wcs);
00082 cpl_propertylist_delete(prop_wcs);
00083 return;
00084 }
00085 cpl_test_nonnull(wcs);
00086
00087
00088 xorig = 1;
00089 yorig = 2;
00090 irplib_wcs_xytoradec(wcs, xorig, yorig, &ra, &dec);
00091 cpl_test_error(CPL_ERROR_NONE);
00092
00093
00094 irplib_wcs_radectoxy(wcs, ra, dec, &xnew, &ynew);
00095 cpl_test_leq((xnew - xorig), 2 * DBL_EPSILON);
00096 cpl_test_leq((ynew - yorig), 2 * DBL_EPSILON);
00097 cpl_test_error(CPL_ERROR_NONE);
00098
00099
00100 cpl_wcs_delete(wcs);
00101 cpl_propertylist_delete(prop_wcs);
00102
00103 }