Mercurial > hg > octave-nkf
annotate liboctave/CNDArray.h @ 9743:26abff55f6fe
optimize bsxfun for common built-in operations
| author | Jaroslav Hajek <highegg@gmail.com> |
|---|---|
| date | Tue, 20 Oct 2009 10:47:22 +0200 |
| parents | b4fdfee405b5 |
| children | a5035bc7fbfb |
| rev | line source |
|---|---|
| 4514 | 1 /* |
| 2 | |
| 8920 | 3 Copyright (C) 2003, 2004, 2005, 2006, 2007, 2008, 2009 John W. Eaton |
| 4514 | 4 |
| 5 This file is part of Octave. | |
| 6 | |
| 7 Octave is free software; you can redistribute it and/or modify it | |
| 8 under the terms of the GNU General Public License as published by the | |
| 7016 | 9 Free Software Foundation; either version 3 of the License, or (at your |
| 10 option) any later version. | |
| 4514 | 11 |
| 12 Octave is distributed in the hope that it will be useful, but WITHOUT | |
| 13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
| 14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
| 15 for more details. | |
| 16 | |
| 17 You should have received a copy of the GNU General Public License | |
| 7016 | 18 along with Octave; see the file COPYING. If not, see |
| 19 <http://www.gnu.org/licenses/>. | |
| 4514 | 20 |
| 21 */ | |
| 22 | |
| 23 #if !defined (octave_ComplexNDArray_h) | |
| 24 #define octave_ComplexNDArray_h 1 | |
| 25 | |
| 26 #include "MArrayN.h" | |
| 27 #include "CMatrix.h" | |
| 28 | |
| 29 #include "mx-defs.h" | |
|
8774
b756ce0002db
split implementation and interface in mx-op-defs and MArray-defs
Jaroslav Hajek <highegg@gmail.com>
parents:
8650
diff
changeset
|
30 #include "mx-op-decl.h" |
|
9743
26abff55f6fe
optimize bsxfun for common built-in operations
Jaroslav Hajek <highegg@gmail.com>
parents:
9732
diff
changeset
|
31 #include "bsxfun-decl.h" |
| 4514 | 32 |
| 33 class | |
| 6108 | 34 OCTAVE_API |
| 4514 | 35 ComplexNDArray : public MArrayN<Complex> |
| 36 { | |
| 37 public: | |
| 5523 | 38 |
|
9656
b29504415a2e
provide NDArray->Matrix->Vector typedef mappers
Jaroslav Hajek <highegg@gmail.com>
parents:
9607
diff
changeset
|
39 typedef ComplexMatrix matrix_type; |
|
b29504415a2e
provide NDArray->Matrix->Vector typedef mappers
Jaroslav Hajek <highegg@gmail.com>
parents:
9607
diff
changeset
|
40 |
| 4514 | 41 ComplexNDArray (void) : MArrayN<Complex> () { } |
| 42 | |
| 4587 | 43 ComplexNDArray (const dim_vector& dv) : MArrayN<Complex> (dv) { } |
| 4514 | 44 |
| 4587 | 45 ComplexNDArray (const dim_vector& dv, const Complex& val) |
| 46 : MArrayN<Complex> (dv, val) { } | |
| 4514 | 47 |
| 48 ComplexNDArray (const ComplexNDArray& a) : MArrayN<Complex> (a) { } | |
| 49 | |
| 50 ComplexNDArray (const ComplexMatrix& a) : MArrayN<Complex> (a) { } | |
| 51 | |
|
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
52 template <class U> |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
53 ComplexNDArray (const MArrayN<U>& a) : MArrayN<Complex> (a) { } |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
54 |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
55 template <class U> |
|
9732
b4fdfee405b5
remove ArrayN<T> + fix nonhom. diag-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents:
9721
diff
changeset
|
56 ComplexNDArray (const Array<U>& a) : MArrayN<Complex> (a) { } |
| 4514 | 57 |
|
8956
d91fa4b20bbb
ensure nonnegative char -> real conversion
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
58 ComplexNDArray (const charNDArray&); |
|
d91fa4b20bbb
ensure nonnegative char -> real conversion
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
59 |
| 4514 | 60 ComplexNDArray& operator = (const ComplexNDArray& a) |
| 61 { | |
| 62 MArrayN<Complex>::operator = (a); | |
| 63 return *this; | |
| 64 } | |
| 65 | |
| 4543 | 66 // unary operations |
| 67 | |
| 68 boolNDArray operator ! (void) const; | |
| 69 | |
| 5775 | 70 // FIXME -- this is not quite the right thing. |
| 4514 | 71 |
|
7922
935be827eaf8
error for NaN values in & and | expressions
John W. Eaton <jwe@octave.org>
parents:
7789
diff
changeset
|
72 bool any_element_is_nan (void) const; |
| 4687 | 73 bool any_element_is_inf_or_nan (void) const; |
| 74 bool all_elements_are_real (void) const; | |
| 75 bool all_integers (double& max_val, double& min_val) const; | |
| 76 bool too_large_for_float (void) const; | |
| 77 | |
| 4556 | 78 boolNDArray all (int dim = -1) const; |
| 79 boolNDArray any (int dim = -1) const; | |
| 4514 | 80 |
| 4584 | 81 ComplexNDArray cumprod (int dim = -1) const; |
| 82 ComplexNDArray cumsum (int dim = -1) const; | |
| 4569 | 83 ComplexNDArray prod (int dim = -1) const; |
| 84 ComplexNDArray sum (int dim = -1) const; | |
|
9721
192d94cff6c1
improve sum & implement the 'extra' option, refactor some code
Jaroslav Hajek <highegg@gmail.com>
parents:
9656
diff
changeset
|
85 ComplexNDArray xsum (int dim = -1) const; |
| 4569 | 86 ComplexNDArray sumsq (int dim = -1) const; |
| 5275 | 87 ComplexNDArray concat (const ComplexNDArray& rb, const Array<octave_idx_type>& ra_idx); |
| 88 ComplexNDArray concat (const NDArray& rb, const Array<octave_idx_type>& ra_idx); | |
| 4569 | 89 |
| 4844 | 90 ComplexNDArray max (int dim = 0) const; |
|
9732
b4fdfee405b5
remove ArrayN<T> + fix nonhom. diag-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents:
9721
diff
changeset
|
91 ComplexNDArray max (Array<octave_idx_type>& index, int dim = 0) const; |
| 4844 | 92 ComplexNDArray min (int dim = 0) const; |
|
9732
b4fdfee405b5
remove ArrayN<T> + fix nonhom. diag-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents:
9721
diff
changeset
|
93 ComplexNDArray min (Array<octave_idx_type>& index, int dim = 0) const; |
|
8777
724c0f46d9d4
implement cummin/cummax functions
Jaroslav Hajek <highegg@gmail.com>
parents:
8774
diff
changeset
|
94 |
|
724c0f46d9d4
implement cummin/cummax functions
Jaroslav Hajek <highegg@gmail.com>
parents:
8774
diff
changeset
|
95 ComplexNDArray cummax (int dim = 0) const; |
|
9732
b4fdfee405b5
remove ArrayN<T> + fix nonhom. diag-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents:
9721
diff
changeset
|
96 ComplexNDArray cummax (Array<octave_idx_type>& index, int dim = 0) const; |
|
8777
724c0f46d9d4
implement cummin/cummax functions
Jaroslav Hajek <highegg@gmail.com>
parents:
8774
diff
changeset
|
97 ComplexNDArray cummin (int dim = 0) const; |
|
9732
b4fdfee405b5
remove ArrayN<T> + fix nonhom. diag-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents:
9721
diff
changeset
|
98 ComplexNDArray cummin (Array<octave_idx_type>& index, int dim = 0) const; |
|
8777
724c0f46d9d4
implement cummin/cummax functions
Jaroslav Hajek <highegg@gmail.com>
parents:
8774
diff
changeset
|
99 |
|
9513
9f870f73ab7d
implement built-in diff
Jaroslav Hajek <highegg@gmail.com>
parents:
9237
diff
changeset
|
100 ComplexNDArray diff (octave_idx_type order = 1, int dim = 0) const; |
|
9f870f73ab7d
implement built-in diff
Jaroslav Hajek <highegg@gmail.com>
parents:
9237
diff
changeset
|
101 |
| 5275 | 102 ComplexNDArray& insert (const NDArray& a, octave_idx_type r, octave_idx_type c); |
| 103 ComplexNDArray& insert (const ComplexNDArray& a, octave_idx_type r, octave_idx_type c); | |
| 104 ComplexNDArray& insert (const ComplexNDArray& a, const Array<octave_idx_type>& ra_idx); | |
| 4765 | 105 |
| 4634 | 106 NDArray abs (void) const; |
|
8998
a48fba01e4ac
optimize isnan/isinf/isfinite mappers
Jaroslav Hajek <highegg@gmail.com>
parents:
8956
diff
changeset
|
107 boolNDArray isnan (void) const; |
|
a48fba01e4ac
optimize isnan/isinf/isfinite mappers
Jaroslav Hajek <highegg@gmail.com>
parents:
8956
diff
changeset
|
108 boolNDArray isinf (void) const; |
|
a48fba01e4ac
optimize isnan/isinf/isfinite mappers
Jaroslav Hajek <highegg@gmail.com>
parents:
8956
diff
changeset
|
109 boolNDArray isfinite (void) const; |
| 4569 | 110 |
|
9237
3c1762c7e787
Add missing xxx_API decoration and remove misplaced ones
Michael Goffioul <michael.goffioul@gmail.com>
parents:
8998
diff
changeset
|
111 friend OCTAVE_API ComplexNDArray conj (const ComplexNDArray& a); |
|
8650
a1ae2aae903e
abs,real,imag,conj: use code from mx-inlines rather than the generic map
Jaroslav Hajek <highegg@gmail.com>
parents:
7922
diff
changeset
|
112 |
| 4773 | 113 ComplexNDArray fourier (int dim = 1) const; |
| 114 ComplexNDArray ifourier (int dim = 1) const; | |
| 115 | |
| 116 ComplexNDArray fourier2d (void) const; | |
| 117 ComplexNDArray ifourier2d (void) const; | |
| 118 | |
| 119 ComplexNDArray fourierNd (void) const; | |
| 120 ComplexNDArray ifourierNd (void) const; | |
| 121 | |
| 4514 | 122 ComplexMatrix matrix_value (void) const; |
| 123 | |
| 4902 | 124 ComplexNDArray squeeze (void) const { return MArrayN<Complex>::squeeze (); } |
| 4532 | 125 |
| 5275 | 126 static void increment_index (Array<octave_idx_type>& ra_idx, |
| 4532 | 127 const dim_vector& dimensions, |
| 128 int start_dimension = 0); | |
| 129 | |
| 5275 | 130 static octave_idx_type compute_index (Array<octave_idx_type>& ra_idx, |
| 4556 | 131 const dim_vector& dimensions); |
| 132 | |
| 4514 | 133 // i/o |
| 134 | |
| 6108 | 135 friend OCTAVE_API std::ostream& operator << (std::ostream& os, const ComplexNDArray& a); |
| 136 friend OCTAVE_API std::istream& operator >> (std::istream& is, ComplexNDArray& a); | |
| 4514 | 137 |
| 138 static Complex resize_fill_value (void) { return Complex (0.0, 0.0); } | |
| 139 | |
| 140 // bool all_elements_are_real (void) const; | |
| 141 // bool all_integers (double& max_val, double& min_val) const; | |
| 142 | |
|
7620
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7503
diff
changeset
|
143 ComplexNDArray diag (octave_idx_type k = 0) const; |
|
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7503
diff
changeset
|
144 |
|
9607
1be3c73ed7b5
reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents:
9546
diff
changeset
|
145 ComplexNDArray& changesign (void) |
|
1be3c73ed7b5
reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents:
9546
diff
changeset
|
146 { |
|
1be3c73ed7b5
reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents:
9546
diff
changeset
|
147 MArrayN<Complex>::changesign (); |
|
1be3c73ed7b5
reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents:
9546
diff
changeset
|
148 return *this; |
|
1be3c73ed7b5
reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents:
9546
diff
changeset
|
149 } |
|
1be3c73ed7b5
reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents:
9546
diff
changeset
|
150 |
|
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
151 typedef double (*dmapper) (const Complex&); |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
152 typedef Complex (*cmapper) (const Complex&); |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
153 typedef bool (*bmapper) (const Complex&); |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
154 |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
155 NDArray map (dmapper fcn) const; |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
156 ComplexNDArray map (cmapper fcn) const; |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
157 boolNDArray map (bmapper fcn) const; |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
158 |
| 4514 | 159 private: |
| 160 | |
| 4634 | 161 ComplexNDArray (Complex *d, const dim_vector& dv) |
| 4587 | 162 : MArrayN<Complex> (d, dv) { } |
| 4514 | 163 }; |
| 164 | |
|
8650
a1ae2aae903e
abs,real,imag,conj: use code from mx-inlines rather than the generic map
Jaroslav Hajek <highegg@gmail.com>
parents:
7922
diff
changeset
|
165 extern OCTAVE_API ComplexNDArray conj (const ComplexNDArray& a); |
|
a1ae2aae903e
abs,real,imag,conj: use code from mx-inlines rather than the generic map
Jaroslav Hajek <highegg@gmail.com>
parents:
7922
diff
changeset
|
166 |
| 6108 | 167 extern OCTAVE_API ComplexNDArray min (const Complex& c, const ComplexNDArray& m); |
| 168 extern OCTAVE_API ComplexNDArray min (const ComplexNDArray& m, const Complex& c); | |
| 169 extern OCTAVE_API ComplexNDArray min (const ComplexNDArray& a, const ComplexNDArray& b); | |
| 4844 | 170 |
| 6108 | 171 extern OCTAVE_API ComplexNDArray max (const Complex& c, const ComplexNDArray& m); |
| 172 extern OCTAVE_API ComplexNDArray max (const ComplexNDArray& m, const Complex& c); | |
| 173 extern OCTAVE_API ComplexNDArray max (const ComplexNDArray& a, const ComplexNDArray& b); | |
| 4844 | 174 |
| 6708 | 175 NDS_CMP_OP_DECLS (ComplexNDArray, Complex, OCTAVE_API) |
| 176 NDS_BOOL_OP_DECLS (ComplexNDArray, Complex, OCTAVE_API) | |
| 4543 | 177 |
| 6708 | 178 SND_CMP_OP_DECLS (Complex, ComplexNDArray, OCTAVE_API) |
| 179 SND_BOOL_OP_DECLS (Complex, ComplexNDArray, OCTAVE_API) | |
| 4543 | 180 |
| 6708 | 181 NDND_CMP_OP_DECLS (ComplexNDArray, ComplexNDArray, OCTAVE_API) |
| 182 NDND_BOOL_OP_DECLS (ComplexNDArray, ComplexNDArray, OCTAVE_API) | |
| 4543 | 183 |
| 4514 | 184 MARRAY_FORWARD_DEFS (MArrayN, ComplexNDArray, Complex) |
| 185 | |
|
9546
1beb23d2b892
optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents:
9513
diff
changeset
|
186 extern OCTAVE_API ComplexNDArray& operator *= (ComplexNDArray& a, double s); |
|
1beb23d2b892
optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents:
9513
diff
changeset
|
187 extern OCTAVE_API ComplexNDArray& operator /= (ComplexNDArray& a, double s); |
|
1beb23d2b892
optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents:
9513
diff
changeset
|
188 |
|
9743
26abff55f6fe
optimize bsxfun for common built-in operations
Jaroslav Hajek <highegg@gmail.com>
parents:
9732
diff
changeset
|
189 BSXFUN_STDOP_DECLS (ComplexNDArray, OCTAVE_API) |
|
26abff55f6fe
optimize bsxfun for common built-in operations
Jaroslav Hajek <highegg@gmail.com>
parents:
9732
diff
changeset
|
190 BSXFUN_STDREL_DECLS (ComplexNDArray, OCTAVE_API) |
|
26abff55f6fe
optimize bsxfun for common built-in operations
Jaroslav Hajek <highegg@gmail.com>
parents:
9732
diff
changeset
|
191 |
| 4514 | 192 #endif |
| 193 | |
| 194 /* | |
| 195 ;;; Local Variables: *** | |
| 196 ;;; mode: C++ *** | |
| 197 ;;; End: *** | |
| 198 */ |