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
