Mercurial > octave-nkf
view src/matherr.c @ 8969:3ecbc236e2e0
Have sparse LU return permutation matrices rather than sparse matrices.
This could well impact user code. It'd be interesting to see if there
is any actual fall-out... Quite often, the permutation matrices are
applied to *dense* vectors. Returning permutation matrices rather
than sparse matrices is a slight performance enhancement, but likely
lost in the noise.
author | Jason Riedy <jason@acm.org> |
---|---|
date | Tue, 10 Mar 2009 21:54:49 -0400 |
parents | 2eb392d058bb |
children | cd96d29c5efa |
line wrap: on
line source
/* Copyright (C) 1997, 2005, 2007 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 3 of the License, 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, see <http://www.gnu.org/licenses/>. */ #ifdef HAVE_CONFIG_H #include <config.h> #endif #if defined (EXCEPTION_IN_MATH) #include "lo-math.h" int matherr (struct exception *x) { /* Possibly print our own message someday. Should probably be user-switchable. */ switch (x->type) { case DOMAIN: case SING: case OVERFLOW: case UNDERFLOW: case TLOSS: case PLOSS: default: break; } /* But don't print the system message. */ return 1; } #endif /* ;;; Local Variables: *** ;;; mode: C *** ;;; End: *** */