Re: error in test beta.m

[Tatsuro MATSUOKA]

Hello

I have a look at beta.m

retval = real (exp (gammaln (a) + gammaln (b) - gammaln (a+b)));

In the case of the test in failure gammaln(0) is the origin of the error,
octave.exe:43> gammaln(0)
 ***MESSAGE FROM ROUTINE DGAMMA IN LIBRARY SLATEC.
 ***FATAL ERROR, PROG ABORTED, TRACEBACK REQUESTED
 *  X IS 0
 *  ERROR NUMBER = 4
 *
 ***END OF MESSAGE

 ***JOB ABORT DUE TO FATAL ERROR.
0          ERROR MESSAGE SUMMARY
 LIBRARY    SUBROUTINE MESSAGE START             NERR     LEVEL     COUNT
 SLATEC     DGAMMA     X IS 0                       4         2         5

error: lgamma: exception encountered in Fortran subroutine dlgams_

Hmmm!

How about gammaln(0) on the other platform ?

Regards

Tatsuro
--- Tatsuro MATSUOKA wrote:

> Hello
> 
> I have encountered the following error in development source.
> (last change  Sun, 25 Jul 2010 22:27:08 -0700 d9147775da54)
> 
> Enviornmen:s MinGW32 GCC 4.5.0, windows xp pro.
> 
> ./run-octave
>  :
> octave.exe:3> test beta verbose
> octave.exe:1> test beta
>  ***MESSAGE FROM ROUTINE DGAMMA IN LIBRARY SLATEC.
>  ***FATAL ERROR, PROG ABORTED, TRACEBACK REQUESTED
>  *  X IS 0
>  *  ERROR NUMBER = 4
>  *
>  ***END OF MESSAGE
> 
>  ***JOB ABORT DUE TO FATAL ERROR.
> 0          ERROR MESSAGE SUMMARY
>  LIBRARY    SUBROUTINE MESSAGE START             NERR     LEVEL     COUNT
>  SLATEC     DGAMMA     X IS 0                       4         2         1
> 
>   ***** test
>  a = 0.25 + (0:5) * 0.5;
>  tol = 10 * max (a) * eps;
>  assert (zeros (size (a)), beta (a, -a), tol)
>  assert (zeros (size (a)), beta (-a, a), tol)
> !!!!! test failed
> exception encountered in Fortran subroutine dlgams_
> *******************
> 
> To make clear deeply I probably have to use ./run-octave -g
> I have tried 'br dlgams_' and tcatch catch.
> 
> I do not have enough knowledge of gdb.
> 
>  
> (gdb) br dlgams_
> No symbol table is loaded.  Use the "file" command.
> Make breakpoint pending on future shared library load? (y or [n]) y
> Breakpoint 1 (dlgams_) pending.
> 
> (gdb) run
> <snip>
> 
> (gdb) br dlgams_
> No symbol table is loaded.  Use the "file" command.
> Make breakpoint pending on future shared library load? (y or [n]) y
> 
> octave.exe:1> a = 0.25 + (0:5) * 0.5;
> octave.exe:2> beta (a, -a)
> Breakpoint 1, 0x00501ebd in dlgams_ ()
>    from 
> d:\usr\Tatsu\mingwhome\octaves\OctBuild\hg-devel\libcruft\.libs\libcruft-3-3-51+.dll
> (gdb) tcatch catch
> Note: breakpoint 2 also set at pc 0x640d9278.
> Temporary catchpoint 3 (catch)
> (gdb) continue
> Continuing.
> <snip>
> Breakpoint 1, 0x00501ebd in dlgams_ ()
>    from 
> d:\usr\Tatsu\mingwhome\octaves\OctBuild\hg-devel\libcruft\.libs\libcruft-3-3-51+.dll
> (gdb) continue
> Continuing.
>  ***MESSAGE FROM ROUTINE DGAMMA IN LIBRARY SLATEC.
>  ***FATAL ERROR, PROG ABORTED, TRACEBACK REQUESTED
>  *  X IS 0
>  *  ERROR NUMBER = 4
>  *
>  ***END OF MESSAGE
> 
>  ***JOB ABORT DUE TO FATAL ERROR.
> 0          ERROR MESSAGE SUMMARY
>  LIBRARY    SUBROUTINE MESSAGE START             NERR     LEVEL     COUNT
>  SLATEC     DGAMMA     X IS 0                       4         2         2
> 
> error: lgamma: exception encountered in Fortran subroutine dlgams_
> error: evaluating argument list element number 1
> error: evaluating argument list element number 1
> error: called from:
> error:  
>
d:\\usr\\Tatsu\\mingwhome\\octaves\\OctBuild\\hg-devel\\..\\..\\hg\\octave-work\\scripts\\specfun\beta.m
> at line 54, column 10
> 
> *******************
> 
> Please give me some suggestions.
>  
> Regards
> 
> Tatsuro
> 
> --------------------------------------
> Get the new Internet Explorer 8 optimized for Yahoo! JAPAN
> http://pr.mail.yahoo.co.jp/ie8/
> 


--------------------------------------
Get the new Internet Explorer 8 optimized for Yahoo! JAPAN
http://pr.mail.yahoo.co.jp/ie8/