This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: gsl_sf_psi_1piy
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Hiroshi Imamura <hima at cmt dot is dot tohoku dot ac dot jp>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Wed, 29 May 2002 19:13:35 +0100 (BST)
- Subject: Re: gsl_sf_psi_1piy
- References: <20020528184449.36C4.HIMA@cmt.is.tohoku.ac.jp>
Hiroshi Imamura writes:
> Hi,
> gsl_sf_psi_1piy() is an even function. (see 6.3.10 of Abramowitz
> and Stegun)
> The following changes should be made in psi.c:
> line 491
> const double lny = log(y); /* FIXME: y < 0 ?? */
> --->
> const double lny = log(ay);
> line 500
> const double lny = log(y);
> --->
> const double lny = log(ay);
Thanks. I'll add those changes for the next release.
regards
Brian Gough