Mercurial > hg > octave-nkf
annotate libinterp/octave-value/ov-re-mat.h @ 20830:b65888ec820e draft default tip gccjit
dmalcom gcc jit import
| author | Stefan Mahr <dac922@gmx.de> |
|---|---|
| date | Fri, 27 Feb 2015 16:59:36 +0100 |
| parents | 09ed6f7538dd |
| children |
| rev | line source |
|---|---|
| 2376 | 1 /* |
| 2 | |
|
19898
4197fc428c7d
maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents:
19795
diff
changeset
|
3 Copyright (C) 1996-2015 John W. Eaton |
| 11523 | 4 Copyright (C) 2009-2010 VZLU Prague |
| 2376 | 5 |
| 6 This file is part of Octave. | |
| 7 | |
| 8 Octave is free software; you can redistribute it and/or modify it | |
| 9 under the terms of the GNU General Public License as published by the | |
| 7016 | 10 Free Software Foundation; either version 3 of the License, or (at your |
| 11 option) any later version. | |
| 2376 | 12 |
| 13 Octave is distributed in the hope that it will be useful, but WITHOUT | |
| 14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
| 15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
| 16 for more details. | |
| 17 | |
| 18 You should have received a copy of the GNU General Public License | |
| 7016 | 19 along with Octave; see the file COPYING. If not, see |
| 20 <http://www.gnu.org/licenses/>. | |
| 2376 | 21 |
| 22 */ | |
| 23 | |
|
17822
ebb3ef964372
maint: Use common #define syntax "octave_filename_h" in h_files.
Rik <rik@octave.org>
parents:
17787
diff
changeset
|
24 #if !defined (octave_ov_re_mat_h) |
|
ebb3ef964372
maint: Use common #define syntax "octave_filename_h" in h_files.
Rik <rik@octave.org>
parents:
17787
diff
changeset
|
25 #define octave_ov_re_mat_h 1 |
| 2376 | 26 |
| 27 #include <cstdlib> | |
| 28 | |
|
8950
d865363208d6
include <iosfwd> instead of <iostream> in header files
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
29 #include <iosfwd> |
| 2376 | 30 #include <string> |
| 31 | |
| 32 #include "mx-base.h" | |
| 33 #include "str-vec.h" | |
| 34 | |
| 35 #include "error.h" | |
| 4944 | 36 #include "oct-stream.h" |
| 2376 | 37 #include "ov-base.h" |
| 3219 | 38 #include "ov-base-mat.h" |
| 2376 | 39 #include "ov-typeinfo.h" |
| 40 | |
| 5785 | 41 #include "MatrixType.h" |
| 42 | |
| 2376 | 43 class octave_value_list; |
| 44 | |
| 45 class tree_walker; | |
| 46 | |
| 47 // Real matrix values. | |
| 48 | |
| 49 class | |
| 6153 | 50 OCTINTERP_API |
| 4513 | 51 octave_matrix : public octave_base_matrix<NDArray> |
| 2376 | 52 { |
| 53 public: | |
| 54 | |
| 55 octave_matrix (void) | |
| 4513 | 56 : octave_base_matrix<NDArray> () { } |
| 2376 | 57 |
| 58 octave_matrix (const Matrix& m) | |
| 4513 | 59 : octave_base_matrix<NDArray> (m) { } |
| 60 | |
| 5785 | 61 octave_matrix (const Matrix& m, const MatrixType& t) |
| 62 : octave_base_matrix<NDArray> (m, t) { } | |
| 63 | |
| 4513 | 64 octave_matrix (const NDArray& nda) |
| 65 : octave_base_matrix<NDArray> (nda) { } | |
| 2376 | 66 |
|
9732
b4fdfee405b5
remove ArrayN<T> + fix nonhom. diag-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents:
9685
diff
changeset
|
67 octave_matrix (const Array<double>& m) |
| 4911 | 68 : octave_base_matrix<NDArray> (NDArray (m)) { } |
| 69 | |
| 2376 | 70 octave_matrix (const DiagMatrix& d) |
| 4513 | 71 : octave_base_matrix<NDArray> (Matrix (d)) { } |
| 2376 | 72 |
| 3418 | 73 octave_matrix (const RowVector& v) |
| 4513 | 74 : octave_base_matrix<NDArray> (Matrix (v)) { } |
| 2376 | 75 |
| 3418 | 76 octave_matrix (const ColumnVector& v) |
| 4513 | 77 : octave_base_matrix<NDArray> (Matrix (v)) { } |
| 2376 | 78 |
| 79 octave_matrix (const octave_matrix& m) | |
| 4513 | 80 : octave_base_matrix<NDArray> (m) { } |
| 2376 | 81 |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
82 octave_matrix (const Array<octave_idx_type>& idx, |
|
9351
e2344f4af0cb
autocache indices returned from find et al.
Jaroslav Hajek <highegg@gmail.com>
parents:
9350
diff
changeset
|
83 bool zero_based = false, bool cache_index = false) |
|
e2344f4af0cb
autocache indices returned from find et al.
Jaroslav Hajek <highegg@gmail.com>
parents:
9350
diff
changeset
|
84 : octave_base_matrix<NDArray> (NDArray (idx, zero_based)) |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
85 { |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
86 // Auto-create cache to speed up subsequent indexing. |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
87 if (zero_based && cache_index) |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
88 set_idx_cache (idx_vector (idx)); |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
89 } |
|
9351
e2344f4af0cb
autocache indices returned from find et al.
Jaroslav Hajek <highegg@gmail.com>
parents:
9350
diff
changeset
|
90 |
|
9479
d9716e3ee0dd
supply optimized compiled sub2ind & ind2sub
Jaroslav Hajek <highegg@gmail.com>
parents:
9358
diff
changeset
|
91 octave_matrix (const NDArray& nda, const idx_vector& cache) |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
92 : octave_base_matrix<NDArray> (nda) |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
93 { |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
94 set_idx_cache (cache); |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
95 } |
|
9479
d9716e3ee0dd
supply optimized compiled sub2ind & ind2sub
Jaroslav Hajek <highegg@gmail.com>
parents:
9358
diff
changeset
|
96 |
| 2376 | 97 ~octave_matrix (void) { } |
| 98 | |
| 5759 | 99 octave_base_value *clone (void) const { return new octave_matrix (*this); } |
| 100 octave_base_value *empty_clone (void) const { return new octave_matrix (); } | |
| 2376 | 101 |
|
8345
c777f3ce02d8
smarter conversion lookup
Jaroslav Hajek <highegg@gmail.com>
parents:
7885
diff
changeset
|
102 type_conv_info numeric_demotion_function (void) const; |
|
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
103 |
| 5759 | 104 octave_base_value *try_narrowing_conversion (void); |
| 2410 | 105 |
|
18129
e473c4853afc
enable non-integer ranges as indices by default and deprecate preference
John W. Eaton <jwe@octave.org>
parents:
17822
diff
changeset
|
106 idx_vector index_vector (bool /* require_integers */ = false) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
107 { return idx_cache ? *idx_cache : set_idx_cache (idx_vector (matrix)); } |
| 2376 | 108 |
| 9685 | 109 builtin_type_t builtin_type (void) const { return btyp_double; } |
| 110 | |
| 2376 | 111 bool is_real_matrix (void) const { return true; } |
| 112 | |
| 113 bool is_real_type (void) const { return true; } | |
| 114 | |
| 5895 | 115 bool is_double_type (void) const { return true; } |
| 116 | |
| 7576 | 117 bool is_float_type (void) const { return true; } |
| 118 | |
| 5533 | 119 int8NDArray |
| 120 int8_array_value (void) const { return int8NDArray (matrix); } | |
| 121 | |
| 122 int16NDArray | |
| 123 int16_array_value (void) const { return int16NDArray (matrix); } | |
| 124 | |
| 125 int32NDArray | |
| 126 int32_array_value (void) const { return int32NDArray (matrix); } | |
| 127 | |
| 128 int64NDArray | |
| 129 int64_array_value (void) const { return int64NDArray (matrix); } | |
| 130 | |
| 131 uint8NDArray | |
| 132 uint8_array_value (void) const { return uint8NDArray (matrix); } | |
| 133 | |
| 134 uint16NDArray | |
| 135 uint16_array_value (void) const { return uint16NDArray (matrix); } | |
| 136 | |
| 137 uint32NDArray | |
| 138 uint32_array_value (void) const { return uint32NDArray (matrix); } | |
| 139 | |
| 140 uint64NDArray | |
| 141 uint64_array_value (void) const { return uint64NDArray (matrix); } | |
| 142 | |
| 2376 | 143 double double_value (bool = false) const; |
| 144 | |
|
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
145 float float_value (bool = false) const; |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
146 |
| 3145 | 147 double scalar_value (bool frc_str_conv = false) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
148 { return double_value (frc_str_conv); } |
| 2916 | 149 |
| 4513 | 150 Matrix matrix_value (bool = false) const; |
| 2376 | 151 |
|
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
152 FloatMatrix float_matrix_value (bool = false) const; |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
153 |
| 2376 | 154 Complex complex_value (bool = false) const; |
| 155 | |
|
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
156 FloatComplex float_complex_value (bool = false) const; |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
157 |
| 4513 | 158 ComplexMatrix complex_matrix_value (bool = false) const; |
| 2376 | 159 |
|
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
160 FloatComplexMatrix float_complex_matrix_value (bool = false) const; |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
161 |
| 4699 | 162 ComplexNDArray complex_array_value (bool = false) const; |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
163 |
|
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
164 FloatComplexNDArray float_complex_array_value (bool = false) const; |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
165 |
| 5943 | 166 boolNDArray bool_array_value (bool warn = false) const; |
| 5898 | 167 |
| 4758 | 168 charNDArray char_array_value (bool = false) const; |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
169 |
| 4550 | 170 NDArray array_value (bool = false) const { return matrix; } |
| 4505 | 171 |
|
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
172 FloatNDArray float_array_value (bool = false) const { return matrix; } |
|
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
173 |
| 5164 | 174 SparseMatrix sparse_matrix_value (bool = false) const; |
| 175 | |
| 176 SparseComplexMatrix sparse_complex_matrix_value (bool = false) const; | |
| 177 | |
|
8366
8b1a2555c4e2
implement diagonal matrix objects
Jaroslav Hajek <highegg@gmail.com>
parents:
8345
diff
changeset
|
178 octave_value diag (octave_idx_type k = 0) const; |
|
8b1a2555c4e2
implement diagonal matrix objects
Jaroslav Hajek <highegg@gmail.com>
parents:
8345
diff
changeset
|
179 |
|
14557
e8e86ae3abbc
make diag (x, m, n) return a proper diagonal matrix object (bug #36099)
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
180 octave_value diag (octave_idx_type m, octave_idx_type n) const; |
|
e8e86ae3abbc
make diag (x, m, n) return a proper diagonal matrix object (bug #36099)
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
181 |
|
10339
de2d43bcb083
optimize some lazy index operations
Jaroslav Hajek <highegg@gmail.com>
parents:
10313
diff
changeset
|
182 octave_value reshape (const dim_vector& new_dims) const; |
|
de2d43bcb083
optimize some lazy index operations
Jaroslav Hajek <highegg@gmail.com>
parents:
10313
diff
changeset
|
183 |
|
de2d43bcb083
optimize some lazy index operations
Jaroslav Hajek <highegg@gmail.com>
parents:
10313
diff
changeset
|
184 octave_value squeeze (void) const; |
|
de2d43bcb083
optimize some lazy index operations
Jaroslav Hajek <highegg@gmail.com>
parents:
10313
diff
changeset
|
185 |
|
10273
3a8c13b71612
implement special-case optimization for sort of index vectors
Jaroslav Hajek <highegg@gmail.com>
parents:
10160
diff
changeset
|
186 octave_value sort (octave_idx_type dim = 0, sortmode mode = ASCENDING) const; |
|
3a8c13b71612
implement special-case optimization for sort of index vectors
Jaroslav Hajek <highegg@gmail.com>
parents:
10160
diff
changeset
|
187 octave_value sort (Array<octave_idx_type> &sidx, octave_idx_type dim = 0, |
| 10313 | 188 sortmode mode = ASCENDING) const; |
|
10273
3a8c13b71612
implement special-case optimization for sort of index vectors
Jaroslav Hajek <highegg@gmail.com>
parents:
10160
diff
changeset
|
189 |
|
10339
de2d43bcb083
optimize some lazy index operations
Jaroslav Hajek <highegg@gmail.com>
parents:
10313
diff
changeset
|
190 sortmode is_sorted (sortmode mode = UNSORTED) const; |
|
de2d43bcb083
optimize some lazy index operations
Jaroslav Hajek <highegg@gmail.com>
parents:
10313
diff
changeset
|
191 |
|
de2d43bcb083
optimize some lazy index operations
Jaroslav Hajek <highegg@gmail.com>
parents:
10313
diff
changeset
|
192 Array<octave_idx_type> sort_rows_idx (sortmode mode = ASCENDING) const; |
|
de2d43bcb083
optimize some lazy index operations
Jaroslav Hajek <highegg@gmail.com>
parents:
10313
diff
changeset
|
193 |
|
de2d43bcb083
optimize some lazy index operations
Jaroslav Hajek <highegg@gmail.com>
parents:
10313
diff
changeset
|
194 sortmode is_sorted_rows (sortmode mode = UNSORTED) const; |
|
de2d43bcb083
optimize some lazy index operations
Jaroslav Hajek <highegg@gmail.com>
parents:
10313
diff
changeset
|
195 |
|
9606
a04352386a6b
clear index cache on ++,-- operators
Jaroslav Hajek <highegg@gmail.com>
parents:
9479
diff
changeset
|
196 // Use matrix_ref here to clear index cache. |
|
a04352386a6b
clear index cache on ++,-- operators
Jaroslav Hajek <highegg@gmail.com>
parents:
9479
diff
changeset
|
197 void increment (void) { matrix_ref () += 1.0; } |
| 2376 | 198 |
|
9606
a04352386a6b
clear index cache on ++,-- operators
Jaroslav Hajek <highegg@gmail.com>
parents:
9479
diff
changeset
|
199 void decrement (void) { matrix_ref () -= 1.0; } |
| 2376 | 200 |
|
9607
1be3c73ed7b5
reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents:
9606
diff
changeset
|
201 void changesign (void) { matrix_ref ().changesign (); } |
|
1be3c73ed7b5
reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents:
9606
diff
changeset
|
202 |
| 5279 | 203 octave_value convert_to_str_internal (bool pad, bool force, char type) const; |
| 2376 | 204 |
| 4643 | 205 void print_raw (std::ostream& os, bool pr_as_read_syntax = false) const; |
| 206 | |
| 6974 | 207 bool save_ascii (std::ostream& os); |
| 4687 | 208 |
| 209 bool load_ascii (std::istream& is); | |
| 210 | |
| 211 bool save_binary (std::ostream& os, bool& save_as_floats); | |
| 212 | |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
213 bool load_binary (std::istream& is, bool swap, |
| 10313 | 214 oct_mach_info::float_format fmt); |
| 4687 | 215 |
|
20070
09ed6f7538dd
avoid needing to include hdf5 in public header files (bug #44370, #43180)
John W. Eaton <jwe@octave.org> and Mike Miller <mtmiller@ieee.org>
parents:
19898
diff
changeset
|
216 bool save_hdf5 (octave_hdf5_id loc_id, const char *name, bool save_as_floats); |
| 4687 | 217 |
|
20070
09ed6f7538dd
avoid needing to include hdf5 in public header files (bug #44370, #43180)
John W. Eaton <jwe@octave.org> and Mike Miller <mtmiller@ieee.org>
parents:
19898
diff
changeset
|
218 bool load_hdf5 (octave_hdf5_id loc_id, const char *name); |
| 4687 | 219 |
| 4944 | 220 int write (octave_stream& os, int block_size, |
| 10313 | 221 oct_data_conv::data_type output_type, int skip, |
| 222 oct_mach_info::float_format flt_fmt) const | |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
223 { return os.write (matrix, block_size, output_type, skip, flt_fmt); } |
| 4944 | 224 |
|
9358
d4b1314a7c31
mex.cc (mxArray_octave_value::get_data): avoid enumerating types that can be handled as foreign
John W. Eaton <jwe@octave.org>
parents:
9351
diff
changeset
|
225 // Unsafe. This function exists to support the MEX interface. |
|
d4b1314a7c31
mex.cc (mxArray_octave_value::get_data): avoid enumerating types that can be handled as foreign
John W. Eaton <jwe@octave.org>
parents:
9351
diff
changeset
|
226 // You should not use it anywhere else. |
|
d4b1314a7c31
mex.cc (mxArray_octave_value::get_data): avoid enumerating types that can be handled as foreign
John W. Eaton <jwe@octave.org>
parents:
9351
diff
changeset
|
227 void *mex_get_data (void) const { return matrix.mex_get_data (); } |
|
d4b1314a7c31
mex.cc (mxArray_octave_value::get_data): avoid enumerating types that can be handled as foreign
John W. Eaton <jwe@octave.org>
parents:
9351
diff
changeset
|
228 |
| 5900 | 229 mxArray *as_mxArray (void) const; |
| 230 | |
|
9812
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
9732
diff
changeset
|
231 octave_value map (unary_mapper_t umap) const; |
|
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
232 |
| 2376 | 233 private: |
| 2477 | 234 |
| 3219 | 235 DECLARE_OV_TYPEID_FUNCTIONS_AND_DATA |
| 2376 | 236 }; |
| 237 | |
| 238 #endif |