discrete wavelet transform

[Ivo Alxneit]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

hi

i have put together functions to calculate the discrete wavelet transform in 
1D (gsl_vector_double) and 2D (gsl_matrix_double). i have already implemented 
different wavelets (beylkin, coiflet, daubechies, haar).

i have verified the output against paul addisons illustrated wavelet transform 
book and matlabs wavekit. the numerical output is the same (minus 
rounding/numerics) but note that the position of the T submatrices occuring 
at each levels is different in the book compared with the output of wavekit. 
the position of Th and Tv (the off-diagonal blocks are exchanged). the output 
of my implemetation is the same as in the book.

is this something that could be included in gsl? if yes i would still need to 
write some tests (suggestions welcome) and documentation.

the attached archive contains the code and illustrations for 1D and 2D.
- -- 
Dr. Ivo Alxneit
Laboratory for Solar Technology   phone: +41 56 310 4092
CH-5232 Villigen                    fax: +41 56 310 2624
Paul Scherrer Institute          http://solar.web.psi.ch
Switzerland                        gnupg key: 0x515E30C7
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.2.1 (GNU/Linux)

iD8DBQFAS0kbAd7CE1FeMMcRAkuhAKCjXhtqWMzGaP5BkQGI+Ngv7CSn1wCeMUEm
+1sbYaxt146BSNkwqjzJdmI=
=08yf
-----END PGP SIGNATURE-----
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gsl_wavelet.tar.gz
Type: application/x-tgz
Size: 6310 bytes
Desc: not available
URL: <http://sourceware.org/pipermail/gsl-discuss/attachments/20040307/04bf777b/attachment.bin>