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

[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.  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:

>>> cdf (k) = 1 - gammainc (l, floor (x (k)) + 1);

-- System Information:
Debian Release: 4.0
  APT prefers testing
  APT policy: (990, 'testing'), (500, 'testing-proposed-updates')
Architecture: amd64 (x86_64)
Shell:  /bin/sh linked to /bin/bash
Kernel: Linux 2.6.18-3-amd64
Locale: LANG=C, LC_CTYPE=it_IT at euro (charmap=ANSI_X3.4-1968) (ignored: LC_ALL set to C)

Versions of packages octave2.9 depends on:
ii  atlas3-base [liblapack.s 3.6.0-20.6      Automatically Tuned Linear Algebra
ii  fftw3                    3.1.2-1         library for computing Fast Fourier
ii  lapack3 [liblapack.so.3] 3.0.20000531a-6 library of linear algebra routines
ii  libc6                    2.3.6.ds1-13    GNU C Library: Shared libraries
ii  libgcc1                  1:4.1.1-21      GCC support library
ii  libgfortran1             4.1.1-21        Runtime library for GNU Fortran ap
ii  libglpk0                 4.11-2          linear programming kit (shared lib
ii  libhdf5-serial-1.6.5-0 [ 1.6.5-3         Hierarchical Data Format 5 (HDF5) 
ii  libncurses5              5.5-5           Shared libraries for terminal hand
ii  libpcre3                 6.7-1           Perl 5 Compatible Regular Expressi
ii  libreadline5             5.2-2           GNU readline and history libraries
ii  libstdc++6               4.1.1-21        The GNU Standard C++ Library v3
ii  libufsparse              1.2-7           collection of libraries for comput
ii  refblas3 [libblas.so.3]  1.2-8           Basic Linear Algebra Subroutines 3
ii  texinfo                  4.8.dfsg.1-4    Documentation system for on-line i
ii  zlib1g                   1:1.2.3-13      compression library - runtime

octave2.9 recommends no packages.

-- no debconf information