Mercurial > octave-nkf
annotate liboctave/array/boolMatrix.cc @ 20603:f61c67865d9f
Don't return A for inv (A) when A is a singular Diagonal matrix (bug #46103).
* dDiagMatrix.cc (inverse), fDiagMatrix.cc (inverse): return octave_Inf for the
inverse of a diagonal element whose value is 0.
* test/diag-perm.tst: Add tests for new behavior.
| author | Rik <rik@octave.org> |
|---|---|
| date | Wed, 07 Oct 2015 07:41:44 -0700 |
| parents | a9574e3c6e9e |
| children |
| rev | line source |
|---|---|
| 2828 | 1 // Matrix manipulations. |
| 2 /* | |
| 3 | |
|
19731
4197fc428c7d
maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents:
19351
diff
changeset
|
4 Copyright (C) 1996-2015 John W. Eaton |
| 11523 | 5 Copyright (C) 2009-2010 VZLU Prague, a.s. |
| 2828 | 6 |
| 7 This file is part of Octave. | |
| 8 | |
| 9 Octave is free software; you can redistribute it and/or modify it | |
| 10 under the terms of the GNU General Public License as published by the | |
| 7016 | 11 Free Software Foundation; either version 3 of the License, or (at your |
| 12 option) any later version. | |
| 2828 | 13 |
| 14 Octave is distributed in the hope that it will be useful, but WITHOUT | |
| 15 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
| 16 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
| 17 for more details. | |
| 18 | |
| 19 You should have received a copy of the GNU General Public License | |
| 7016 | 20 along with Octave; see the file COPYING. If not, see |
| 21 <http://www.gnu.org/licenses/>. | |
| 2828 | 22 |
| 23 */ | |
| 24 | |
| 25 #ifdef HAVE_CONFIG_H | |
| 26 #include <config.h> | |
| 27 #endif | |
| 28 | |
| 3503 | 29 #include <iostream> |
| 2828 | 30 |
| 4669 | 31 #include "Array-util.h" |
|
19303
65554f5847ac
don't include oct-locbuf.h in header files unnecessarily
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
32 #include "boolMatrix.h" |
| 2828 | 33 #include "lo-error.h" |
| 34 #include "str-vec.h" | |
| 35 #include "mx-inlines.cc" | |
|
8774
b756ce0002db
split implementation and interface in mx-op-defs and MArray-defs
Jaroslav Hajek <highegg@gmail.com>
parents:
8743
diff
changeset
|
36 #include "mx-op-defs.h" |
| 2828 | 37 |
| 38 // boolMatrix class. | |
| 39 | |
| 40 bool | |
| 41 boolMatrix::operator == (const boolMatrix& a) const | |
| 42 { | |
| 43 if (rows () != a.rows () || cols () != a.cols ()) | |
| 44 return 0; | |
| 45 | |
|
20267
a9574e3c6e9e
Deprecate Array::length() and Sparse::length() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents:
19731
diff
changeset
|
46 return mx_inline_equal (numel (), data (), a.data ()); |
| 2828 | 47 } |
| 48 | |
| 49 bool | |
| 50 boolMatrix::operator != (const boolMatrix& a) const | |
| 51 { | |
| 52 return !(*this == a); | |
| 53 } | |
| 54 | |
| 55 boolMatrix& | |
| 5275 | 56 boolMatrix::insert (const boolMatrix& a, octave_idx_type r, octave_idx_type c) |
| 2828 | 57 { |
|
10351
5150ceb4dbb4
base charMatrix and boolMatrix on Array<char>
Jaroslav Hajek <highegg@gmail.com>
parents:
10158
diff
changeset
|
58 Array<bool>::insert (a, r, c); |
| 2828 | 59 return *this; |
| 60 } | |
| 61 | |
| 3203 | 62 // unary operations |
| 63 | |
| 64 boolMatrix | |
| 65 boolMatrix::operator ! (void) const | |
| 66 { | |
| 5275 | 67 octave_idx_type nr = rows (); |
| 68 octave_idx_type nc = cols (); | |
| 3203 | 69 |
| 70 boolMatrix b (nr, nc); | |
| 71 | |
| 5275 | 72 for (octave_idx_type j = 0; j < nc; j++) |
| 73 for (octave_idx_type i = 0; i < nr; i++) | |
| 3203 | 74 b.elem (i, j) = ! elem (i, j); |
| 75 | |
| 76 return b; | |
| 77 } | |
| 78 | |
| 79 // other operations | |
| 80 | |
| 6979 | 81 boolMatrix |
| 82 boolMatrix::diag (octave_idx_type k) const | |
| 83 { | |
|
10351
5150ceb4dbb4
base charMatrix and boolMatrix on Array<char>
Jaroslav Hajek <highegg@gmail.com>
parents:
10158
diff
changeset
|
84 return Array<bool>::diag (k); |
| 6979 | 85 } |
| 86 | |
|
9550
3d6a9aea2aea
refactor binary & bool ops in liboctave
Jaroslav Hajek <highegg@gmail.com>
parents:
9227
diff
changeset
|
87 MM_BOOL_OPS (boolMatrix, boolMatrix) |
|
3d6a9aea2aea
refactor binary & bool ops in liboctave
Jaroslav Hajek <highegg@gmail.com>
parents:
9227
diff
changeset
|
88 MS_BOOL_OPS (boolMatrix, bool) |
|
3d6a9aea2aea
refactor binary & bool ops in liboctave
Jaroslav Hajek <highegg@gmail.com>
parents:
9227
diff
changeset
|
89 SM_BOOL_OPS (bool, boolMatrix) |
|
9578
7dafdb8b062f
refactor comparison ops implementations
Jaroslav Hajek <highegg@gmail.com>
parents:
9550
diff
changeset
|
90 MM_CMP_OPS (boolMatrix, boolMatrix) |