Igor Schein on Thu, 18 Nov 1999 17:10:32 -0500


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

Re: lngamma() bug(s)


On Mon, Nov 15, 1999 at 03:01:43AM -0500, Igor Schein wrote:
> Hi,
> 
> ? \p9
>    realprecision = 9 significant digits
> ? lngamma(lngamma(lngamma(1/4)))
>   ***   bug in wr_float, please report
> ? \p19
>    realprecision = 19 significant digits
> ? lngamma(lngamma(lngamma(1/4)))
>   ***   lbot>ltop in gerepile.
> 
> ? n=1/4;while(n=lngamma(n),print(n))
> 1.288022524698077457
> -0.1060666245741884610
>   ***   segmentation fault: bug in GP (please report).

All these have been fixed.  However, here's another old bug:

? \p38
   realprecision = 38 significant digits
? lngamma(-7.492741256813943588561666778-0.4185643959353916901644072949*I)
  ***   bug in GP (Segmentation Fault), please report

Thanks

Igor