[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] Confluent Hypergeometric problems
From: |
Raymond Rogers |
Subject: |
[Bug-gsl] Confluent Hypergeometric problems |
Date: |
Fri, 19 Sep 2014 12:08:00 -0400 |
User-agent: |
Mozilla/5.0 (X11; Linux x86_64; rv:31.0) Gecko/20100101 Thunderbird/31.1.1 |
If you will check you will find that
hyperg_U.c fails for U(a,a-1,1) a>5
hyperg_U.c fails for U(2,b,1) for b< -2 (with occasional holes due to
rounding)
These can be checked through gsl-shell, gsl-shell-gui, or I have a
command line interface.
I believe there is a choke point for the error at the start of
hyperg_U_series. There are a couple of sources to it and a couple of
crashes afterwards. A fix should be 13.2.8 in DLMF. I would like to
discuss this with somebody who has experience with this code before I
spend much time on the fix.
This should also fix a couple of "unimplemented cases"
I also have a spreadsheet and program comparing test cases from GSL,
mpmath, and DAMath; for Confluent Hypergeometric tests. Not complete
yet but it 700 (or so) instances.
Ray
--
The primary use of conversation is to satisfy the impulse to talk
George Santanyana
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Bug-gsl] Confluent Hypergeometric problems,
Raymond Rogers <=