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

[Patrick Alken]

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/