36   float *fimg = 0, *temp, thold, mul[2], avg, diff;  in wavelet_denoise()  local
48     fimg = (float *)malloc((size * 3 + iheight + iwidth) * sizeof *fimg);  in wavelet_denoise()
49   merror(fimg, "wavelet_denoise()");  in wavelet_denoise()
50   temp = fimg + size * 3;  in wavelet_denoise()
56       fimg[i] = 256 * sqrt((double)(image[i][c] << scale));  in wavelet_denoise()
62         hat_transform(temp, fimg + hpass + row * iwidth, 1, iwidth, 1 << lev);  in wavelet_denoise()
64           fimg[lpass + row * iwidth + col] = temp[col] * 0.25;  in wavelet_denoise()
68         hat_transform(temp, fimg + lpass + col, iwidth, iheight, 1 << lev);  in wavelet_denoise()
70           fimg[lpass + row * iwidth + col] = temp[row] * 0.25;  in wavelet_denoise()
75         fimg[hpass + i] -= fimg[lpass + i];  in wavelet_denoise()
76         if (fimg[hpass + i] < -thold)  in wavelet_denoise()
77           fimg[hpass + i] += thold;  in wavelet_denoise()
78         else if (fimg[hpass + i] > thold)  in wavelet_denoise()
79           fimg[hpass + i] -= thold;  in wavelet_denoise()
81           fimg[hpass + i] = 0;  in wavelet_denoise()
83           fimg[i] += fimg[hpass + i];  in wavelet_denoise()
88       image[i][c] = CLIP(SQR(fimg[i] + fimg[lpass + i]) / 0x10000);  in wavelet_denoise()
98       window[i] = (ushort *)fimg + width * i;  in wavelet_denoise()
127   free(fimg);  in wavelet_denoise()
132   float *fimg = 0, *temp, thold, mul[2], avg, diff;  in wavelet_denoise()  local
144     fimg = (float *)malloc((size * 3 + iheight + iwidth) * sizeof *fimg);  in wavelet_denoise()
145   merror(fimg, "wavelet_denoise()");  in wavelet_denoise()
146   temp = fimg + size * 3;  in wavelet_denoise()
153     temp = (float *)malloc((iheight + iwidth) * sizeof *fimg);  in wavelet_denoise()
158         fimg[i] = 256 * sqrt((double)(image[i][c] << scale));  in wavelet_denoise()
165           hat_transform(temp, fimg + hpass + row * iwidth, 1, iwidth, 1 << lev);  in wavelet_denoise()
167             fimg[lpass + row * iwidth + col] = temp[col] * 0.25;  in wavelet_denoise()
172           hat_transform(temp, fimg + lpass + col, iwidth, iheight, 1 << lev);  in wavelet_denoise()
174             fimg[lpass + row * iwidth + col] = temp[row] * 0.25;  in wavelet_denoise()
180           fimg[hpass + i] -= fimg[lpass + i];  in wavelet_denoise()
181           if (fimg[hpass + i] < -thold)  in wavelet_denoise()
182             fimg[hpass + i] += thold;  in wavelet_denoise()
183           else if (fimg[hpass + i] > thold)  in wavelet_denoise()
184             fimg[hpass + i] -= thold;  in wavelet_denoise()
186             fimg[hpass + i] = 0;  in wavelet_denoise()
188             fimg[i] += fimg[hpass + i];  in wavelet_denoise()
194         image[i][c] = CLIP(SQR(fimg[i] + fimg[lpass + i]) / 0x10000);  in wavelet_denoise()
211       window[i] = (ushort *)fimg + width * i;  in wavelet_denoise()
240   free(fimg);  in wavelet_denoise()