Mercurial > octave-nkf
annotate liboctave/CNDArray.h @ 7789:82be108cc558
First attempt at single precision tyeps
* * *
corrections to qrupdate single precision routines
* * *
prefer demotion to single over promotion to double
* * *
Add single precision support to log2 function
* * *
Trivial PROJECT file update
* * *
Cache optimized hermitian/transpose methods
* * *
Add tests for tranpose/hermitian and ChangeLog entry for new transpose code
author | David Bateman <dbateman@free.fr> |
---|---|
date | Sun, 27 Apr 2008 22:34:17 +0200 |
parents | 36594d5bbe13 |
children | 935be827eaf8 |
rev | line source |
---|---|
4514 | 1 /* |
2 | |
7017 | 3 Copyright (C) 2003, 2004, 2005, 2006, 2007 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" | |
30 #include "mx-op-defs.h" | |
31 | |
32 class | |
6108 | 33 OCTAVE_API |
4514 | 34 ComplexNDArray : public MArrayN<Complex> |
35 { | |
36 public: | |
5523 | 37 |
4514 | 38 ComplexNDArray (void) : MArrayN<Complex> () { } |
39 | |
4587 | 40 ComplexNDArray (const dim_vector& dv) : MArrayN<Complex> (dv) { } |
4514 | 41 |
4587 | 42 ComplexNDArray (const dim_vector& dv, const Complex& val) |
43 : MArrayN<Complex> (dv, val) { } | |
4514 | 44 |
45 ComplexNDArray (const ComplexNDArray& a) : MArrayN<Complex> (a) { } | |
46 | |
47 ComplexNDArray (const ComplexMatrix& a) : MArrayN<Complex> (a) { } | |
48 | |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
49 template <class U> |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
50 ComplexNDArray (const MArrayN<U>& a) : MArrayN<Complex> (a) { } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7620
diff
changeset
|
51 |
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 ArrayN<U>& a) : MArrayN<Complex> (a) { } |
4514 | 54 |
55 ComplexNDArray& operator = (const ComplexNDArray& a) | |
56 { | |
57 MArrayN<Complex>::operator = (a); | |
58 return *this; | |
59 } | |
60 | |
4543 | 61 // unary operations |
62 | |
63 boolNDArray operator ! (void) const; | |
64 | |
5775 | 65 // FIXME -- this is not quite the right thing. |
4514 | 66 |
4687 | 67 bool any_element_is_inf_or_nan (void) const; |
68 bool all_elements_are_real (void) const; | |
69 bool all_integers (double& max_val, double& min_val) const; | |
70 bool too_large_for_float (void) const; | |
71 | |
4556 | 72 boolNDArray all (int dim = -1) const; |
73 boolNDArray any (int dim = -1) const; | |
4514 | 74 |
4584 | 75 ComplexNDArray cumprod (int dim = -1) const; |
76 ComplexNDArray cumsum (int dim = -1) const; | |
4569 | 77 ComplexNDArray prod (int dim = -1) const; |
78 ComplexNDArray sum (int dim = -1) const; | |
79 ComplexNDArray sumsq (int dim = -1) const; | |
5275 | 80 ComplexNDArray concat (const ComplexNDArray& rb, const Array<octave_idx_type>& ra_idx); |
81 ComplexNDArray concat (const NDArray& rb, const Array<octave_idx_type>& ra_idx); | |
4569 | 82 |
4844 | 83 ComplexNDArray max (int dim = 0) const; |
5275 | 84 ComplexNDArray max (ArrayN<octave_idx_type>& index, int dim = 0) const; |
4844 | 85 ComplexNDArray min (int dim = 0) const; |
5275 | 86 ComplexNDArray min (ArrayN<octave_idx_type>& index, int dim = 0) const; |
87 ComplexNDArray& insert (const NDArray& a, octave_idx_type r, octave_idx_type c); | |
88 ComplexNDArray& insert (const ComplexNDArray& a, octave_idx_type r, octave_idx_type c); | |
89 ComplexNDArray& insert (const ComplexNDArray& a, const Array<octave_idx_type>& ra_idx); | |
4765 | 90 |
4634 | 91 NDArray abs (void) const; |
4569 | 92 |
4773 | 93 ComplexNDArray fourier (int dim = 1) const; |
94 ComplexNDArray ifourier (int dim = 1) const; | |
95 | |
96 ComplexNDArray fourier2d (void) const; | |
97 ComplexNDArray ifourier2d (void) const; | |
98 | |
99 ComplexNDArray fourierNd (void) const; | |
100 ComplexNDArray ifourierNd (void) const; | |
101 | |
4514 | 102 ComplexMatrix matrix_value (void) const; |
103 | |
4902 | 104 ComplexNDArray squeeze (void) const { return MArrayN<Complex>::squeeze (); } |
4532 | 105 |
5275 | 106 static void increment_index (Array<octave_idx_type>& ra_idx, |
4532 | 107 const dim_vector& dimensions, |
108 int start_dimension = 0); | |
109 | |
5275 | 110 static octave_idx_type compute_index (Array<octave_idx_type>& ra_idx, |
4556 | 111 const dim_vector& dimensions); |
112 | |
4514 | 113 // i/o |
114 | |
6108 | 115 friend OCTAVE_API std::ostream& operator << (std::ostream& os, const ComplexNDArray& a); |
116 friend OCTAVE_API std::istream& operator >> (std::istream& is, ComplexNDArray& a); | |
4514 | 117 |
118 static Complex resize_fill_value (void) { return Complex (0.0, 0.0); } | |
119 | |
120 // bool all_elements_are_real (void) const; | |
121 // bool all_integers (double& max_val, double& min_val) const; | |
122 | |
7620
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7503
diff
changeset
|
123 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
|
124 |
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
125 typedef double (*dmapper) (const Complex&); |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
126 typedef Complex (*cmapper) (const Complex&); |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
127 typedef bool (*bmapper) (const Complex&); |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
128 |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
129 NDArray map (dmapper fcn) const; |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
130 ComplexNDArray map (cmapper fcn) const; |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
131 boolNDArray map (bmapper fcn) const; |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
132 |
4514 | 133 private: |
134 | |
4634 | 135 ComplexNDArray (Complex *d, const dim_vector& dv) |
4587 | 136 : MArrayN<Complex> (d, dv) { } |
4514 | 137 }; |
138 | |
6108 | 139 extern OCTAVE_API ComplexNDArray min (const Complex& c, const ComplexNDArray& m); |
140 extern OCTAVE_API ComplexNDArray min (const ComplexNDArray& m, const Complex& c); | |
141 extern OCTAVE_API ComplexNDArray min (const ComplexNDArray& a, const ComplexNDArray& b); | |
4844 | 142 |
6108 | 143 extern OCTAVE_API ComplexNDArray max (const Complex& c, const ComplexNDArray& m); |
144 extern OCTAVE_API ComplexNDArray max (const ComplexNDArray& m, const Complex& c); | |
145 extern OCTAVE_API ComplexNDArray max (const ComplexNDArray& a, const ComplexNDArray& b); | |
4844 | 146 |
6708 | 147 NDS_CMP_OP_DECLS (ComplexNDArray, Complex, OCTAVE_API) |
148 NDS_BOOL_OP_DECLS (ComplexNDArray, Complex, OCTAVE_API) | |
4543 | 149 |
6708 | 150 SND_CMP_OP_DECLS (Complex, ComplexNDArray, OCTAVE_API) |
151 SND_BOOL_OP_DECLS (Complex, ComplexNDArray, OCTAVE_API) | |
4543 | 152 |
6708 | 153 NDND_CMP_OP_DECLS (ComplexNDArray, ComplexNDArray, OCTAVE_API) |
154 NDND_BOOL_OP_DECLS (ComplexNDArray, ComplexNDArray, OCTAVE_API) | |
4543 | 155 |
4514 | 156 MARRAY_FORWARD_DEFS (MArrayN, ComplexNDArray, Complex) |
157 | |
158 #endif | |
159 | |
160 /* | |
161 ;;; Local Variables: *** | |
162 ;;; mode: C++ *** | |
163 ;;; End: *** | |
164 */ |