#include <stdio.h>
#include "wslift.c"

int main(void) {
  float u0[100]={0}, u1[100]={0}, u2[100]={0};
  int i, index[3]={32,36,41}, level[3]={3,4,5};
  u0[index[0]]=1.0;
  u1[index[1]]=1.0;
  u2[index[2]]=1.0;

  wsl97idwt( u0, 100, 1, 3 );
  printf("# 100 point 9,7 wavelet, index %d, level %d:\n", index[0], level[0]);
  for(i=0; i<100; i++) printf("%10.6f\n", u0[i]);

  wsl97idwt( u1, 100, 1, 4 );
  printf("# 100 point 9,7 wavelet, index %d, level %d:\n", index[1], level[1]);
  for(i=0; i<100; i++) printf("%10.6f\n", u1[i]);

  wsl97idwt( u2, 100, 1, 5 );
  printf("# 100 point 9,7 wavelet, index %d, level %d:\n", index[2], level[2]);
  for(i=0; i<100; i++) printf("%10.6f\n", u2[i]);


  return 0;
}