Mercurial > hg > octave-nkf
annotate liboctave/Array2.h @ 10050:dc88a0b6472c
support old style jacobian for fsolve
| author | Jaroslav Hajek <highegg@gmail.com> |
|---|---|
| date | Sat, 02 Jan 2010 21:55:17 +0100 |
| parents | c0b54271904b |
| children | 4c0cdbe0acca |
| 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 Array2 (T *d, octave_idx_type r, octave_idx_type c) : Array<T> (d, dim_vector (r, c)) { } |
| 1988 | 47 |
| 48 public: | |
| 49 | |
| 4513 | 50 Array2 (void) : Array<T> (dim_vector (0, 0)) { } |
| 1988 | 51 |
| 5275 | 52 Array2 (octave_idx_type r, octave_idx_type c) : Array<T> (dim_vector (r, c)) { } |
| 1988 | 53 |
| 5275 | 54 Array2 (octave_idx_type r, octave_idx_type c, const T& val) |
| 4513 | 55 : Array<T> (dim_vector (r, c), val) { } |
| 1988 | 56 |
| 6990 | 57 Array2 (const dim_vector& dv) : Array<T> (dv) |
| 58 { | |
| 59 if (dv.length () != 2) | |
| 60 (*current_liboctave_error_handler) ("too many dimensions"); | |
| 61 } | |
| 6979 | 62 |
| 6990 | 63 Array2 (const dim_vector& dv, const T& val) : Array<T> (dv) |
| 64 { | |
| 65 if (dv.length () != 2) | |
| 66 (*current_liboctave_error_handler) ("too many dimensions"); | |
| 67 else | |
| 68 Array<T>::fill (val); | |
| 69 } | |
| 6979 | 70 |
| 4513 | 71 Array2 (const Array2<T>& a) : Array<T> (a, a.dims ()) { } |
| 1988 | 72 |
| 5275 | 73 Array2 (const Array<T>& a, octave_idx_type r, octave_idx_type c) |
| 4513 | 74 : Array<T> (a, dim_vector (r, c)) { } |
| 1988 | 75 |
|
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
76 template <class U> |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
77 Array2 (const Array<U>& a) : Array<T> (a) { } |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
78 |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
79 template <class U> |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
80 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
|
81 : Array<T> (a, dv) { } |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
82 |
| 1988 | 83 ~Array2 (void) { } |
| 84 | |
| 85 Array2<T>& operator = (const Array2<T>& a) | |
| 86 { | |
| 3832 | 87 if (this != &a) |
| 4645 | 88 Array<T>::operator = (a); |
| 1988 | 89 |
| 90 return *this; | |
| 91 } | |
| 92 | |
|
8290
7cbe01c21986
improve dense array indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
7789
diff
changeset
|
93 void resize (octave_idx_type r, octave_idx_type c) |
|
7cbe01c21986
improve dense array indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
7789
diff
changeset
|
94 { Array<T>::resize_fill (r, c, Array<T>::resize_fill_value ()); } |
| 3665 | 95 |
| 5275 | 96 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
|
97 { Array<T>::resize_fill (r, c, val); } |
| 2109 | 98 |
| 5275 | 99 Array2<T>& insert (const Array2<T>& a, octave_idx_type r, octave_idx_type c) |
| 2006 | 100 { |
| 4513 | 101 Array<T>::insert (a, r, c); |
| 102 return *this; | |
| 2006 | 103 } |
| 1988 | 104 |
| 4513 | 105 Array2<T> transpose (void) const |
| 106 { | |
| 107 Array<T> tmp = Array<T>::transpose (); | |
| 108 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); | |
| 109 } | |
| 2382 | 110 |
|
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
111 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
|
112 { |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
113 Array<T> tmp = Array<T>::hermitian (fcn); |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
114 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
115 } |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
116 |
|
8579
7e0f36dfefbe
implement octave_value_list using Array
Jaroslav Hajek <highegg@gmail.com>
parents:
8367
diff
changeset
|
117 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
|
118 const T& rfv = Array<T>::resize_fill_value ()) const |
| 4513 | 119 { |
| 120 Array<T> tmp = Array<T>::index (i, resize_ok, rfv); | |
| 121 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); | |
| 122 } | |
| 2382 | 123 |
|
8579
7e0f36dfefbe
implement octave_value_list using Array
Jaroslav Hajek <highegg@gmail.com>
parents:
8367
diff
changeset
|
124 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
|
125 const T& rfv = Array<T>::resize_fill_value ()) const |
| 4513 | 126 { |
| 127 Array<T> tmp = Array<T>::index (i, j, resize_ok, rfv); | |
| 128 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); | |
| 129 } | |
| 7433 | 130 |
|
7463
2467639bd8c0
eliminate UNDEFINED sort mode
John W. Eaton <jwe@octave.org>
parents:
7433
diff
changeset
|
131 Array2<T> sort (octave_idx_type dim = 0, sortmode mode = ASCENDING) const |
| 7433 | 132 { |
| 133 Array<T> tmp = Array<T>::sort (dim, mode); | |
| 134 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); | |
| 135 } | |
| 136 | |
| 137 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
|
138 sortmode mode = ASCENDING) const |
| 7433 | 139 { |
| 140 Array<T> tmp = Array<T>::sort (sidx, dim, mode); | |
| 141 return Array2<T> (tmp, tmp.rows (), tmp.columns ()); | |
| 142 } | |
|
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
143 |
|
7620
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7503
diff
changeset
|
144 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
|
145 { |
|
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7503
diff
changeset
|
146 return Array<T>::diag (k); |
|
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 |
|
9812
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
149 // FIXME: should go away. |
|
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
150 template <class U> |
|
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
151 Array2<U> |
|
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
152 map (U (&fcn) (T)) const |
|
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
153 { return Array<T>::template map<U> (fcn); } |
|
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
154 |
|
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
155 template <class U> |
|
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
156 Array2<U> |
|
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
157 map (U (&fcn) (const T&)) const |
|
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
158 { return Array<T>::template map<U> (fcn); } |
| 1988 | 159 }; |
| 160 | |
| 161 #endif | |
| 162 | |
| 163 /* | |
| 164 ;;; Local Variables: *** | |
| 165 ;;; mode: C++ *** | |
| 166 ;;; End: *** | |
| 167 */ |