[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
## [Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values

**From**: |
Daniel Schury |

**Subject**: |
[Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values |

**Date**: |
Fri, 10 May 2013 11:43:59 +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
Hi,
I've got some problems with the usage of the above mentioned function.
My usage is
"double part22 = gsl_sf_hyperg_1F1(-p/2.0, 1.0/2.0,
- -1.0/2.0*((E1-x)/dE)^2.0);"
This works for most values, but fails for some. With p = 2.09, E1 = 8, x
= 26.8 and dE = 0.5 my program crashes.
I'm using the library shipped with Origin 9.0 in C, Windows XP, 32 bit.
I'm not completely into the mathematics about the hypergeometric
function, but mathematica for example can compute this function with the
same values, so for me it seems not like a mathematical, but a
programming problem.
Can somebody confirm this problem maybe on a different platform or has
an idea how to trace the problem?
Many thanks in advance.
- --
Freundliche Grüße
Daniel Schury
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.12 (GNU/Linux)
Comment: Using GnuPG with undefined - http://www.enigmail.net/
iEYEARECAAYFAlGMwVwACgkQk7XSfqVgB0oIsQCfU42V6cW45UbT+o20pdehy+Yc
70QAoLh3P/AHi3KAPXYlLM/PYU71FwAs
=8vxP
-----END PGP SIGNATURE-----

**[Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values**,
*Daniel Schury* **<=**
**Re: [Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values**, *Sam Mason*, `2013/05/10`
**Re: [Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values**, *Daniel Schury*, `2013/05/10`
**Re: [Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values**, *Sam Mason*, `2013/05/10`
**Re: [Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values**, *Daniel Schury*, `2013/05/10`
**Re: [Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values**, *Sam Mason*, `2013/05/10`
**Re: [Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values**, *Daniel Schury*, `2013/05/10`
**[Help-gsl] gsl_sf_hyperg_1F1 need better(?) values of the exp()-function**, *Daniel Schury*, `2013/05/13`