Mercurial > octave-nkf
changeset 5700:67118c88cee7
[project @ 2006-03-21 17:31:45 by jwe]
| author | jwe |
|---|---|
| date | Tue, 21 Mar 2006 17:31:45 +0000 |
| parents | 8459ed2e2668 |
| children | 252b6754e545 |
| files | liboctave/ChangeLog liboctave/lo-specfun.cc |
| diffstat | 2 files changed, 20 insertions(+), 12 deletions(-) [+] |
line wrap: on
line diff
--- a/liboctave/ChangeLog Tue Mar 21 16:42:54 2006 +0000 +++ b/liboctave/ChangeLog Tue Mar 21 17:31:45 2006 +0000 @@ -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>
--- a/liboctave/lo-specfun.cc Tue Mar 21 16:42:54 2006 +0000 +++ b/liboctave/lo-specfun.cc Tue Mar 21 17:31:45 2006 +0000 @@ -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; }