Re: [Bug-gsl] Discrete Hankel Transform - unexpected behaviour

[Brian Gough]

At Wed, 20 Aug 2008 16:15:55 +1200,
Luca Fascione wrote:
> I'm using your dht code to transform some data I need to convolve.
> In trying to understand how many samples of my function I needed I 
> conducted the following experiment:
>  - I generated a dht class with fixed maximum value and various sample 
> counts (128, 256, 512, 1024, 2048)
>  - I transformed and plotted the various datasets
> 
> I was very surprised when I realized the transformed dataset doesn't 
> seem to converge to a defined target value,
> but keeps lowering in magnitude instead.

I'm not familiar with the Hankel transform myself, but have you
checked that the values agree with the definition and that the inverse
returns the original function (within a constant factor)?  I imagine
that the normalisation factor varies with the number of sample points.

-- 
Brian Gough

GNU Scientific Library -
http://www.gnu.org/software/gsl/