revised version of wavelet transform

Brian Gough bjg@network-theory.co.uk
Sun Apr 25 20:23:00 GMT 2004


I've had a look at this now.  It is a good improvement.

A couple of points I noticed,

- in the initialisation of h1, h2, g1, g2 in daubechies and bspline
there are non-constant static arrays, which are not thread-safe.  To
avoid this, I'd recommend supplying all the coefficients as
precomputed tables (even though they only differ in sign).

- there could be some value in putting the allowed wavelet lengths in
an array in the wavelet_type struct, to allow a calling routine to
inspect them. The array could be terminated with 0.

-- 
Brian



More information about the Gsl-discuss mailing list