From: Stefan Gerlach <stefan.gerlach@uni-konstanz.de>
To: gsl-discuss@sourceware.org
Subject: possible bug in gsl_sf_hyperg_U
Date: Thu, 03 Apr 2008 12:45:00 -0000 [thread overview]
Message-ID: <200804031110.22962.stefan.gerlach@uni-konstanz.de> (raw)
[-- Attachment #1: Type: text/plain, Size: 452 bytes --]
Hi,
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
AFAIK the hypergeometricU is well defined for z> 0.0.
I get the same error with gsl-1.11. Any ideas?
--
Stefan Gerlach
[-- Attachment #2: Type: application/pgp-signature, Size: 189 bytes --]
reply other threads:[~2008-04-03 12:45 UTC|newest]
Thread overview: [no followups] expand[flat|nested] mbox.gz Atom feed
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=200804031110.22962.stefan.gerlach@uni-konstanz.de \
--to=stefan.gerlach@uni-konstanz.de \
--cc=gsl-discuss@sourceware.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).