help-gsl
[Top][All Lists]
Advanced

[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-----



reply via email to

[Prev in Thread] Current Thread [Next in Thread]