Re: [Bug-gsl] [bug #43258] confluent hypergeometric functions

[Raymond Rogers]

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/

--
The primary use of conversation is to satisfy the impulse to talk
George Santanyana