PreviousIndexLast

[bug-gsl] Re: possible bug in gsl_sf_hyperg_U (2008-04-06)

From: Brian Gough <bjg@network-theory.co.uk>
Subject: Re: possible bug in gsl_sf_hyperg_U
Date: Sun, 06 Apr 2008 20:49:04 +0100
To: Stefan.Gerlach@...............
Cc: bug-gsl@gnu.org

At Thu, 3 Apr 2008 11:22:01 +0200,
Stefan Gerlach wrote:
> I have a problem using gsl_sf_hyperg_U with certain arguments.
> running this example
>
> for(z=3.0;z>0.0;z-=0.1)
> printf("%g %g\n",z,gsl_sf_hyperg_U(-2.0,4.0,z));
>
> i get
>
> ...
> 2.2 2.84
> 2.1 3.41
> 2 4
> gsl: psi.c:381: ERROR: domain error
> Default GSL error handler invoked.
> Aborted
>

Thanks for the bug report. I have logged it in the bug database and
added a test case for it to the regression tests.

--
Brian Gough

Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/