33 #include <casu_utils.h>
34 #include <casu_filt.h>
38 static void initdata(
float *data);
48 cpl_test_init(PACKAGE_BUGREPORT,CPL_MSG_WARNING);
52 data = cpl_malloc(DATSZ*DATSZ*
sizeof(
float));
53 bpm = cpl_calloc(DATSZ*DATSZ,
sizeof(
float));
55 testpt = nbpm + DATSZ;
60 casu_bfilt(data,bpm,DATSZ,DATSZ,5,MEANCALC,1);
62 cpl_test_rel(data[testpt],111.0,1.0e-6);
67 casu_bfilt(data,bpm,DATSZ,DATSZ,5,MEDIANCALC,1);
69 cpl_test_rel(data[testpt],111.0,1.0e-6);
76 casu_bfilt(data,bpm,DATSZ,DATSZ,5,MEANCALC,1);
78 cpl_test_rel(data[testpt],111.25,1.0e-6);
83 memset(bpm,0,DATSZ*DATSZ);
85 casu_bfilt(data,bpm,DATSZ,DATSZ,5,MEDIANCALC,1);
86 cpl_test_rel(data[testpt],111.5,1.0e-6);
93 return(cpl_test_end(0));
96 static void initdata(
float *data) {
103 for (j = 0; j <= DATSZ-1; j++) {
104 for (i = 0; i <= DATSZ-1; i++) {
106 if (i >= mid-3 && i <= mid+3 && j >= mid-3 && j <= mid+3) {