[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values
|
From: |
Daniel Schury |
|
Subject: |
Re: [Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values |
|
Date: |
Fri, 10 May 2013 15:23:58 +0200 |
|
User-agent: |
Mozilla/5.0 (X11; Linux x86_64; rv:17.0) Gecko/20130404 Thunderbird/17.0.5 |
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
Well thanks a lot! So the problem is the calculation of the exponential
function, namely with gsl_sf_exp_mult_err_e().
Am 10.05.2013 14:38, schrieb Sam Mason:
> On 10 May 2013 13:07, Daniel Schury <address@hidden> wrote:
>> Am 10.05.2013 13:58, schrieb Sam Mason:
>>> gsl_sf_hyperg_1F1(-1.045, 0.5, -706.88)
>> Ok, now I know what you mean. Does not change anything, program crashes.
>
> Huh, I should have tried that! I get:
>
> gsl: exp.c:254: ERROR: overflow
>
> not sure if I can help much more as I don't know about the
> Hypergeometric functions.
>
> Sorry,
>
> Sam
>
- --
Freundliche Grüße
Daniel Schury
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.12 (GNU/Linux)
Comment: Using GnuPG with undefined - http://www.enigmail.net/
iEYEARECAAYFAlGM9O4ACgkQk7XSfqVgB0pEnwCfbO/tuKr9EMWbIRS6ZFfuFZr0
I5cAoK9oaWcUlMK0136Ohyjoz7r23WGC
=GE25
-----END PGP SIGNATURE-----