#include <stdio.h>
#define L 4
static const float h[L] = {  0.48296291314453416,  0.83651630373780794,
                             0.22414386804201339, -0.12940952255126037};
static const float g[L] = { -0.12940952255126037, -0.22414386804201339,
                             0.83651630373780794, -0.48296291314453416};
#include "ipcqfilt.c"
#include "ipdwt.c"
int main(void) {
  float scaling[128]={0}, wavelet[128]={0};  int i;
  scaling[1]=1.0;  wavelet[5]=1.0;

  ipdwt( scaling, 128, 5 );
  puts("# 128 point Daubechies 4 scaling function, level 5:");
  for(i=0; i<128; i++) printf("%10.6f\n", scaling[i]);

  ipdwt( wavelet, 128, 5 );
  puts("# 128 point Daubechies 4 wavelet, level 5:");
  for(i=0; i<128; i++) printf("%10.6f\n", wavelet[i]);

  return 0;
}