Re: error in test beta.m

[Tatsuro MATSUOKA]

Sorry my poor Engish 

--- Tatsuro MATSUOKA  wrote:

> 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 ?
==> What happens when gammaln(0) on other platforms ?



> 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/
> 


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