[Bug-gsl] Bug in dht

[Benno Rem]

Dear all,

I noticed that, as soon as I try to do a Hankel Transform (nu = 0) on a 
gaussian of the form exp( - pi * r^2 ), transform it back and multiply by the 
factor Jzero0(0,N)^2, that I don't get the exactly same function back.
Actually, the function that is returned is scaled by a factor 1.00391 with 
respect to the original function.

For just one transform it can be considered as just a numerical error, but as 
soon as I try to do some more complicated stuff I get strong abbreviations from 
the real result.

I'm not sure if this is a bug, or just a misunderstanding of the concept, but 
it would help me a lot if someone could give me an idea.

Regards,

Benno REM
address@hidden
------------------------------------------
Département de Physique de l'Ecole Normale Supérieure
Laboratoire Kastler Brossel
24 rue Lhomond 75005 Paris

Phone office:    +33 1 44 32 33 03
Phone lab:       +33 1 44 32 33 08
Fax office:      +33 1 44 32 34 34