Karim Belabas on Tue, 11 Dec 2018 15:45:35 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: bug in function airy_i in file hypergeom.c


* Peter [2018-12-11 15:13]:
> Dear developers,
> 
> I am absolutely new to the internals of the code. Therefore I can just
> describe the misbehavior of airy_i:
> 
> the test file for hypergeometric functions contains the instruction
> "airy(-100)[2]" which should return
> 
> airybi(-100) ~= 0.02427... but it returns 0.01213...
> 
> This might be caused by a malfunction of lines 74 - 80 in hypergeom.c

This is a known issue in the new code for Bessel(nu,z) function when
|z| -> oo committed yesterday. It's been caught up during the night by
our CI code, see

  http://pari.math.u-bordeaux.fr/buildlog.html

("BUG"'s all over the place); the cause is known as well, I'm working on
a fix right now. Stay tuned :-).

Cheers,

    K.B.
--
Karim Belabas, IMB (UMR 5251)  Tel: (+33) (0)5 40 00 26 17
Universite de Bordeaux         Fax: (+33) (0)5 40 00 21 23
351, cours de la Liberation    http://www.math.u-bordeaux.fr/~kbelabas/
F-33405 Talence (France)       http://pari.math.u-bordeaux.fr/  [PARI/GP]
`