This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: Bug in gsl_sf_gamma_inc_P_impl
- To: Teemu Ikonen <tpikonen at pcu dot helsinki dot fi>
- Subject: Re: Bug in gsl_sf_gamma_inc_P_impl
- From: Gerard Jungman <jungman at lanl dot gov>
- Date: Tue, 10 Oct 2000 12:33:30 -0600
- CC: gsl-discuss at sourceware dot cygnus dot com
- Organization: LANL T-8
- References: <Pine.OSF.4.20.0010082305090.19771-100000@tekno.helsinki.fi>
Teemu Ikonen wrote:
>
> I get GSL_EMAXITER when trying to calculate
> gsl_sf_gamma_inc_P_impl(100.0, 99.0+-eps, r)
> with gsl-0.6.
>
> Other x values work ok, only values very close to 99.0 cause the routine
> to fail.
This turned out to be a more general problem for the case a=x+1,
but only on some paths. I fixed the continued-fraction
evaluation on the failure path, so it should be ok now.
I just checked in the fixes.
Thanks.
--
G. Jungman