This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: BUG in gsl_sf_bessel_zero_J1_impl
- To: "Gerard Vermeulen (CRTBT)" <gvermeul at labs dot polycnrs-gre dot fr>
- Subject: Re: BUG in gsl_sf_bessel_zero_J1_impl
- From: Gerard Jungman <jungman at lanl dot gov>
- Date: Wed, 25 Oct 2000 12:16:35 -0600
- CC: gsl-discuss at sources dot redhat dot com
- Organization: LANL T-8
- References: <39F70B19.C333DABA@labs.polycnrs-gre.fr>
"Gerard Vermeulen (CRTBT)" wrote:
>
> Concerns version 0.6:
>
> result.val is about 0.3 after the call
> gsl_sf_bessel_zero_J1_impl(0, &result);
> instead of 0, as I would expect.
>
> Higher zeros are OK.
Yes, this was never accounted for, under the assumption
that the behaviour for s=0 was "undefined" or something.
Anyway, I agree that was silly, so I fixed it to
return 0.0 for s=0.
> I also point out that gsl_sf_bessel_zero_J1_e & friends are mentioned
> in the documentation, but are not implemented, yet.
gsl_sf_bessel_zero_J1_e() was inadvertently left unimplemented;
I fixed that. I'm not sure what you mean by "friends", as
both the gsl_sf_bessel_zero_J0_e() and gsl_sf_bessel_zero_Jnu_e()
implementations are present.
The fixes are now checked into CVS.
Thanks for the report!
--
G. Jungman