[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Bug report: gsl_sf_hyperg_2F1
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Bug report: gsl_sf_hyperg_2F1 |
Date: |
Tue, 20 Sep 2011 17:27:05 +0100 |
User-agent: |
Wanderlust/2.15.6 (Almost Unreal) Emacs/23.3 Mule/6.0 (HANACHIRUSATO) |
At Mon, 5 Sep 2011 16:46:44 -0400,
Robert L. Wolpert wrote:
>
> Hi, GNU.
>
> The Gauss hypergeometric function 2F1 is described in detail in Abramowitz &
> Stegan (1964) eqn. 15.1.1, page 556. It is implemented in the gsl library as
> "gsl_sf_hyperg_2F1". The implementation frequently returns "NaN" or "Inf"
> for domain values where 2F1 is well-defined and easy to calculate. I hope
> you will address this in an upcoming version, because your implementations
> are widely used and respected.
Thanks for the bug report. I've added a test case for the values in
your example program, although I don't have a fix for it right now.
>
> Specifically:
>
> GSL Version: 1.3
Not sure if the above version is a typo, but 1.3 is very old. The
current release is 1.15.
--
Brian Gough