[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] gsl_sf_beta_inc underflow error
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] gsl_sf_beta_inc underflow error |
Date: |
Tue, 31 Aug 2010 19:16:47 +0100 |
User-agent: |
Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) |
At Sat, 28 Aug 2010 17:41:53 +0200,
Thomas Tanner wrote:
>
> Hello,
>
> calling gsl_sf_beta_inc(0.5,101.5,0.999457);
>
> with GSL version 1.14 results in
>
> gsl: exp.c:545: ERROR: underflow
> Default GSL error handler invoked.
> Abort trap
>
> Other libraries (e.g., cephes) correctly return 1.0 as result.
>
Thanks for the bug report. I've removed the incorrect error code,
which came from an intermediate underflow that does not affect the
final result.
--
Brian Gough