fwt.h

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

/* 棧嶶僂僃乕僽儗僢僩曄姺 僿僢僟乕僼傽僀儖

  Copyright 1998,1999 H. Nakano 
*/

#define L  64    /* 1師尦怣崋偺挿偝 */
#define L2 32    /* L/2 */ 
#define N  64    /* 2師尦怣崋(惓曽攝楍)偺堦曈偺挿偝 */
#define N2 32    /* N/2 */
#define NN 4096  /* N*N */

#define K  4    /* 僂僃乕僽儗僢僩悢楍 p_k 偺挿偝 */

#define PI  3.14159265358979

#ifndef max
 #define	max(A, B)	((A) > (B) ? (A) : (B))
#endif
/*****************************************************************************/

void fwt2d (
	       double s0[N][N],
		unsigned int s_len,
		double p[],
	       double q[],
		unsigned int sup,
	       double s1[N2][N2],
	       double w1v[N2][N2],
		double w1h[N2][N2],
		double w1d[N2][N2]);

void ifwt2d (
	       double s1[N2][N2],
	       double w1v[N2][N2],
		double w1h[N2][N2],
		double w1d[N2][N2],
		unsigned int s_len,
		double p[],
	       double q[],
		unsigned int sup,
   	       double s0[N][N]
		   );


void fwt1d (
	       double s0[],
		unsigned int s_len,
		double p[],
	       double q[],
		unsigned int sup,
	       double s1[],
	       double w1[]);

void ifwt1d (
	       double s1[],
	       double w1[],	  
		unsigned int s_len,
		double p[],
	       double q[],
		unsigned int sup,
		double s0[]
);

void main(void);