Mercurial > hg > octave-nkf
annotate libinterp/octave-value/ov-base-mat.h @ 20811:a22d8a2eb0e5
fix adaptive strategy in ode solvers.
* script/ode/ode45.m: remove unused option OutputSave
* script/ode/private/integrate_adaptive.m: rewrite algorithm
to be more compatible.
* script/ode/private/runge_kutta_45_dorpri.m: use kahan summation
for time increment.
| author | Carlo de Falco <carlo.defalco@polimi.it> |
|---|---|
| date | Sun, 11 Oct 2015 18:44:58 +0200 |
| parents | 4197fc428c7d |
| children |
| rev | line source |
|---|---|
| 3219 | 1 /* |
| 2 | |
|
19898
4197fc428c7d
maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents:
18484
diff
changeset
|
3 Copyright (C) 1998-2015 John W. Eaton |
| 11523 | 4 Copyright (C) 2009-2010 VZLU Prague |
| 3219 | 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. | |
| 3219 | 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/>. | |
| 3219 | 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_base_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_base_mat_h 1 |
| 3219 | 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> |
| 3219 | 30 #include <string> |
| 31 | |
| 32 #include "mx-base.h" | |
| 33 #include "str-vec.h" | |
| 6376 | 34 #include "MatrixType.h" |
| 3219 | 35 |
| 36 #include "error.h" | |
| 4649 | 37 #include "oct-obj.h" |
| 3219 | 38 #include "ov-base.h" |
| 39 #include "ov-typeinfo.h" | |
| 40 | |
| 41 class tree_walker; | |
| 42 | |
| 43 // Real matrix values. | |
| 44 | |
| 3223 | 45 template <class MT> |
| 46 class | |
| 3219 | 47 octave_base_matrix : public octave_base_value |
| 48 { | |
| 49 public: | |
| 50 | |
| 51 octave_base_matrix (void) | |
|
11584
cda4aa780d58
Another round of initialising members in the constructor initialisation list
Pascal Dupuis <Pascal.Dupuis@uclouvain.be>
parents:
11523
diff
changeset
|
52 : octave_base_value (), matrix (), typ (), idx_cache () { } |
| 3219 | 53 |
| 6376 | 54 octave_base_matrix (const MT& m, const MatrixType& t = MatrixType ()) |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11584
diff
changeset
|
55 : octave_base_value (), matrix (m), |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14557
diff
changeset
|
56 typ (t.is_known () ? new MatrixType (t) : 0), idx_cache () |
| 4676 | 57 { |
| 58 if (matrix.ndims () == 0) | |
| 59 matrix.resize (dim_vector (0, 0)); | |
| 60 } | |
| 3219 | 61 |
| 62 octave_base_matrix (const octave_base_matrix& m) | |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11584
diff
changeset
|
63 : octave_base_value (), matrix (m.matrix), |
|
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11584
diff
changeset
|
64 typ (m.typ ? new MatrixType (*m.typ) : 0), |
|
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11584
diff
changeset
|
65 idx_cache (m.idx_cache ? new idx_vector (*m.idx_cache) : 0) |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
66 { } |
| 3219 | 67 |
|
9350
16a5f9e1fdb3
cache idx_vector result in matrices once used for indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8994
diff
changeset
|
68 ~octave_base_matrix (void) { clear_cached_info (); } |
| 3219 | 69 |
| 4901 | 70 size_t byte_size (void) const { return matrix.byte_size (); } |
| 71 | |
| 5147 | 72 octave_value squeeze (void) const { return MT (matrix.squeeze ()); } |
| 4532 | 73 |
|
8458
d254a21e0120
reimplement full as method of octave_base_value
Jaroslav Hajek <highegg@gmail.com>
parents:
8150
diff
changeset
|
74 octave_value full_value (void) const { return matrix; } |
|
d254a21e0120
reimplement full as method of octave_base_value
Jaroslav Hajek <highegg@gmail.com>
parents:
8150
diff
changeset
|
75 |
|
8523
ad3afaaa19c1
implement non-copying contiguous range indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8458
diff
changeset
|
76 void maybe_economize (void) { matrix.maybe_economize (); } |
|
ad3afaaa19c1
implement non-copying contiguous range indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8458
diff
changeset
|
77 |
| 4247 | 78 octave_value subsref (const std::string& type, |
| 10313 | 79 const std::list<octave_value_list>& idx); |
| 3219 | 80 |
|
7651
443a8f5a50fd
require both subsref variants to be defined in octave_value subclasses
John W. Eaton <jwe@octave.org>
parents:
7620
diff
changeset
|
81 octave_value_list subsref (const std::string& type, |
| 10313 | 82 const std::list<octave_value_list>& idx, int) |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
83 { return subsref (type, idx); } |
| 4271 | 84 |
| 4247 | 85 octave_value subsasgn (const std::string& type, |
| 10313 | 86 const std::list<octave_value_list>& idx, |
| 87 const octave_value& rhs); | |
| 3933 | 88 |
| 5885 | 89 octave_value do_index_op (const octave_value_list& idx, |
| 10313 | 90 bool resize_ok = false); |
| 3220 | 91 |
|
13696
d6118a2c0644
fix indexing bug for matrices inside anonymous functions
John W. Eaton <jwe@octave.org>
parents:
12153
diff
changeset
|
92 octave_value_list do_multi_index_op (int, const octave_value_list& idx) |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
93 { return do_index_op (idx); } |
|
13696
d6118a2c0644
fix indexing bug for matrices inside anonymous functions
John W. Eaton <jwe@octave.org>
parents:
12153
diff
changeset
|
94 |
| 3928 | 95 void assign (const octave_value_list& idx, const MT& rhs); |
| 96 | |
|
8679
280fae940bb0
optimize scalar indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8523
diff
changeset
|
97 void assign (const octave_value_list& idx, typename MT::element_type rhs); |
|
280fae940bb0
optimize scalar indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8523
diff
changeset
|
98 |
|
8150
283989f2da9b
make null assignment matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
7651
diff
changeset
|
99 void delete_elements (const octave_value_list& idx); |
|
283989f2da9b
make null assignment matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
7651
diff
changeset
|
100 |
| 4513 | 101 dim_vector dims (void) const { return matrix.dims (); } |
| 102 | |
| 8994 | 103 octave_idx_type numel (void) const { return matrix.numel (); } |
| 104 | |
|
10653
ec5fa46e0e45
override ndims and numel for scalars and matrices
Jaroslav Hajek <highegg@gmail.com>
parents:
10544
diff
changeset
|
105 int ndims (void) const { return matrix.ndims (); } |
|
ec5fa46e0e45
override ndims and numel for scalars and matrices
Jaroslav Hajek <highegg@gmail.com>
parents:
10544
diff
changeset
|
106 |
| 5602 | 107 octave_idx_type nnz (void) const { return matrix.nnz (); } |
| 108 | |
| 4567 | 109 octave_value reshape (const dim_vector& new_dims) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
110 { return MT (matrix.reshape (new_dims)); } |
| 4567 | 111 |
| 4593 | 112 octave_value permute (const Array<int>& vec, bool inv = false) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
113 { return MT (matrix.permute (vec, inv)); } |
| 4593 | 114 |
| 5731 | 115 octave_value resize (const dim_vector& dv, bool fill = false) const; |
| 4915 | 116 |
| 4017 | 117 octave_value all (int dim = 0) const { return matrix.all (dim); } |
| 118 octave_value any (int dim = 0) const { return matrix.any (dim); } | |
| 3221 | 119 |
|
9350
16a5f9e1fdb3
cache idx_vector result in matrices once used for indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8994
diff
changeset
|
120 MatrixType matrix_type (void) const { return typ ? *typ : MatrixType (); } |
|
16a5f9e1fdb3
cache idx_vector result in matrices once used for indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8994
diff
changeset
|
121 MatrixType matrix_type (const MatrixType& _typ) const; |
| 5785 | 122 |
|
7620
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
123 octave_value diag (octave_idx_type k = 0) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
124 { return octave_value (matrix.diag (k)); } |
|
7620
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7463
diff
changeset
|
125 |
|
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
|
126 octave_value diag (octave_idx_type m, octave_idx_type n) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
127 { return octave_value (matrix.diag (m, n)); } |
|
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
|
128 |
|
7463
2467639bd8c0
eliminate UNDEFINED sort mode
John W. Eaton <jwe@octave.org>
parents:
7433
diff
changeset
|
129 octave_value sort (octave_idx_type dim = 0, sortmode mode = ASCENDING) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
130 { return octave_value (matrix.sort (dim, mode)); } |
| 7433 | 131 octave_value sort (Array<octave_idx_type> &sidx, octave_idx_type dim = 0, |
| 10313 | 132 sortmode mode = ASCENDING) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
133 { return octave_value (matrix.sort (sidx, dim, mode)); } |
| 7433 | 134 |
|
8734
767ed8cc6634
rename internal issorted and issorted_rows functions to is_sorted and is_sorted_rows
John W. Eaton <jwe@octave.org>
parents:
8733
diff
changeset
|
135 sortmode is_sorted (sortmode mode = UNSORTED) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
136 { return matrix.is_sorted (mode); } |
|
8721
e9cb742df9eb
imported patch sort3.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
8679
diff
changeset
|
137 |
|
8733
3ef774603887
rename all uses of sortrows_idx to sort_rows_idx
John W. Eaton <jwe@octave.org>
parents:
8721
diff
changeset
|
138 Array<octave_idx_type> sort_rows_idx (sortmode mode = ASCENDING) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
139 { return matrix.sort_rows_idx (mode); } |
|
8721
e9cb742df9eb
imported patch sort3.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
8679
diff
changeset
|
140 |
|
8734
767ed8cc6634
rename internal issorted and issorted_rows functions to is_sorted and is_sorted_rows
John W. Eaton <jwe@octave.org>
parents:
8733
diff
changeset
|
141 sortmode is_sorted_rows (sortmode mode = UNSORTED) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
142 { return matrix.is_sorted_rows (mode); } |
|
8721
e9cb742df9eb
imported patch sort3.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
8679
diff
changeset
|
143 |
| 3221 | 144 bool is_matrix_type (void) const { return true; } |
| 145 | |
| 146 bool is_numeric_type (void) const { return true; } | |
| 147 | |
| 3219 | 148 bool is_defined (void) const { return true; } |
| 149 | |
| 150 bool is_constant (void) const { return true; } | |
| 151 | |
| 3220 | 152 bool is_true (void) const; |
| 153 | |
| 4604 | 154 bool print_as_scalar (void) const; |
| 3219 | 155 |
|
18484
bcd71a2531d3
Support disp/display overloading in classdef
Michael Goffioul <michael.goffioul@gmail.com>
parents:
17870
diff
changeset
|
156 void print (std::ostream& os, bool pr_as_read_syntax = false); |
| 3219 | 157 |
| 3933 | 158 void print_info (std::ostream& os, const std::string& prefix) const; |
| 159 | |
| 17870 | 160 void short_disp (std::ostream& os) const; |
|
17866
ea0ecbe2eaf5
display matrix values in GUI workspace viewer (bug #40499)
John W. Eaton <jwe@octave.org>
parents:
17822
diff
changeset
|
161 |
|
9546
1beb23d2b892
optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents:
9358
diff
changeset
|
162 MT& matrix_ref (void) |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
163 { |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
164 clear_cached_info (); |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
165 return matrix; |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
166 } |
|
9546
1beb23d2b892
optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents:
9358
diff
changeset
|
167 |
|
10665
0f310fce905d
implement octave_base_matrix::matrix_ref const
Jaroslav Hajek <highegg@gmail.com>
parents:
10653
diff
changeset
|
168 const MT& matrix_ref (void) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
169 { |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
170 return matrix; |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
171 } |
|
10665
0f310fce905d
implement octave_base_matrix::matrix_ref const
Jaroslav Hajek <highegg@gmail.com>
parents:
10653
diff
changeset
|
172 |
|
10670
654fbde5dceb
make cellfun's fast scalar collection mechanism public
Jaroslav Hajek <highegg@gmail.com>
parents:
10665
diff
changeset
|
173 octave_value |
|
654fbde5dceb
make cellfun's fast scalar collection mechanism public
Jaroslav Hajek <highegg@gmail.com>
parents:
10665
diff
changeset
|
174 fast_elem_extract (octave_idx_type n) const; |
|
654fbde5dceb
make cellfun's fast scalar collection mechanism public
Jaroslav Hajek <highegg@gmail.com>
parents:
10665
diff
changeset
|
175 |
|
654fbde5dceb
make cellfun's fast scalar collection mechanism public
Jaroslav Hajek <highegg@gmail.com>
parents:
10665
diff
changeset
|
176 bool |
|
654fbde5dceb
make cellfun's fast scalar collection mechanism public
Jaroslav Hajek <highegg@gmail.com>
parents:
10665
diff
changeset
|
177 fast_elem_insert (octave_idx_type n, const octave_value& x); |
|
654fbde5dceb
make cellfun's fast scalar collection mechanism public
Jaroslav Hajek <highegg@gmail.com>
parents:
10665
diff
changeset
|
178 |
| 3219 | 179 protected: |
| 180 | |
| 181 MT matrix; | |
| 5785 | 182 |
|
9350
16a5f9e1fdb3
cache idx_vector result in matrices once used for indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8994
diff
changeset
|
183 idx_vector set_idx_cache (const idx_vector& idx) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
184 { |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
185 delete idx_cache; |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
186 idx_cache = idx ? new idx_vector (idx) : 0; |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
187 return idx; |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
188 } |
|
9350
16a5f9e1fdb3
cache idx_vector result in matrices once used for indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8994
diff
changeset
|
189 |
|
16a5f9e1fdb3
cache idx_vector result in matrices once used for indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8994
diff
changeset
|
190 void clear_cached_info (void) const |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
191 { |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
192 delete typ; typ = 0; |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
193 delete idx_cache; idx_cache = 0; |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
194 } |
|
9350
16a5f9e1fdb3
cache idx_vector result in matrices once used for indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8994
diff
changeset
|
195 |
|
16a5f9e1fdb3
cache idx_vector result in matrices once used for indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8994
diff
changeset
|
196 mutable MatrixType *typ; |
|
16a5f9e1fdb3
cache idx_vector result in matrices once used for indexing
Jaroslav Hajek <highegg@gmail.com>
parents:
8994
diff
changeset
|
197 mutable idx_vector *idx_cache; |
|
12153
e0e50f48df37
Explicitly disallow copying in some classes
Pascal Dupuis <Pascal.Dupuis@uclouvain.be>
parents:
11586
diff
changeset
|
198 |
|
e0e50f48df37
Explicitly disallow copying in some classes
Pascal Dupuis <Pascal.Dupuis@uclouvain.be>
parents:
11586
diff
changeset
|
199 private: |
|
e0e50f48df37
Explicitly disallow copying in some classes
Pascal Dupuis <Pascal.Dupuis@uclouvain.be>
parents:
11586
diff
changeset
|
200 |
|
e0e50f48df37
Explicitly disallow copying in some classes
Pascal Dupuis <Pascal.Dupuis@uclouvain.be>
parents:
11586
diff
changeset
|
201 // No assignment. |
|
e0e50f48df37
Explicitly disallow copying in some classes
Pascal Dupuis <Pascal.Dupuis@uclouvain.be>
parents:
11586
diff
changeset
|
202 |
|
e0e50f48df37
Explicitly disallow copying in some classes
Pascal Dupuis <Pascal.Dupuis@uclouvain.be>
parents:
11586
diff
changeset
|
203 octave_base_matrix& operator = (const octave_base_matrix&); |
| 3219 | 204 }; |
| 205 | |
| 206 #endif |