octave: liboctave/array/dDiagMatrix.cc annotate

annotate liboctave/array/dDiagMatrix.cc @ 27919:1891570abac8

update Octave Project Developers copyright for the new year In files that have the "Octave Project Developers" copyright notice, update for 2020.

author	John W. Eaton <jwe@octave.org>
date	Mon, 06 Jan 2020 22:29:51 -0500
parents	b442ec6dda5c
children	bd51beb6205e

rev	line source
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	1 /*
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	2
27919 1891570abac8 update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 27918 diff changeset	3 Copyright (C) 1994-2020 The Octave Project Developers
27918 b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 27093 diff changeset	4
b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 27093 diff changeset	5 See the file COPYRIGHT.md in the top-level directory of this distribution
b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 27093 diff changeset	6 or <https://octave.org/COPYRIGHT.html/>.
b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 27093 diff changeset	7
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	8
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	9 This file is part of Octave.
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	10
24534 194eb4bd202b maint: Update punctuation for GPL v3 license text. Rik <rik@octave.org> parents: 23807 diff changeset	11 Octave is free software: you can redistribute it and/or modify it
22755 3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	12 under the terms of the GNU General Public License as published by
24534 194eb4bd202b maint: Update punctuation for GPL v3 license text. Rik <rik@octave.org> parents: 23807 diff changeset	13 the Free Software Foundation, either version 3 of the License, or
22755 3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	14 (at your option) any later version.
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	15
22755 3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	16 Octave is distributed in the hope that it will be useful, but
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	17 WITHOUT ANY WARRANTY; without even the implied warranty of
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	18 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	19 GNU General Public License for more details.
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	20
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	21 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5307 diff changeset	22 along with Octave; see the file COPYING. If not, see
24534 194eb4bd202b maint: Update punctuation for GPL v3 license text. Rik <rik@octave.org> parents: 23807 diff changeset	23 <https://www.gnu.org/licenses/>.
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	24
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	25 */
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	26
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21723 diff changeset	27 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21202 diff changeset	28 # include "config.h"
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	29 #endif
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	30
25438 cb1606f78f6b prefer <istream>, <ostream>, or <iosfwd> to <iostream> where possible John W. Eaton <jwe@octave.org> parents: 25054 diff changeset	31 #include <ostream>
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	32
4669 334a27c8f453 [project @ 2003-11-26 07:02:42 by jwe] jwe parents: 4513 diff changeset	33 #include "Array-util.h"
1368 33bb7975f866 [project @ 1995-09-06 00:12:48 by jwe] jwe parents: 1360 diff changeset	34 #include "lo-error.h"
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	35 #include "mx-base.h"
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	36 #include "mx-inlines.cc"
1650 23aa282707e8 [project @ 1995-12-20 06:53:12 by jwe] jwe parents: 1627 diff changeset	37 #include "oct-cmplx.h"
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	38
1360 7eb93d12654c [project @ 1995-09-05 21:51:54 by jwe] jwe parents: 1315 diff changeset	39 // Diagonal Matrix class.
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	40
2385 170053c0f75e [project @ 1996-10-12 18:23:01 by jwe] jwe parents: 1993 diff changeset	41 bool
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	42 DiagMatrix::operator == (const DiagMatrix& a) const
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	43 {
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	44 if (rows () != a.rows () \|\| cols () != a.cols ())
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	45 return 0;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	46
9550 3d6a9aea2aea refactor binary & bool ops in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	47 return mx_inline_equal (length (), data (), a.data ());
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	48 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	49
2385 170053c0f75e [project @ 1996-10-12 18:23:01 by jwe] jwe parents: 1993 diff changeset	50 bool
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	51 DiagMatrix::operator != (const DiagMatrix& a) const
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	52 {
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	53 return !(*this == a);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	54 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	55
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	56 DiagMatrix&
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	57 DiagMatrix::fill (double val)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	58 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	59 for (octave_idx_type i = 0; i < length (); i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	60 elem (i, i) = val;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	61 return *this;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	62 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	63
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	64 DiagMatrix&
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	65 DiagMatrix::fill (double val, octave_idx_type beg, octave_idx_type end)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	66 {
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	67 if (beg < 0 \|\| end >= length () \|\| end < beg)
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	68 (*current_liboctave_error_handler) ("range error for fill");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	69
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	70 for (octave_idx_type i = beg; i <= end; i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	71 elem (i, i) = val;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	72
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	73 return *this;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	74 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	75
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	76 DiagMatrix&
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	77 DiagMatrix::fill (const ColumnVector& a)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	78 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	79 octave_idx_type len = length ();
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	80 if (a.numel () != len)
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	81 (*current_liboctave_error_handler) ("range error for fill");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	82
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	83 for (octave_idx_type i = 0; i < len; i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	84 elem (i, i) = a.elem (i);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	85
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	86 return *this;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	87 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	88
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	89 DiagMatrix&
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	90 DiagMatrix::fill (const RowVector& a)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	91 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	92 octave_idx_type len = length ();
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	93 if (a.numel () != len)
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	94 (*current_liboctave_error_handler) ("range error for fill");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	95
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	96 for (octave_idx_type i = 0; i < len; i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	97 elem (i, i) = a.elem (i);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	98
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	99 return *this;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	100 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	101
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	102 DiagMatrix&
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	103 DiagMatrix::fill (const ColumnVector& a, octave_idx_type beg)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	104 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	105 octave_idx_type a_len = a.numel ();
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	106 if (beg < 0 \|\| beg + a_len >= length ())
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	107 (*current_liboctave_error_handler) ("range error for fill");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	108
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	109 for (octave_idx_type i = 0; i < a_len; i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	110 elem (i+beg, i+beg) = a.elem (i);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	111
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	112 return *this;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	113 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	114
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	115 DiagMatrix&
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	116 DiagMatrix::fill (const RowVector& a, octave_idx_type beg)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	117 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	118 octave_idx_type a_len = a.numel ();
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	119 if (beg < 0 \|\| beg + a_len >= length ())
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	120 (*current_liboctave_error_handler) ("range error for fill");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	121
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	122 for (octave_idx_type i = 0; i < a_len; i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	123 elem (i+beg, i+beg) = a.elem (i);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	124
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	125 return *this;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	126 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	127
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	128 DiagMatrix
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7924 diff changeset	129 DiagMatrix::abs (void) const
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7924 diff changeset	130 {
15448 0a0912a9ab6e Replace deprecated DiagArray2<T>::diag calls with DiagArray2<T>::extract_diag Jordi Gutiérrez Hermoso <jordigh@octave.org> parents: 15271 diff changeset	131 return DiagMatrix (extract_diag ().abs (), rows (), columns ());
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7924 diff changeset	132 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7924 diff changeset	133
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7924 diff changeset	134 DiagMatrix
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	135 real (const ComplexDiagMatrix& a)
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	136 {
15448 0a0912a9ab6e Replace deprecated DiagArray2<T>::diag calls with DiagArray2<T>::extract_diag Jordi Gutiérrez Hermoso <jordigh@octave.org> parents: 15271 diff changeset	137 return DiagMatrix (real (a.extract_diag ()), a.rows (), a.cols ());
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	138 }
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	139
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	140 DiagMatrix
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	141 imag (const ComplexDiagMatrix& a)
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	142 {
15448 0a0912a9ab6e Replace deprecated DiagArray2<T>::diag calls with DiagArray2<T>::extract_diag Jordi Gutiérrez Hermoso <jordigh@octave.org> parents: 15271 diff changeset	143 return DiagMatrix (imag (a.extract_diag ()), a.rows (), a.cols ());
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	144 }
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1192 diff changeset	145
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	146 Matrix
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	147 DiagMatrix::extract (octave_idx_type r1, octave_idx_type c1,
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	148 octave_idx_type r2, octave_idx_type c2) const
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	149 {
17663 7975d75f933c Use std::swap in liboctave instead of temporary variable. Rik <rik@octave.org> parents: 15448 diff changeset	150 if (r1 > r2) { std::swap (r1, r2); }
7975d75f933c Use std::swap in liboctave instead of temporary variable. Rik <rik@octave.org> parents: 15448 diff changeset	151 if (c1 > c2) { std::swap (c1, c2); }
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	152
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	153 octave_idx_type new_r = r2 - r1 + 1;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	154 octave_idx_type new_c = c2 - c1 + 1;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	155
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	156 Matrix result (new_r, new_c);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	157
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	158 for (octave_idx_type j = 0; j < new_c; j++)
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	159 for (octave_idx_type i = 0; i < new_r; i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	160 result.elem (i, j) = elem (r1+i, c1+j);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	161
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	162 return result;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	163 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	164
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	165 // extract row or column i.
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	166
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	167 RowVector
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	168 DiagMatrix::row (octave_idx_type i) const
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	169 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	170 octave_idx_type r = rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	171 octave_idx_type c = cols ();
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	172 if (i < 0 \|\| i >= r)
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	173 (*current_liboctave_error_handler) ("invalid row selection");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	174
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	175 RowVector retval (c, 0.0);
26503 a534d1be2aa0 dDiagMatrix.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	176 if (r <= c \|\| i < c)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	177 retval.elem (i) = elem (i, i);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	178
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	179 return retval;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	180 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	181
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	182 RowVector
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	183 DiagMatrix::row (char *s) const
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	184 {
533 883197c5ad75 [project @ 1994-07-20 20:02:51 by jwe] jwe parents: 458 diff changeset	185 if (! s)
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	186 (*current_liboctave_error_handler) ("invalid row selection");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	187
26503 a534d1be2aa0 dDiagMatrix.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	188 char c = s[0];
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	189 if (c == 'f' \|\| c == 'F')
25688 b2917b7858ba maint: Use Octave convention for spacing of C++ cast statements. Rik <rik@octave.org> parents: 25438 diff changeset	190 return row (static_cast<octave_idx_type> (0));
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	191 else if (c == 'l' \|\| c == 'L')
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	192 return row (rows () - 1);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	193 else
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	194 (*current_liboctave_error_handler) ("invalid row selection");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	195 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	196
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	197 ColumnVector
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	198 DiagMatrix::column (octave_idx_type i) const
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	199 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	200 octave_idx_type r = rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	201 octave_idx_type c = cols ();
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	202 if (i < 0 \|\| i >= c)
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	203 (*current_liboctave_error_handler) ("invalid column selection");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	204
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	205 ColumnVector retval (r, 0.0);
26503 a534d1be2aa0 dDiagMatrix.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	206 if (r >= c \|\| i < r)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	207 retval.elem (i) = elem (i, i);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	208
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	209 return retval;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	210 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	211
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	212 ColumnVector
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	213 DiagMatrix::column (char *s) const
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	214 {
533 883197c5ad75 [project @ 1994-07-20 20:02:51 by jwe] jwe parents: 458 diff changeset	215 if (! s)
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	216 (*current_liboctave_error_handler) ("invalid column selection");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	217
26503 a534d1be2aa0 dDiagMatrix.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	218 char c = s[0];
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	219 if (c == 'f' \|\| c == 'F')
25688 b2917b7858ba maint: Use Octave convention for spacing of C++ cast statements. Rik <rik@octave.org> parents: 25438 diff changeset	220 return column (static_cast<octave_idx_type> (0));
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	221 else if (c == 'l' \|\| c == 'L')
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	222 return column (cols () - 1);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	223 else
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	224 (*current_liboctave_error_handler) ("invalid column selection");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	225 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	226
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	227 DiagMatrix
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	228 DiagMatrix::inverse (void) const
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	229 {
8811 20dfb885f877 int -> octave_idx fixes John W. Eaton <jwe@octave.org> parents: 8375 diff changeset	230 octave_idx_type info;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	231 return inverse (info);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	232 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	233
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	234 DiagMatrix
23433 c9fab0bc983e maint: Use convention 'int& x' for naming references. Rik <rik@octave.org> parents: 23220 diff changeset	235 DiagMatrix::inverse (octave_idx_type& info) const
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	236 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	237 octave_idx_type r = rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	238 octave_idx_type c = cols ();
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	239 if (r != c)
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	240 (*current_liboctave_error_handler) ("inverse requires square matrix");
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	241
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	242 DiagMatrix retval (r, c);
1627 8f798f963805 [project @ 1995-11-13 14:37:05 by jwe] jwe parents: 1368 diff changeset	243
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	244 info = 0;
27093 6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	245 octave_idx_type len = r; // alias for readability
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	246 octave_idx_type z_count = 0; // zeros
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	247 octave_idx_type nz_count = 0; // non-zeros
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	248 for (octave_idx_type i = 0; i < len; i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	249 {
27093 6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	250 if (xelem (i, i) == 0.0)
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	251 {
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	252 z_count++;
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	253 if (nz_count > 0)
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	254 break;
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	255 }
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	256 else
27093 6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	257 {
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	258 nz_count++;
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	259 if (z_count > 0)
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	260 break;
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	261 retval.elem (i, i) = 1.0 / xelem (i, i);
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	262 }
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	263 }
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	264 if (nz_count == 0)
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	265 {
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	266 (*current_liboctave_error_handler)
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	267 ("inverse of the null matrix not defined");
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	268 }
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	269 else if (z_count > 0)
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	270 {
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	271 info = -1;
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	272 element_type *data = retval.fortran_vec ();
6e18f0ce268c Inverse of a sparse/diagonal singular matrix should be a sparse/diagonal matrix of Infs. marco.caliari@univr.it parents: 26503 diff changeset	273 std::fill (data, data + len, octave::numeric_limits<double>::Inf ());
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	274 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	275
1627 8f798f963805 [project @ 1995-11-13 14:37:05 by jwe] jwe parents: 1368 diff changeset	276 return retval;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	277 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	278
8840 c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	279 DiagMatrix
18467 c5a101de2d88 Allow pinv to work on Diagonal Matrices with a tolerance (bug #41546). Rik <rik@octave.org> parents: 18084 diff changeset	280 DiagMatrix::pseudo_inverse (double tol) const
8840 c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	281 {
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	282 octave_idx_type r = rows ();
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	283 octave_idx_type c = cols ();
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	284 octave_idx_type len = length ();
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	285
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	286 DiagMatrix retval (c, r);
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	287
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	288 for (octave_idx_type i = 0; i < len; i++)
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	289 {
18481 16b0cd465ecd Handle special case of 0 for pinv with Diagonal matrices. Rik <rik@octave.org> parents: 18467 diff changeset	290 double val = std::abs (elem (i, i));
18496 80b8873c90ba Fix pinv bug with 0 values introduced in cset 16b0cd465ecd. Rik <rik@octave.org> parents: 18494 diff changeset	291 if (val < tol \|\| val == 0.0)
18467 c5a101de2d88 Allow pinv to work on Diagonal Matrices with a tolerance (bug #41546). Rik <rik@octave.org> parents: 18084 diff changeset	292 retval.elem (i, i) = 0.0;
18496 80b8873c90ba Fix pinv bug with 0 values introduced in cset 16b0cd465ecd. Rik <rik@octave.org> parents: 18494 diff changeset	293 else
18494 1fa5bba16218 [mq]: pinv Rik <rik@octave.org> parents: 18481 diff changeset	294 retval.elem (i, i) = 1.0 / elem (i, i);
8840 c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	295 }
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	296
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	297 return retval;
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	298 }
c690e3772583 support diagonal matrices in pinv Jaroslav Hajek <highegg@gmail.com> parents: 8811 diff changeset	299
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	300 // diagonal matrix by diagonal matrix -> diagonal matrix operations
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	301
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	302 // diagonal matrix by diagonal matrix -> diagonal matrix operations
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	303
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	304 DiagMatrix
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	305 operator * (const DiagMatrix& a, const DiagMatrix& b)
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	306 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	307 octave_idx_type a_nr = a.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	308 octave_idx_type a_nc = a.cols ();
2385 170053c0f75e [project @ 1996-10-12 18:23:01 by jwe] jwe parents: 1993 diff changeset	309
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	310 octave_idx_type b_nr = b.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	311 octave_idx_type b_nc = b.cols ();
2385 170053c0f75e [project @ 1996-10-12 18:23:01 by jwe] jwe parents: 1993 diff changeset	312
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	313 if (a_nc != b_nr)
22327 d0562b3159c7 move more classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	314 octave::err_nonconformant ("operator *", a_nr, a_nc, b_nr, b_nc);
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	315
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	316 DiagMatrix c (a_nr, b_nc);
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	317
18084 8e056300994b Follow coding convention of defining and initializing only 1 variable per line in liboctave. Rik <rik@octave.org> parents: 17769 diff changeset	318 octave_idx_type len = c.length ();
23450 855122b993da maint: Wrap tertiary operator in parentheses "(COND ? x : y)". Rik <rik@octave.org> parents: 23433 diff changeset	319 octave_idx_type lenm = (len < a_nc ? len : a_nc);
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	320
9698 7c6d5d8c8d37 fix diagdiag multiplication Jaroslav Hajek <highegg@gmail.com>* parents: 9601 diff changeset	321 for (octave_idx_type i = 0; i < lenm; i++)
7c6d5d8c8d37 fix diagdiag multiplication Jaroslav Hajek <highegg@gmail.com>* parents: 9601 diff changeset	322 c.dgxelem (i) = a.dgelem (i) * b.dgelem (i);
7c6d5d8c8d37 fix diagdiag multiplication Jaroslav Hajek <highegg@gmail.com>* parents: 9601 diff changeset	323 for (octave_idx_type i = lenm; i < len; i++)
7c6d5d8c8d37 fix diagdiag multiplication Jaroslav Hajek <highegg@gmail.com>* parents: 9601 diff changeset	324 c.dgxelem (i) = 0.0;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	325
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	326 return c;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	327 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	328
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	329 // other operations
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	330
8371 c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	331 DET
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	332 DiagMatrix::determinant (void) const
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	333 {
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	334 DET det (1.0);
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	335 if (rows () != cols ())
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	336 (*current_liboctave_error_handler) ("determinant requires square matrix");
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	337
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	338 octave_idx_type len = length ();
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	339 for (octave_idx_type i = 0; i < len; i++)
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	340 det *= elem (i, i);
8371 c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	341
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	342 return det;
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	343 }
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	344
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	345 double
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	346 DiagMatrix::rcond (void) const
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	347 {
15448 0a0912a9ab6e Replace deprecated DiagArray2<T>::diag calls with DiagArray2<T>::extract_diag Jordi Gutiérrez Hermoso <jordigh@octave.org> parents: 15271 diff changeset	348 ColumnVector av = extract_diag (0).map<double> (fabs);
18084 8e056300994b Follow coding convention of defining and initializing only 1 variable per line in liboctave. Rik <rik@octave.org> parents: 17769 diff changeset	349 double amx = av.max ();
8e056300994b Follow coding convention of defining and initializing only 1 variable per line in liboctave. Rik <rik@octave.org> parents: 17769 diff changeset	350 double amn = av.min ();
8371 c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	351 return amx == 0 ? 0.0 : amn / amx;
c3f7e2549abb make det & inv aware of diagonal & permutation matrices Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	352 }
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	353
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	354 std::ostream&
5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	355 operator << (std::ostream& os, const DiagMatrix& a)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	356 {
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	357 // int field_width = os.precision () + 7;
1360 7eb93d12654c [project @ 1995-09-05 21:51:54 by jwe] jwe parents: 1315 diff changeset	358
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	359 for (octave_idx_type i = 0; i < a.rows (); i++)
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	360 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	361 for (octave_idx_type j = 0; j < a.cols (); j++)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	362 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	363 if (i == j)
23807 336f89b6208b Use character literals 'c' rather than string literals "c" when possible. Rik <rik@octave.org> parents: 23461 diff changeset	364 os << ' ' /* setw (field_width) */ << a.elem (i, i);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	365 else
23807 336f89b6208b Use character literals 'c' rather than string literals "c" when possible. Rik <rik@octave.org> parents: 23461 diff changeset	366 os << ' ' /* setw (field_width) */ << 0.0;
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	367 }
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	368 os << "\n";
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	369 }
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	370 return os;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	371 }

Mercurial > octave

annotate liboctave/array/dDiagMatrix.cc @ 27919:1891570abac8