fwtest.c

2D weblet programming2D weblet programming 2D weblet programming 2D weblet programming 2D weblet

/******************************************************
 僒儞僾儖僾儘僌儔儉 丗 侾師尦崅懍僂僃乕僽儗僢僩曄姺
                           儊僀儞儖乕僠儞
 Copyright 1998,1999 H. Nakano
********************************************************/

#include <stdio.h>
#include <math.h>
#include "fwt.h"


void main(void){

   double s0[L];                    /* 擖椡怣崋 S0 */
   double s0r[L];                   /* 嵞峔惉怣崋 */
   double s1[L2],w1[L2];            /* 暘夝怣崋 s1, w1 */
   double p[K]={0.482962913145, 0.836516303738, 0.224143868042, -0.129409522551}; 
            /* 僪儀僔僀偺悢楍 p_k (N=2) */
   double q[K];  /* 僪儀僔僀偺悢楍 q_k (N=2) */
 
   int i;
   int sup=K,s_len=L;

	for(i=0;i<L;i++){  /* s0 偺弶婜抣愝掕 */
		s0[i]=0.0;
	}
	for(i=0;i<16;i++){  
		s0[i]=(float)((i+1)*(i+1))/256.0;
	}
	for(i=16;i<32;i++){
		s0[i]=.2;
	}
	for(i=32;i<48;i++){
		s0[i]=(float)((47-i)*(47-i))/256.0-.5;
	}


   for(i=0;i<sup;i++){  /* p_k 偐傜 q_k 傪惗惉 */
      q[i]=pow(-1,i)*p[sup-i-1];
	}

	printf("擖椡怣崋 \n");
	for(i=0;i<L;i++){
		printf("%7.4f ",s0[i]);
	}
	printf("\n");
	
   fwt1d(s0,s_len,p,q,sup,s1,w1);  /* 1師尦崅懍僂僃乕僽儗僢僩曄姺 */
	
   ifwt1d(s1,w1,s_len/2,p,q,sup,s0r);  /* 1師尦崅懍僂僃乕僽儗僢僩媡曄姺 */

	printf("嵞峔惉怣崋 \n");
	for(i=0;i<L;i++){
		printf("%7.4f ",s0r[i]);
	}
	printf("\n");
}