Mercurial > octave-nkf
view liboctave/boolMatrix.cc @ 6980:e5d51f9fa1c9
[project @ 2007-10-08 21:52:04 by jwe]
author | jwe |
---|---|
date | Mon, 08 Oct 2007 21:52:04 +0000 |
parents | 2883ea1c5c18 |
children | 93c65f2a5668 |
line wrap: on
line source
// Matrix manipulations. /* Copyright (C) 1996, 1997 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, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. */ #ifdef HAVE_CONFIG_H #include <config.h> #endif #include <iostream> #include "Array-util.h" #include "lo-error.h" #include "str-vec.h" #include "mx-base.h" #include "mx-inlines.cc" // boolMatrix class. bool boolMatrix::operator == (const boolMatrix& a) const { if (rows () != a.rows () || cols () != a.cols ()) return 0; return mx_inline_equal (data (), a.data (), length ()); } bool boolMatrix::operator != (const boolMatrix& a) const { return !(*this == a); } boolMatrix& boolMatrix::insert (const boolMatrix& a, octave_idx_type r, octave_idx_type c) { Array2<bool>::insert (a, r, c); return *this; } // unary operations boolMatrix boolMatrix::operator ! (void) const { octave_idx_type nr = rows (); octave_idx_type nc = cols (); boolMatrix b (nr, nc); for (octave_idx_type j = 0; j < nc; j++) for (octave_idx_type i = 0; i < nr; i++) b.elem (i, j) = ! elem (i, j); return b; } // other operations boolMatrix boolMatrix::diag (void) const { return diag (0); } boolMatrix boolMatrix::diag (octave_idx_type k) const { octave_idx_type nnr = rows (); octave_idx_type nnc = cols (); if (k > 0) nnc -= k; else if (k < 0) nnr += k; boolMatrix d; if (nnr > 0 && nnc > 0) { octave_idx_type ndiag = (nnr < nnc) ? nnr : nnc; d.resize (ndiag, 1); if (k > 0) { for (octave_idx_type i = 0; i < ndiag; i++) d.xelem (i) = elem (i, i+k); } else if (k < 0) { for (octave_idx_type i = 0; i < ndiag; i++) d.xelem (i) = elem (i-k, i); } else { for (octave_idx_type i = 0; i < ndiag; i++) d.xelem (i) = elem (i, i); } } else (*current_liboctave_error_handler) ("diag: requested diagonal out of range"); return d; } // FIXME Do these really belong here? Maybe they should be // in a base class? boolMatrix boolMatrix::all (int dim) const { MX_ALL_OP (dim); } boolMatrix boolMatrix::any (int dim) const { MX_ANY_OP (dim); } MM_CMP_OPS (boolMatrix, , boolMatrix, ) /* ;;; Local Variables: *** ;;; mode: C++ *** ;;; End: *** */