22#include "../hdrl_cat_seeing.h"
27 cpl_test_init(PACKAGE_BUGREPORT,CPL_MSG_WARNING);
32 double ell[] = {0.009, 0.050, 0.025, 0.033, 0.010, 0.012, 0.014, 0.005, 0.011, 0.011};
34 double pk[] = {1007.065, 101.293, 204.195, 493.250, 559.111, 609.799, 642.603, 698.117, 740.227, 797.222};
36 double areal[][10] = { {120., 53., 73., 104., 97., 109., 104., 107., 110., 110.},
37 {100., 39., 60., 75., 80., 81., 86., 90., 94., 95.},
38 { 75., 21., 42., 63., 66., 64., 68., 69., 69., 70.},
39 { 63., 2., 21., 45., 45., 45., 45., 47., 56., 56.},
40 { 45., 0., 4., 25., 29., 29., 37., 37., 37., 37.},
41 { 25., 0., 0., 9., 9., 13., 14., 21., 21., 21.},
42 { 9., 0., 0., 0., 0., 0., 0., 1., 1. },
43 { 0., 0., 0., 0., 0., 0., 0., 0. } };
52 double **areals = cpl_malloc(8 *
sizeof(
double *));
53 for (cpl_size i = 0; i < 8; i++) {
54 areals[i] = cpl_malloc(10 *
sizeof(
double));
56 for (cpl_size j = 0; j < 10; j++) {
57 areals[i][j] = areal[i][j];
62 double *work = cpl_malloc(10 *
sizeof(
double));
65 cpl_test_rel(fwhm, 4.50384, 0.01);
69 for (cpl_size i = 0; i < 8; i++) {
75 return cpl_test_end(0);
cpl_error_code hdrl_seeing(ap_t *ap, cpl_size nrows, double *ellipt, double *pkht, double **areal, double *work, double *fwhm)
Work out the median seeing.