[Pkg-octave-devel] Bug#493869: octave3.0: gaminv and gamcdf crash for big values of parameter A
Sebastien Villemot
sebastien.villemot at ens.fr
Tue Aug 5 13:49:06 UTC 2008
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
More information about the Pkg-octave-devel
mailing list