[Help-gsl] GSL Hankel normalization factor

[Darshana L Weerawarne]

Hi all,

I'm trying to use gsl dht to calculate hankel transform of a function. But
I'm having some troubles in determining the normalization factor. Following
is my psudo code.

    gsl_dht *tf;
    tf =  gsl_dht_alloc (4096);
    gsl_dht_init (tf, 1, 0.5);
    norm = (1/gsl_sf_bessel_zero_Jnu(1, 4096))^2;
    gsl_dht_apply (tf, uin1, uout);
    gsl_dht_apply (tf, uout, uin2);

Normalization does not satisfy the following condition.
    uin2 = uin1*norm

Could someone help?

Kind regards,
-- 
Darshana Lakmal Weerawarne