# HG changeset patch
# User jwe
# Date 1142962305 0
# Node ID 67118c88cee7d688410fcf06e7a51b4b2b0065ca
# Parent  8459ed2e2668b1d1333398cdd23ea1d71beb6c29
[project @ 2006-03-21 17:31:45 by jwe]

diff --git a/liboctave/ChangeLog b/liboctave/ChangeLog
--- a/liboctave/ChangeLog
+++ b/liboctave/ChangeLog
@@ -1,13 +1,21 @@
+2006-03-21  John W. Eaton  <jwe@octave.org>
+
+	* lo-specfun.cc (betainc): Use F77_XFCN instead of F77_FUNC for
+	call to xdbetai.
+
 2006-03-21  David Bateman  <dbateman@free.fr>
 
-	    * dSparse.cc (solve): Add argument singular_fallback, to allow
-	    fallback to QR solvers to be optional.
-	    * CSparse.cc (solve): ditto.
-	    * dSparse.h (solve): update declaration for new argument.
-	    * CSparse.h (solve): ditto.
-	    * sparse-dmsolve.cc (dmsolve): Use singular_fallback argument
-	    to bypass QR solvers when solving the well determined part of
-	    the problem.
+	* lo-specfun.cc (xlgamma): Require arg strictly greater than 0.
+	Use F77_XFCN instead of F77_FUNC for call to dlgams.
+
+	* dSparse.cc (solve): Add argument singular_fallback, to allow
+	fallback to QR solvers to be optional.
+	* CSparse.cc (solve): Ditto.
+	* dSparse.h (solve): update declaration for new argument.
+	* CSparse.h (solve): Ditto.
+	* sparse-dmsolve.cc (dmsolve): Use singular_fallback argument
+	to bypass QR solvers when solving the well determined part of
+	the problem.
 
 2006-03-17  John W. Eaton  <jwe@octave.org>
 
diff --git a/liboctave/lo-specfun.cc b/liboctave/lo-specfun.cc
--- a/liboctave/lo-specfun.cc
+++ b/liboctave/lo-specfun.cc
@@ -170,11 +170,11 @@
   double result;
   double sgngam;
 
-  if (x < 0)
+  if (x <= 0)
     (*current_liboctave_error_handler)
       ("xlgamma: argument must be nonnegative");
-
-  F77_XFCN (dlgams, DLGAMS, (x, result, sgngam));
+  else
+    F77_XFCN (dlgams, DLGAMS, (x, result, sgngam));
 
   return result;
 }
@@ -947,7 +947,7 @@
 betainc (double x, double a, double b)
 {
   double retval;
-  F77_FUNC (xdbetai, XDBETAI) (x, a, b, retval);
+  F77_XFCN (xdbetai, XDBETAI, (x, a, b, retval));
   return retval;
 }