[Pkg-octave-devel] Bug#493869: Bug#493869: octave3.0: gaminv and gamcdf crash for big values of parameter A

[Thomas Weber]

Am Dienstag, den 05.08.2008, 15:49 +0200 schrieb Sebastien Villemot:
> Package: octave3.0
> Version: 1:3.0.1-6
> Severity: normal
> 
> gaminv and gamcdf function crash with the following parameters:

> octave> gaminv(0.5,600,1)
>  ***MESSAGE FROM ROUTINE D9LGIT IN LIBRARY SLATEC.
>  ***FATAL ERROR, PROG ABORTED, TRACEBACK REQUESTED
>  *  NO CONVERGENCE IN 200 TERMS OF CONTINUED FRACTION
>  *  ERROR NUMBER = 3
>  *   
>  ***END OF MESSAGE

> octave> gamcdf(600, 600, 1)
>  ***MESSAGE FROM ROUTINE D9LGIT IN LIBRARY SLATEC.
>  ***FATAL ERROR, PROG ABORTED, TRACEBACK REQUESTED
>  *  NO CONVERGENCE IN 200 TERMS OF CONTINUED FRACTION
>  *  ERROR NUMBER = 3
>  *   
>  ***END OF MESSAGE

That's not a crash, is it? The functions inform you that they didn't
converge. 


> Note that if you replace 600 by 500, the problem doesn't occur.

See
http://bugs.debian.org/cgi-bin/bugreport.cgi?bug=418158

d9lgit.f needs an improvement, either directly or by completely
replacing it with something better.

	Thomas