---

polytest.cpp File Reference

Test for polynomial interpolation wavelets. More...

#include "waveletDebug.h"
#include "polyHaar.h"

Go to the source code of this file.

Functions
bool	compare (const double *v1, const double *v2, const size_t N)
	Compare two integer arrays of size N. More...
double*	copy (double *vec, const size_t N)
	Make a copy of an integer array. More...
	main ()
	\function. More...
Variables
double	data []

---

Detailed Description

Test for polynomial interpolation wavelets.

Definition in file polytest.cpp.

---

Function Documentation

bool compare ( const double * v1, const double * v2, const size_t N )

Compare two integer arrays of size N. Return true if they are equal, false otherwise. Definition at line 18 of file polytest.cpp. { bool rslt = true; for (size_t i = 0; i < N; i++) { if (v1[i] != v2[i]) { rslt = false; break; } } return rslt; } // compare

double * copy ( double * vec, const size_t N )

Make a copy of an integer array. Note that the function allocates memory, but does not deallocate it. Definition at line 40 of file polytest.cpp. { double *newVec = new double[ N ]; for (size_t i = 0; i < N; i++) { newVec[i] = vec[i]; } return newVec; } // copy

main ( )

\function. The entry point for code to test the wavelet packet transform. The code in main provides a simple example of how to call the wavelet packet transform code. This test case can be used with various wavelet functions (e.g., Haar, linear interpolation, Daubechies D4, polynomial interpolation). In the context of filters, the polynomial interpolation wavelets (see poly.h and polyHaar.h) are very bad choices. The polynomial interpolation wavelets do not divide the spectrum into high and low frequency components. Instead, the power spectrum shows one peak. I'm not sure if this means that the polynomial interpolation wavelets are bad for all applications. For example, another way to look at wavelets is as approximation functions. For example, wavelet compression relies approximation functions. In this case polynomial wavelets may be a good choice for some data sets. Definition at line 53 of file polytest.cpp. { waveletDebug d; polyHaar<double *> w; size_t N = sizeof( data) / sizeof(double); double *copyVec = copy( data, N ); double *ptr = data; w.forwardTrans( ptr, N ); d.pr_ordered( data, N ); w.inverseTrans( ptr, N ); if (compare(data, copyVec, N)) printf("inverse transform succeeded\n"); else printf("inverse transform failed\n"); return 0; }

---

Variable Documentation

double data[]

Initial value: { 32.0, 10.0, 20.0, 38.0, 37.0, 28.0, 38.0, 34.0, 18.0, 24.0, 18.0, 9.0, 23.0, 24.0, 28.0, 34.0 } Definition at line 9 of file polytest.cpp.

---