octave-dspies: liboctave/DiagArray2.h annotate

annotate liboctave/DiagArray2.h @ 10312:cbc402e64d83

untabify liboctave header files

author	John W. Eaton <jwe@octave.org>
date	Thu, 11 Feb 2010 12:14:48 -0500
parents	4c0cdbe0acca
children	12884915a8e4

rev	line source
1993 1b57120c997b [project @ 1996-03-03 01:16:15 by jwe] jwe parents: 1988 diff changeset	1 // Template array classes
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	2 /*
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	3
7017 a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe] jwe parents: 7016 diff changeset	4 Copyright (C) 1996, 1997, 2000, 2002, 2003, 2004, 2005, 2006, 2007
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe] jwe parents: 7016 diff changeset	5 John W. Eaton
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	6 Copyright (C) 2008, 2009 Jaroslav Hajek
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	7
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	8 This file is part of Octave.
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	9
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	10 Octave is free software; you can redistribute it and/or modify it
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	11 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5775 diff changeset	12 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5775 diff changeset	13 option) any later version.
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	14
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	15 Octave is distributed in the hope that it will be useful, but WITHOUT
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	16 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	17 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	18 for more details.
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	19
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	20 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: 5775 diff changeset	21 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5775 diff changeset	22 <http://www.gnu.org/licenses/>.
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	23
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	24 */
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	25
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	26 #if !defined (octave_DiagArray2_h)
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	27 #define octave_DiagArray2_h 1
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	28
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	29 #include <cassert>
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	30 #include <cstdlib>
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	31
4645 bd2067547b40 [project @ 2003-11-23 08:07:52 by jwe] jwe parents: 4513 diff changeset	32 #include "Array.h"
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	33 #include "Array2.h"
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	34 #include "lo-error.h"
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	35
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	36 // A two-dimensional array with diagonal elements only.
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	37 // Idea and example code for Proxy class and functions from:
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	38 //
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	39 // From: kanze@us-es.sel.de (James Kanze)
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	40 // Subject: Re: How to overload [] to do READ/WRITE differently ?
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	41 // Message-ID: <KANZE.93Nov29151407@slsvhdt.us-es.sel.de>
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	42 // Sender: news@us-es.sel.de
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	43 // Date: 29 Nov 1993 14:14:07 GMT
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	44 // --
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	45 // James Kanze email: kanze@us-es.sel.de
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	46 // GABI Software, Sarl., 8 rue du Faisan, F-67000 Strasbourg, France
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	47
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	48 // Array<T> is inherited privately so that some methods, like index, don't
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	49 // produce unexpected results.
8375 e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	50
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	51 template <class T>
3585 d9803711e047 [project @ 2000-02-08 04:35:39 by jwe] jwe parents: 2847 diff changeset	52 class
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	53 DiagArray2 : protected Array<T>
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	54 {
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	55 private:
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	56
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	57 T get (octave_idx_type i) { return Array<T>::xelem (i); }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	58
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	59 void set (const T& val, octave_idx_type i) { Array<T>::xelem (i) = val; }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	60
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	61 class Proxy
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	62 {
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	63 public:
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	64
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	65 Proxy (DiagArray2<T> *ref, octave_idx_type r, octave_idx_type c)
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	66 : i (r), j (c), object (ref) { }
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	67
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	68 const Proxy& operator = (const T& val) const;
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	69
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	70 operator T () const;
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	71
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	72 private:
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	73
5775 ace8d8d26933 [project @ 2006-04-24 19:13:06 by jwe] jwe parents: 5760 diff changeset	74 // FIXME -- this is declared private to keep the user from
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	75 // taking the address of a Proxy. Maybe it should be implemented
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	76 // by means of a companion function in the DiagArray2 class.
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	77
5760 8d7162924bd3 [project @ 2006-04-14 04:01:37 by jwe] jwe parents: 5752 diff changeset	78 T *operator& () const { assert (0); return 0; }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	79
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	80 octave_idx_type i;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	81 octave_idx_type j;
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	82
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	83 DiagArray2<T> *object;
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	84
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	85 };
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	86
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	87 friend class Proxy;
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	88
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	89 protected:
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	90 octave_idx_type d1, d2;
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	91
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	92 DiagArray2 (T *d, octave_idx_type r, octave_idx_type c)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	93 : Array<T> (d, std::min (r, c)), d1 (r), d2 (c) { }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	94
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	95 public:
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	96
8918 f5408862892f Consistently use element_type in the array classes. Jason Riedy <jason@acm.org> parents: 8524 diff changeset	97 using Array<T>::element_type;
8375 e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	98
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	99 DiagArray2 (void)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	100 : Array<T> (), d1 (0), d2 (0) { }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	101
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	102 DiagArray2 (octave_idx_type r, octave_idx_type c)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	103 : Array<T> (std::min (r, c)), d1 (r), d2 (c) { }
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4192 diff changeset	104
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	105 DiagArray2 (octave_idx_type r, octave_idx_type c, const T& val)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	106 : Array<T> (std::min (r, c), val), d1 (r), d2 (c) { }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	107
9612 66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 8918 diff changeset	108 DiagArray2 (const dim_vector& dv)
66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 8918 diff changeset	109 : Array<T> (std::min (dv(0), dv(1))), d1 (dv(0)), d2 (dv(0))
66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 8918 diff changeset	110 {
66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 8918 diff changeset	111 if (dv.length () != 2)
10312 cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	112 (*current_liboctave_error_handler) ("too many dimensions");
9612 66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 8918 diff changeset	113 }
66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 8918 diff changeset	114
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	115 DiagArray2 (const Array<T>& a)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	116 : Array<T> (a), d1 (a.numel ()), d2 (a.numel ()) { }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	117
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	118 DiagArray2 (const DiagArray2<T>& a)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	119 : Array<T> (a), d1 (a.d1), d2 (a.d2) { }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	120
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	121 template <class U>
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	122 DiagArray2 (const DiagArray2<U>& a)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	123 : Array<T> (a.diag ()), d1 (a.dim1 ()), d2 (a.dim2 ()) { }
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	124
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	125 ~DiagArray2 (void) { }
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	126
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	127 DiagArray2<T>& operator = (const DiagArray2<T>& a)
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	128 {
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	129 if (this != &a)
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	130 {
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	131 Array<T>::operator = (a);
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	132 d1 = a.d1;
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	133 d2 = a.d2;
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	134 }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	135
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	136 return *this;
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	137 }
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	138
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	139 octave_idx_type dim1 (void) const { return d1; }
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	140 octave_idx_type dim2 (void) const { return d2; }
8375 e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	141
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	142 octave_idx_type rows (void) const { return dim1 (); }
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	143 octave_idx_type cols (void) const { return dim2 (); }
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	144 octave_idx_type columns (void) const { return dim2 (); }
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	145
9698 7c6d5d8c8d37 fix diagdiag multiplication Jaroslav Hajek <highegg@gmail.com>* parents: 9612 diff changeset	146 octave_idx_type diag_length (void) const { return Array<T>::length (); }
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	147 // FIXME: a dangerous ambiguity?
8375 e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	148 octave_idx_type length (void) const { return Array<T>::length (); }
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	149 octave_idx_type nelem (void) const { return dim1 () * dim2 (); }
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	150 octave_idx_type numel (void) const { return nelem (); }
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	151
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	152 size_t byte_size (void) const { return length () * sizeof (T); }
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	153
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	154 dim_vector dims (void) const { return dim_vector (d1, d2); }
8375 e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	155
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	156 Array<T> diag (octave_idx_type k = 0) const;
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	157
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	158 // Warning: the non-const two-index versions will silently ignore assignments
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	159 // to off-diagonal elements.
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	160
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	161 T elem (octave_idx_type r, octave_idx_type c) const
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	162 {
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	163 return (r == c) ? Array<T>::elem (r) : T (0);
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	164 }
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	165
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	166 T& elem (octave_idx_type r, octave_idx_type c)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	167 {
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	168 static T zero (0);
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	169 return (r == c) ? Array<T>::elem (r) : zero;
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	170 }
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	171
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	172 T dgelem (octave_idx_type i) const
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	173 { return Array<T>::elem (i); }
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	174
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	175 T& dgelem (octave_idx_type i)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	176 { return Array<T>::elem (i); }
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	177
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	178 T checkelem (octave_idx_type r, octave_idx_type c) const;
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	179 Proxy checkelem (octave_idx_type r, octave_idx_type c);
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	180
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	181 T operator () (octave_idx_type r, octave_idx_type c) const
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	182 {
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	183 #if defined (BOUNDS_CHECKING)
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	184 return checkelem (r, c);
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	185 #else
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	186 return elem (r, c);
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	187 #endif
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	188 }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	189
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	190 // FIXME: can this cause problems?
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	191 #if defined (BOUNDS_CHECKING)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	192 Proxy operator () (octave_idx_type r, octave_idx_type c)
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	193 {
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	194 return checkelem (r, c);
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	195 }
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	196 #else
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	197 T& operator () (octave_idx_type r, octave_idx_type c)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	198 {
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	199 return elem (r, c);
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	200 }
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	201 #endif
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	202
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	203 // No checking.
8366 8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	204
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	205 T xelem (octave_idx_type r, octave_idx_type c) const
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	206 {
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	207 return (r == c) ? Array<T>::xelem (r) : T (0);
8b1a2555c4e2 implement diagonal matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8314 diff changeset	208 }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	209
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	210 T& dgxelem (octave_idx_type i)
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	211 { return Array<T>::xelem (i); }
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	212
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	213 T dgxelem (octave_idx_type i) const
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	214 { return Array<T>::xelem (i); }
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	215
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	216 void resize (octave_idx_type n, octave_idx_type m);
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	217 void resize_fill (octave_idx_type n, octave_idx_type m, const T& val);
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	218
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	219 DiagArray2<T> transpose (void) const;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	220 DiagArray2<T> hermitian (T (*fcn) (const T&) = 0) const;
8375 e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	221
8524 937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	222 operator Array2<T> (void) const;
937921654627 clean up Array and DiagArray2 Jaroslav Hajek <highegg@gmail.com> parents: 8375 diff changeset	223
8375 e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	224 const T *data (void) const { return Array<T>::data (); }
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	225
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	226 const T *fortran_vec (void) const { return Array<T>::fortran_vec (); }
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	227
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	228 T *fortran_vec (void) { return Array<T>::fortran_vec (); }
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	229
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	230 void print_info (std::ostream& os, const std::string& prefix) const
e3c9102431a9 fix design problems of diag & perm matrix classes Jaroslav Hajek <highegg@gmail.com> parents: 8366 diff changeset	231 { Array<T>::print_info (os, prefix); }
1988 7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	232 };
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	233
7b56630a1e05 [project @ 1996-03-02 00:33:22 by jwe] jwe parents: diff changeset	234 #endif

Mercurial > octave-dspies

annotate liboctave/DiagArray2.h @ 10312:cbc402e64d83