#include<stdio.h>
#include<math.h>
void fwt1d(
	int j, /*現在のレベル*/
	double s0[],	/*入力信号*/
	unsigned int s_len,	/*入力信号長*/
	double p[],		/*スケーリング数列 pk*/
	double q[],		/*ウェーブレット数列 qk*/
	unsigned int sup,	/*数列長*/
	double s1[],	/*分解信号 s1*/
	double w1[])	/*分解信号 w1*/
{
	unsigned int n,k;
	unsigned int index;
	printf("ok!\n");
	for(k = 0; k < s_len/2; k ++){
	printf("%d\n",k);
		s1[k] = 0.0;
		w1[k] = 0.0;
		for(n = 0;n < sup; n++){
			index = (n + 2 * k) %s_len ;
			s1[k] += p[n] * s0[index];
			w1[k] += q[n] * s0[index];
		}
		w1[k] = 1/sqrt(pow(2,j))*w1[k];
	}
	
}