Mercurial > octave-nkf
view src/chol.cc @ 636:fae2bd91c027
[project @ 1994-08-23 18:39:50 by jwe]
author | jwe |
---|---|
date | Tue, 23 Aug 1994 18:39:50 +0000 |
parents | 8e4e7e5f307e |
children | 0a81458ef677 |
line wrap: on
line source
// f-chol.cc -*- C++ -*- /* Copyright (C) 1993, 1994 John W. Eaton This file is part of Octave. Octave is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2, or (at your option) any later version. Octave is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with Octave; see the file COPYING. If not, write to the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. */ #ifdef HAVE_CONFIG_H #include "config.h" #endif #include "dbleCHOL.h" #include "CmplxCHOL.h" #include "tree-const.h" #include "user-prefs.h" #include "gripes.h" #include "error.h" #include "utils.h" #include "help.h" #include "defun-dld.h" DEFUN_DLD ("chol", Fchol, Schol, 2, 1, "R = chol (X): cholesky factorization") { Octave_object retval; if (args.length () != 2 || nargout > 1) { print_usage ("chol"); return retval; } tree_constant arg = args(1); int nr = arg.rows (); int nc = arg.columns (); if (empty_arg ("chol", nr, nc) < 0) return retval; if (arg.is_real_type ()) { Matrix m = arg.matrix_value (); if (! error_state) { int info; CHOL fact (m, info); if (info != 0) error ("chol: matrix not positive definite"); else retval = fact.chol_matrix (); } } else if (arg.is_complex_type ()) { ComplexMatrix m = arg.complex_matrix_value (); if (! error_state) { int info; ComplexCHOL fact (m, info); if (info != 0) error ("chol: matrix not positive definite"); else retval = fact.chol_matrix (); } } else { gripe_wrong_type_arg ("chol", arg); } return retval; } /* ;;; Local Variables: *** ;;; mode: C++ *** ;;; page-delimiter: "^/\\*" *** ;;; End: *** */