Mercurial > hg > octave-lyh
annotate liboctave/Array2.h @ 8920:eb63fbe60fab
update copyright notices
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Sat, 07 Mar 2009 10:41:27 -0500 |
| parents | 7e0f36dfefbe |
| children | f80c566bc751 |
| rev | line source |
|---|---|
| 1993 | 1 // Template array classes |
| 1988 | 2 /* |
| 3 | |
| 7017 | 4 Copyright (C) 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, |
| 8920 | 5 2005, 2007, 2008, 2009 John W. Eaton |
| 1988 | 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. | |
| 1988 | 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/>. | |
| 1988 | 22 |
| 23 */ | |
| 24 | |
| 25 #if !defined (octave_Array2_h) | |
| 26 #define octave_Array2_h 1 | |
| 27 | |
| 28 #include <cassert> | |
| 3473 | 29 #include <climits> |
| 1988 | 30 #include <cstdlib> |
| 31 | |
| 32 #include "Array.h" | |
| 33 #include "lo-error.h" | |
| 7231 | 34 #include "lo-math.h" |
| 1988 | 35 |
| 36 class idx_vector; | |
| 37 | |
| 38 // Two dimensional array class. | |
| 39 | |
| 40 template <class T> | |
| 3585 | 41 class |
| 42 Array2 : public Array<T> | |
| 1988 | 43 { |
| 3492 | 44 protected: |
| 3473 | 45 |
| 5275 | 46 static octave_idx_type get_size (octave_idx_type r, octave_idx_type c) { return Array<T>::get_size (r, c); } |
| 3473 | 47 |
| 5275 | 48 Array2 (T *d, octave_idx_type r, octave_idx_type c) : Array<T> (d, dim_vector (r, c)) { } |
| 1988 | 49 |
| 50 public: | |
| 51 | |
| 4513 | 52 Array2 (void) : Array<T> (dim_vector (0, 0)) { } |
| 1988 | 53 |
| 5275 | 54 Array2 (octave_idx_type r, octave_idx_type c) : Array<T> (dim_vector (r, c)) { } |
| 1988 | 55 |
| 5275 | 56 Array2 (octave_idx_type r, octave_idx_type c, const T& val) |
| 4513 | 57 : Array<T> (dim_vector (r, c), val) { } |
| 1988 | 58 |
| 6990 | 59 Array2 (const dim_vector& dv) : Array<T> (dv) |
| 60 { | |
| 61 if (dv.length () != 2) | |
| 62 (*current_liboctave_error_handler) ("too many dimensions"); | |
| 63 } | |
| 6979 | 64 |
| 6990 | 65 Array2 (const dim_vector& dv, const T& val) : Array<T> (dv) |
| 66 { | |
| 67 if (dv.length () != 2) | |
| 68 (*current_liboctave_error_handler) ("too many dimensions"); | |
| 69 else | |
| 70 Array<T>::fill (val); | |
| 71 } | |
| 6979 | 72 |
| 4513 | 73 Array2 (const Array2<T>& a) : Array<T> (a, a.dims ()) { } |
| 1988 | 74 |
| 5275 | 75 Array2 (const Array<T>& a, octave_idx_type r, octave_idx_type c) |
| 4513 | 76 : Array<T> (a, dim_vector (r, c)) { } |
| 1988 | 77 |
|
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
78 template <class U> |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
79 Array2 (const Array<U>& a) : Array<T> (a) { } |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
80 |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
81 template <class U> |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
82 Array2 (const Array<U>& a, const dim_vector& dv) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
83 : Array<T> (a, dv) { } |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
84 |
| 1988 | 85 ~Array2 (void) { } |
| 86 | |
| 87 Array2<T>& operator = (const Array2<T>& a) | |
| 88 { | |
| 3832 | 89 if (this != &a) |
| 4645 | 90 Array<T>::operator = (a); |
| 1988 | 91 |
| 92 return *this; | |
| 93 } | |
| 94 | |
|
8290
7cbe01c21986
improve dense array indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
7789
diff
changeset
|
95 void resize (octave_idx_type r, octave_idx_type c) |
|
7cbe01c21986
improve dense array indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
7789
diff
changeset
|
96 { Array<T>::resize_fill (r, c, Array<T>::resize_fill_value ()); } |
| 3665 | 97 |
| 5275 | 98 void resize (octave_idx_type r, octave_idx_type c, const T& val) |
|
8290
7cbe01c21986
improve dense array indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
7789
diff
changeset
|
99 { Array<T>::resize_fill (r, c, val); } |
| 2109 | 100 |
| 5275 | 101 Array2<T>& insert (const Array2<T>& a, octave_idx_type r, octave_idx_type c) |
| 2006 | 102 { |
| 4513 | 103 Array<T>::insert (a, r, c); |
| 104 return *this; | |
| 2006 | 105 } |
| 1988 | 106 |
| 4513 | 107 Array2<T> transpose (void) const |
| 108 { | |
| 109 Array<T> tmp = Array<T>::transpose (); | |
| 110 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); | |
| 111 } | |
| 2382 | 112 |
|
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
113 Array2<T> hermitian (T (*fcn) (const T&) = 0) const |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
114 { |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
115 Array<T> tmp = Array<T>::hermitian (fcn); |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
116 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
117 } |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
118 |
|
8579
7e0f36dfefbe
implement octave_value_list using Array
Jaroslav Hajek <highegg@gmail.com>
parents:
8367
diff
changeset
|
119 Array2<T> index (const idx_vector& i, bool resize_ok = false, |
|
8323
655193554ad2
Correct use of resize_fill_value in Array2<T>::index methods.
David Bateman <dbateman@free.fr>
parents:
8290
diff
changeset
|
120 const T& rfv = Array<T>::resize_fill_value ()) const |
| 4513 | 121 { |
| 122 Array<T> tmp = Array<T>::index (i, resize_ok, rfv); | |
| 123 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); | |
| 124 } | |
| 2382 | 125 |
|
8579
7e0f36dfefbe
implement octave_value_list using Array
Jaroslav Hajek <highegg@gmail.com>
parents:
8367
diff
changeset
|
126 Array2<T> index (const idx_vector& i, const idx_vector& j, bool resize_ok = false, |
|
8323
655193554ad2
Correct use of resize_fill_value in Array2<T>::index methods.
David Bateman <dbateman@free.fr>
parents:
8290
diff
changeset
|
127 const T& rfv = Array<T>::resize_fill_value ()) const |
| 4513 | 128 { |
| 129 Array<T> tmp = Array<T>::index (i, j, resize_ok, rfv); | |
| 130 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); | |
| 131 } | |
| 7433 | 132 |
|
7463
2467639bd8c0
eliminate UNDEFINED sort mode
John W. Eaton <jwe@octave.org>
parents:
7433
diff
changeset
|
133 Array2<T> sort (octave_idx_type dim = 0, sortmode mode = ASCENDING) const |
| 7433 | 134 { |
| 135 Array<T> tmp = Array<T>::sort (dim, mode); | |
| 136 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); | |
| 137 } | |
| 138 | |
| 139 Array2<T> sort (Array<octave_idx_type> &sidx, octave_idx_type dim = 0, | |
|
7463
2467639bd8c0
eliminate UNDEFINED sort mode
John W. Eaton <jwe@octave.org>
parents:
7433
diff
changeset
|
140 sortmode mode = ASCENDING) const |
| 7433 | 141 { |
| 142 Array<T> tmp = Array<T>::sort (sidx, dim, mode); | |
| 143 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); | |
| 144 } | |
|
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
145 |
|
7620
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7503
diff
changeset
|
146 Array2<T> diag (octave_idx_type k) const |
|
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7503
diff
changeset
|
147 { |
|
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7503
diff
changeset
|
148 return Array<T>::diag (k); |
|
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7503
diff
changeset
|
149 } |
|
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7503
diff
changeset
|
150 |
|
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
151 template <class U, class F> |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
152 Array2<U> map (F fcn) const |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
153 { |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
154 return Array<T>::template map<U> (fcn); |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
155 } |
| 1988 | 156 }; |
| 157 | |
| 158 #endif | |
| 159 | |
| 160 /* | |
| 161 ;;; Local Variables: *** | |
| 162 ;;; mode: C++ *** | |
| 163 ;;; End: *** | |
| 164 */ |