[Pkg-octave-devel] Bug#418158: Bug#418158: octave2.9: poissinv error

[Thomas Weber]

Am Samstag, den 07.04.2007, 14:42 +0200 schrieb Francesco Potorti`:
> Package: octave2.9
> Version: 2.9.9-8etch1
> Severity: normal
> 
> octave:70> poisscdf(10,10)
> ans =  0.58304
> octave:71> poissinv(10,10)
> ans = NaN
> 
> The above should be an error.  

Housekeeping:
Is my understanding correct that the above is actually correct and the
remaining error is the problem below? See also the linked discussion on
the bug at octave mailingllist.

> Additionally, poisscdf(x,x) works for
> x=100, but does not for x=1000:
> 
> octave:72> poisscdf(1000,1000)
>  ***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
>  
>  ***JOB ABORT DUE TO FATAL ERROR.
> 0          ERROR MESSAGE SUMMARY
>  LIBRARY    SUBROUTINE MESSAGE START             NERR     LEVEL     COUNT
>  SLATEC     D9LGIT     NO CONVERGENCE IN 20         3         2         1
> 
> error: exception encountered in Fortran subroutine xgammainc_
> error: in /usr/share/octave/2.9.9/m/statistics/distributions/poisscdf.m near line 58, column 14:

Thanks
	Thomas