[Octave-bug-tracker] [bug #50919] Beta function throws FORTRAN errors (Windows only)

[John W. Eaton]

Follow-up Comment #6, bug #50919 (project octave):

I don't think it's fp87 vs sse that's causing the trouble here.  The beta
function ultimately calls the rc_gamma function in
liboctave/numeric/lo-specfun.cc.  That function uses lgamma_r if it is
available, otherwise it uses the SLATEC function DLGAMS (or ALGAMS).  I can
reproduce the problem on my Debian system if I force Octave to use the SLATEC
function.

On Windows systems, lgamma is available, but lgamma_r is not, so it uses the
SLATEC function.  I pushed the following changeset so that Octave will use
lgamma if it is available, and only if both lgamma_r and lgamma are
unavailable, will it fall back to using the SLATEC function.  I'll test with
an mxe-octave build when buildbot does its next build.

http://hg.savannah.gnu.org/hgweb/octave/rev/5296657cde44

Maybe it would be good to have a proper replacement for lgamma_r.  If so, I
think the right place for it is gnulib, not Octave itself.

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?50919>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/