You can send me the two new source files
On 09/22/2014 03:13 PM, Raymond Rogers wrote:
Okay Patrick,
I cleared up my problems and the hyperg_U fixme's . There are
three small (10 lines/each or so) blocks of inserted code and the
modified test file with the fixme's commented out.
It needs to be reviewed for style and error calculations. I am not
a numerical analyst so I just allowed the error margins fall through.
How would you like these delivered? Two new source files, a
pointer to one of my repositories, or just a diff file?
Ray
On 09/19/2014 12:19 PM, Patrick Alken wrote:
Thanks for your report. If you manage to fix some of these cases, we
are happy to accept a patch
On 09/19/2014 10:17 AM, Patrick Alken wrote:
URL:
<http://savannah.gnu.org/bugs/?43258>
Summary: confluent hypergeometric functions
Project: GNU Scientific Library
Submitted by: psa
Submitted on: Fri 19 Sep 2014 04:17:55 PM GMT
Category: Runtime error
Severity: 3 - Normal
Operating System:
Status: None
Assigned to: None
Open/Closed: Open
Release:
Discussion Lock: Any
_______________________________________________________
Details:
From raymond dot rogers72 =at= gmail dot com
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.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?43258>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
