Igor Schein on Mon, 15 Nov 1999 03:01:43 -0500


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

lngamma() bug(s)


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).

Curious thing here is that you see 3 abnormal indications
by varying the same thing.

Thanks

Igor