# HG changeset patch # User Rik # Date 1341772506 25200 # Node ID bcf86cc2f1ee16c04a99fdd111dcd7f0131a28f7 # Parent 460a3c6d8bf100be01f2879253e2bd9ced7b6067 Use xlgamma instead of lgamma in betaincinv for portability across systems. * lo-specfun.cc (betaincinv): Use xlgamma instead of lgamma. diff -r 460a3c6d8bf1 -r bcf86cc2f1ee liboctave/lo-specfun.cc --- a/liboctave/lo-specfun.cc Sun Jul 08 11:28:50 2012 -0700 +++ b/liboctave/lo-specfun.cc Sun Jul 08 11:35:06 2012 -0700 @@ -3287,7 +3287,7 @@ bool indx; double pp, prev, qq, r, s, sae = -37.0, sq, t, tx, value, w, xin, ycur, yprev; - double beta = lgamma (p) + lgamma (q) - lgamma (p + q); + double beta = xlgamma (p) + xlgamma (q) - xlgamma (p + q); bool err = false; fpu = pow (10.0, sae); value = y;