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

[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
 
 ***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         8

error: exception encountered in Fortran subroutine xgammainc_
error: evaluating assignment expression near line 54, column 15
error: evaluating if command near line 53, column 5
error: evaluating if command near line 52, column 3
error: called from `gamcdf' in file `/usr/share/octave/3.0.1/m/statistics/distributions/gamcdf.m'
error: evaluating binary operator `-' near line 73, column 37
error: evaluating binary operator `./' near line 73, column 42
error: evaluating assignment expression near line 73, column 13
error: evaluating for command near line 72, column 5
error: evaluating if command near line 57, column 3
error: called from `gaminv' in file `/usr/share/octave/3.0.1/m/statistics/distributions/gaminv.m'
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
 
 ***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         9

error: exception encountered in Fortran subroutine xgammainc_
error: evaluating assignment expression near line 54, column 15
error: evaluating if command near line 53, column 5
error: evaluating if command near line 52, column 3
error: called from `gamcdf' in file `/usr/share/octave/3.0.1/m/statistics/distributions/gamcdf.m'
octave> 

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

-- System Information:
Debian Release: lenny/sid
  APT prefers testing
  APT policy: (500, 'testing')
Architecture: i386 (i686)

Kernel: Linux 2.6.25-2-686 (SMP w/1 CPU core)
Locale: LANG=fr_FR.UTF-8, LC_CTYPE=fr_FR.UTF-8 (charmap=UTF-8)
Shell: /bin/sh linked to /bin/bash

Versions of packages octave3.0 depends on:
ii  libatlas3gf-base [libl 3.6.0-21.5        Automatically Tuned Linear Algebra
ii  libblas3gf [libblas.so 1.2-2             Basic Linear Algebra Subroutines 3
ii  libc6                  2.7-10            GNU C Library: Shared libraries
ii  libcurl3-gnutls        7.18.2-5          Multi-protocol file transfer libra
ii  libfftw3-3             3.1.2-3.1         library for computing Fast Fourier
ii  libgcc1                1:4.3.1-2         GCC support library
ii  libgfortran3           4.3.1-2           Runtime library for GNU Fortran ap
ii  libglpk0               4.29-2            linear programming kit with intege
ii  libhdf5-serial-1.6.6-0 1.6.6-4           Hierarchical Data Format 5 (HDF5) 
ii  liblapack3gf [liblapac 3.1.1-1           library of linear algebra routines
ii  libncurses5            5.6+20080713-1    shared libraries for terminal hand
ii  libpcre3               7.6-2.1           Perl 5 Compatible Regular Expressi
ii  libqhull5              2003.1-9          calculate convex hulls and related
ii  libreadline5           5.2-3             GNU readline and history libraries
ii  libstdc++6             4.3.1-2           The GNU Standard C++ Library v3
ii  libsuitesparse-3.1.0   3.1.0-3           collection of libraries for comput
ii  texinfo                4.11.dfsg.1-4     Documentation system for on-line i
ii  zlib1g                 1:1.2.3.3.dfsg-12 compression library - runtime

Versions of packages octave3.0 recommends:
ii  gnuplot                       4.2.2-1.1  A command-line driven interactive 
ii  libatlas3gf-base              3.6.0-21.5 Automatically Tuned Linear Algebra

Versions of packages octave3.0 suggests:
pn  octave3.0-doc                 <none>     (no description available)
ii  octave3.0-emacsen             1:3.0.1-6  Emacs support for the GNU Octave l
ii  octave3.0-headers             1:3.0.1-6  header files for the GNU Octave la
ii  octave3.0-htmldoc             1:3.0.1-6  HTML documentation on the GNU Octa
ii  octave3.0-info                1:3.0.1-6  GNU Info documentation on the GNU 

-- no debconf information