Mercurial > octave
view liboctave/util/sparse-util.cc @ 21213:f7d1050b9b53
maint: Clean up various usages of #ifdef.
* randmtzig.c: Use #ifdef rather than just #if.
* EditControl.h : Use '#if ! defined' rather than '#ifndef' in
guard block to match Octave style.
* dialog.cc, settings-dialog.cc: Add FIXME notes about questionable use of #if
mechanism.
* file-editor-tab.cc, Array-b.cc, Array-ch.cc, Array-d.cc, Array-f.cc,
Array-i.cc, Array-s.cc, MatrixType.h, quit.h, lo-sysdep.cc, lo-cutils.c,
oct-alloc.h, oct-shlib.cc, sparse-sort.h: #define HAVE_XXX macros to 1, not
just empty but defined.
* octave-txt-lexer.h, octave-cmd.h, octave-preserve-stream-state.h, txt-eng.h,
zfstream.h, oct-conf-features.h, oct-conf.h: Use octave_ namespace prefix on
name of #define used to prevent multiple inclusion of headers.
* parser.h, webinfo.h, ov-oncleanup.h, op-int.h, display-available.h,
shared-fcns.h: Add #define guard to prevent multiple inclusion.
* quadcc.cc: use all capitals for #define MIN_CQUAD_HEAPSIZE.
* ov-intx.h: Add note that this file must not use guard #define.
* eigs-base.h, randmtzig.h: Write '! defined' rather than '!defined'.
* file-ops.cc, oct-sparse.h: Use parentheses around complex #if tests.
* oct-syscalls.cc, oct-base64.cc, statdefs.h: Indent #ifdef blocks correctly.
* oct-conf-post.in.h: Use "! defined". Define macros to 1, not just empty but
defined.
author | Rik <rik@octave.org> |
---|---|
date | Sat, 06 Feb 2016 18:22:32 -0800 |
parents | f7121e111991 |
children | 40de9f8f23a6 |
line wrap: on
line source
/* Copyright (C) 2005-2015 David Bateman Copyright (C) 1998-2005 Andy Adler 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 #include <stdio.h> #include <stdarg.h> #include "lo-error.h" #include "oct-sparse.h" #include "sparse-util.h" // FIXME: this overload is here due to API change in SuiteSparse (3.1 -> 3.2) void SparseCholError (int status, char *file, int line, char *message) { SparseCholError (status, file, line, message); } #ifdef HAVE_CHOLMOD void SparseCholError (int status, const char *file, int line, const char *message) { // Ignore CHOLMOD_NOT_POSDEF, since we handle that in Fchol as an // error or exit status. if (status != CHOLMOD_NOT_POSDEF) (*current_liboctave_warning_with_id_handler) ("Octave:cholmod-message", "warning %i, at line %i in file %s: %s", status, line, file, message); } #else void SparseCholError (int status, const char *, int, const char *) { } #endif int SparseCholPrint (const char *fmt, ...) { va_list args; va_start (args, fmt); int ret = gnulib::vfprintf (stderr, fmt, args); gnulib::fflush (stderr); va_end (args); return ret; } bool sparse_indices_ok (octave_idx_type *r, octave_idx_type *c, octave_idx_type nrows, octave_idx_type ncols, octave_idx_type nnz) { if (nnz > 0) { if (c[0] != 0) (*current_liboctave_error_handler) ("invalid sparse matrix: cidx[0] must be zero"); octave_idx_type jold = 0; for (octave_idx_type j = 1; j < ncols+1; j++) { if (c[j] < c[j-1]) (*current_liboctave_error_handler) ("invalid sparse matrix: cidx elements must appear in ascending order"); if (c[j] > nnz) (*current_liboctave_error_handler) ("invalid sparse matrix: cidx[%d] = %d " "exceeds number of nonzero elements", j, c[j]+1); if (c[j] != jold) { for (octave_idx_type i = jold+1; i < c[j]; i++) { if (r[i] < r[i-1]) (*current_liboctave_error_handler) ("invalid sparse matrix: ridx elements must appear " "in ascending order for each column"); if (r[i] >= nrows) (*current_liboctave_error_handler) ("invalid sparse matrix: ridx[%d] = %d out of range", i, r[i]+1); } jold = c[j]; } } } return true; }