Mercurial > octave
annotate liboctave/array/dim-vector.cc @ 22213:1723063e65f3
dim-vector: add missing, and remove unecessary, #includes.
author | Carnë Draug <carandraug@octave.org> |
---|---|
date | Sat, 06 Aug 2016 05:18:34 +0100 |
parents | 2b487dbe419d |
children | bac0d6f07a3e |
rev | line source |
---|---|
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
1 /* |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
2 |
19697
4197fc428c7d
maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents:
18084
diff
changeset
|
3 Copyright (C) 2003-2015 John W. Eaton |
10686
8675858ba1e2
fix copyright stmt in 10683:4fdb36ffa066
Jaroslav Hajek <highegg@gmail.com>
parents:
10683
diff
changeset
|
4 Copyirght (C) 2009, 2010 VZLU Prague |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
5 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
6 This file is part of Octave. |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
7 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
8 Octave is free software; you can redistribute it and/or modify it |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
9 under the terms of the GNU General Public License as published by the |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
10 Free Software Foundation; either version 3 of the License, or (at your |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
11 option) any later version. |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
12 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
13 Octave is distributed in the hope that it will be useful, but WITHOUT |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
16 for more details. |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
17 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
18 You should have received a copy of the GNU General Public License |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
19 along with Octave; see the file COPYING. If not, see |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
20 <http://www.gnu.org/licenses/>. |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
21 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
22 */ |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
23 |
21724
aba2e6293dd8
use "#if ..." consistently instead of "#ifdef" and "#ifndef"
John W. Eaton <jwe@octave.org>
parents:
21660
diff
changeset
|
24 #if defined (HAVE_CONFIG_H) |
21301
40de9f8f23a6
Use '#include "config.h"' rather than <config.h>.
Rik <rik@octave.org>
parents:
21202
diff
changeset
|
25 # include "config.h" |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
26 #endif |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
27 |
22213
1723063e65f3
dim-vector: add missing, and remove unecessary, #includes.
Carnë Draug <carandraug@octave.org>
parents:
22212
diff
changeset
|
28 #include <sstream> |
1723063e65f3
dim-vector: add missing, and remove unecessary, #includes.
Carnë Draug <carandraug@octave.org>
parents:
22212
diff
changeset
|
29 #include <limits> |
1723063e65f3
dim-vector: add missing, and remove unecessary, #includes.
Carnë Draug <carandraug@octave.org>
parents:
22212
diff
changeset
|
30 #include <string> |
1723063e65f3
dim-vector: add missing, and remove unecessary, #includes.
Carnë Draug <carandraug@octave.org>
parents:
22212
diff
changeset
|
31 #include <new> |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
32 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
33 #include "dim-vector.h" |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
34 |
21573
f3f8e1d3e399
avoid mulitple definitions of static function-scope vars (bug #47372)
John W. Eaton <jwe@octave.org>
parents:
19697
diff
changeset
|
35 octave_idx_type * |
f3f8e1d3e399
avoid mulitple definitions of static function-scope vars (bug #47372)
John W. Eaton <jwe@octave.org>
parents:
19697
diff
changeset
|
36 dim_vector::nil_rep (void) |
f3f8e1d3e399
avoid mulitple definitions of static function-scope vars (bug #47372)
John W. Eaton <jwe@octave.org>
parents:
19697
diff
changeset
|
37 { |
f3f8e1d3e399
avoid mulitple definitions of static function-scope vars (bug #47372)
John W. Eaton <jwe@octave.org>
parents:
19697
diff
changeset
|
38 static dim_vector zv (0, 0); |
f3f8e1d3e399
avoid mulitple definitions of static function-scope vars (bug #47372)
John W. Eaton <jwe@octave.org>
parents:
19697
diff
changeset
|
39 return zv.rep; |
f3f8e1d3e399
avoid mulitple definitions of static function-scope vars (bug #47372)
John W. Eaton <jwe@octave.org>
parents:
19697
diff
changeset
|
40 } |
f3f8e1d3e399
avoid mulitple definitions of static function-scope vars (bug #47372)
John W. Eaton <jwe@octave.org>
parents:
19697
diff
changeset
|
41 |
21751
b571fc85953f
maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents:
21724
diff
changeset
|
42 // The maximum allowed value for a dimension extent. This will normally be a |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
43 // tiny bit off the maximum value of octave_idx_type. |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
44 // Currently 1 is subtracted to allow safe conversion of any 2D Array into |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
45 // Sparse, but this offset may change in the future. |
10830
b4ebfd675321
avoid static initialization disaster in dim_vector
Jaroslav Hajek <highegg@gmail.com>
parents:
10810
diff
changeset
|
46 octave_idx_type |
b4ebfd675321
avoid static initialization disaster in dim_vector
Jaroslav Hajek <highegg@gmail.com>
parents:
10810
diff
changeset
|
47 dim_vector::dim_max (void) |
b4ebfd675321
avoid static initialization disaster in dim_vector
Jaroslav Hajek <highegg@gmail.com>
parents:
10810
diff
changeset
|
48 { |
b4ebfd675321
avoid static initialization disaster in dim_vector
Jaroslav Hajek <highegg@gmail.com>
parents:
10810
diff
changeset
|
49 return std::numeric_limits<octave_idx_type>::max () - 1; |
b4ebfd675321
avoid static initialization disaster in dim_vector
Jaroslav Hajek <highegg@gmail.com>
parents:
10810
diff
changeset
|
50 } |
10810
6683f0c9d742
make the maximum extent externally accessible
Jaroslav Hajek <highegg@gmail.com>
parents:
10715
diff
changeset
|
51 |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
52 void |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
53 dim_vector::chop_all_singletons (void) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
54 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
55 make_unique (); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
56 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
57 int j = 0; |
14846
460a3c6d8bf1
maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists.
Rik <octave@nomad.inbox5.com>
parents:
14624
diff
changeset
|
58 int l = ndims (); |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
59 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
60 for (int i = 0; i < l; i++) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
61 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
62 if (rep[i] != 1) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
63 rep[j++] = rep[i]; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
64 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
65 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
66 if (j == 1) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
67 rep[1] = 1; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
68 |
22212
2b487dbe419d
dim-vector.h (ndims): do not return a reference to the rep.
Carnë Draug <carandraug@octave.org>
parents:
21751
diff
changeset
|
69 rep[-1] = j > 2 ? j : 2; |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
70 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
71 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
72 std::string |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
73 dim_vector::str (char sep) const |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
74 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
75 std::ostringstream buf; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
76 |
21660
53728df3e4c9
maint: for clarity, call ndims() rather than length() on dimension vectors.
Rik <rik@octave.org>
parents:
21574
diff
changeset
|
77 for (int i = 0; i < ndims (); i++) |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
78 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
79 buf << elem (i); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
80 |
21660
53728df3e4c9
maint: for clarity, call ndims() rather than length() on dimension vectors.
Rik <rik@octave.org>
parents:
21574
diff
changeset
|
81 if (i < ndims () - 1) |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
82 buf << sep; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
83 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
84 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
85 std::string retval = buf.str (); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
86 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
87 return retval; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
88 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
89 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
90 int |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
91 dim_vector::num_ones (void) const |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
92 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
93 int retval = 0; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
94 |
21660
53728df3e4c9
maint: for clarity, call ndims() rather than length() on dimension vectors.
Rik <rik@octave.org>
parents:
21574
diff
changeset
|
95 for (int i = 0; i < ndims (); i++) |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
96 if (elem (i) == 1) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
97 retval++; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
98 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
99 return retval; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
100 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
101 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
102 octave_idx_type |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
103 dim_vector::safe_numel (void) const |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
104 { |
10830
b4ebfd675321
avoid static initialization disaster in dim_vector
Jaroslav Hajek <highegg@gmail.com>
parents:
10810
diff
changeset
|
105 octave_idx_type idx_max = dim_max (); |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
106 octave_idx_type n = 1; |
21660
53728df3e4c9
maint: for clarity, call ndims() rather than length() on dimension vectors.
Rik <rik@octave.org>
parents:
21574
diff
changeset
|
107 int n_dims = ndims (); |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
108 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
109 for (int i = 0; i < n_dims; i++) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
110 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
111 n *= rep[i]; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
112 if (rep[i] != 0) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
113 idx_max /= rep[i]; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
114 if (idx_max <= 0) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
115 throw std::bad_alloc (); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
116 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
117 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
118 return n; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
119 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
120 |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
121 dim_vector |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
122 dim_vector::squeeze (void) const |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
123 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
124 dim_vector new_dims = *this; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
125 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
126 bool dims_changed = 1; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
127 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
128 int k = 0; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
129 |
21660
53728df3e4c9
maint: for clarity, call ndims() rather than length() on dimension vectors.
Rik <rik@octave.org>
parents:
21574
diff
changeset
|
130 for (int i = 0; i < ndims (); i++) |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
131 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
132 if (elem (i) == 1) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
133 dims_changed = true; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
134 else |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
135 new_dims(k++) = elem (i); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
136 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
137 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
138 if (dims_changed) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
139 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
140 if (k == 0) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
141 new_dims = dim_vector (1, 1); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
142 else if (k == 1) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
143 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
144 // There is one non-singleton dimension, so we need |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
145 // to decide the correct orientation. |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
146 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
147 if (elem (0) == 1) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
148 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
149 // The original dimension vector had a leading |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
150 // singleton dimension. |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
151 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
152 octave_idx_type tmp = new_dims(0); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
153 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
154 new_dims.resize (2); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
155 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
156 new_dims(0) = 1; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
157 new_dims(1) = tmp; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
158 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
159 else |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
160 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
161 // The first element of the original dimension vector |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
162 // was not a singleton dimension. |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
163 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
164 new_dims.resize (2); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
165 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
166 new_dims(1) = 1; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
167 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
168 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
169 else |
15018
3d8ace26c5b4
maint: Use Octave coding conventions for cuddled parentheses in liboctave/.
Rik <rik@octave.org>
parents:
14846
diff
changeset
|
170 new_dims.resize (k); |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
171 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
172 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
173 return new_dims; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
174 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
175 |
21751
b571fc85953f
maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents:
21724
diff
changeset
|
176 // This is the rule for cat(). cat (dim, A, B) works if one |
10715
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
177 // of the following holds, in this order: |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
178 // |
15018
3d8ace26c5b4
maint: Use Octave coding conventions for cuddled parentheses in liboctave/.
Rik <rik@octave.org>
parents:
14846
diff
changeset
|
179 // 1. size (A, k) == size (B, k) for all k != dim. |
10715
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
180 // In this case, size (C, dim) = size (A, dim) + size (B, dim) and |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
181 // other sizes remain intact. |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
182 // |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
183 // 2. A is 0x0, in which case B is the result |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
184 // 3. B is 0x0, in which case A is the result |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
185 |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
186 bool |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
187 dim_vector::concat (const dim_vector& dvb, int dim) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
188 { |
18084
8e056300994b
Follow coding convention of defining and initializing only 1 variable per line in liboctave.
Rik <rik@octave.org>
parents:
17769
diff
changeset
|
189 int orig_nd = ndims (); |
8e056300994b
Follow coding convention of defining and initializing only 1 variable per line in liboctave.
Rik <rik@octave.org>
parents:
17769
diff
changeset
|
190 int ndb = dvb.ndims (); |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
191 int new_nd = dim < ndb ? ndb : dim + 1; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
192 if (new_nd > orig_nd) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
193 resize (new_nd, 1); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
194 else |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
195 new_nd = orig_nd; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
196 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
197 make_unique (); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
198 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
199 bool match = true; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
200 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
201 for (int i = 0; i < ndb; i++) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
202 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
203 if (i != dim && rep[i] != dvb(i)) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
204 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
205 match = false; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
206 break; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
207 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
208 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
209 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
210 for (int i = ndb; i < new_nd; i++) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
211 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
212 if (i != dim && rep[i] != 1) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
213 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
214 match = false; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
215 break; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
216 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
217 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
218 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
219 if (match) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
220 rep[dim] += (dim < ndb ? dvb(dim) : 1); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
221 else |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
222 { |
21751
b571fc85953f
maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents:
21724
diff
changeset
|
223 // Dimensions don't match. The only allowed fix is to omit 0x0. |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
224 if (ndb == 2 && dvb(0) == 0 && dvb(1) == 0) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
225 match = true; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
226 else if (orig_nd == 2 && rep[0] == 0 && rep[1] == 0) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
227 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
228 *this = dvb; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
229 match = true; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
230 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
231 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
232 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
233 chop_trailing_singletons (); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
234 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
235 return match; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
236 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
237 |
10715
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
238 // Rules for horzcat/vertcat are yet looser. |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
239 // two arrays A, B can be concatenated |
10715
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
240 // horizontally (dim = 2) or vertically (dim = 1) if one of the |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
241 // following holds, in this order: |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
242 // |
15018
3d8ace26c5b4
maint: Use Octave coding conventions for cuddled parentheses in liboctave/.
Rik <rik@octave.org>
parents:
14846
diff
changeset
|
243 // 1. cat (dim, A, B) works |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
244 // |
10715
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
245 // 2. A, B are 2D and one of them is an empty vector, in which |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
246 // case the result is the other one except if both of them |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
247 // are empty vectors, in which case the result is 0x0. |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
248 |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
249 bool |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
250 dim_vector::hvcat (const dim_vector& dvb, int dim) |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
251 { |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
252 if (concat (dvb, dim)) |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
253 return true; |
21660
53728df3e4c9
maint: for clarity, call ndims() rather than length() on dimension vectors.
Rik <rik@octave.org>
parents:
21574
diff
changeset
|
254 else if (ndims () == 2 && dvb.ndims () == 2) |
10715
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
255 { |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
256 bool e2dv = rep[0] + rep[1] == 1; |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
257 bool e2dvb = dvb(0) + dvb(1) == 1; |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
258 if (e2dvb) |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
259 { |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
260 if (e2dv) |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
261 *this = dim_vector (); |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
262 return true; |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
263 } |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
264 else if (e2dv) |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
265 { |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
266 *this = dvb; |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
267 return true; |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
268 } |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
269 } |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
270 |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
271 return false; |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
272 } |
53253f796351
make [] (hopefully) more Matlab compatible
Jaroslav Hajek <highegg@gmail.com>
parents:
10686
diff
changeset
|
273 |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
274 dim_vector |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
275 dim_vector::redim (int n) const |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
276 { |
21660
53728df3e4c9
maint: for clarity, call ndims() rather than length() on dimension vectors.
Rik <rik@octave.org>
parents:
21574
diff
changeset
|
277 int n_dims = ndims (); |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
278 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
279 if (n_dims == n) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
280 return *this; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
281 else if (n_dims < n) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
282 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
283 dim_vector retval = alloc (n); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
284 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
285 for (int i = 0; i < n_dims; i++) |
14624
edf9ca8a92a8
when redimensioning, always pad dim_vector objects with 1 (bug #33216)
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
286 retval.rep[i] = rep[i]; |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
287 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
288 for (int i = n_dims; i < n; i++) |
14624
edf9ca8a92a8
when redimensioning, always pad dim_vector objects with 1 (bug #33216)
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
289 retval.rep[i] = 1; |
10683
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
290 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
291 return retval; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
292 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
293 else |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
294 { |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
295 if (n < 1) n = 1; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
296 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
297 dim_vector retval = alloc (n); |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
298 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
299 retval.rep[1] = 1; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
300 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
301 for (int i = 0; i < n-1; i++) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
302 retval.rep[i] = rep[i]; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
303 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
304 int k = rep[n-1]; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
305 for (int i = n; i < n_dims; i++) |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
306 k *= rep[i]; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
307 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
308 retval.rep[n-1] = k; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
309 |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
310 return retval; |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
311 } |
4fdb36ffa066
omitted file from 10681:0ba9bd294421
Jaroslav Hajek <highegg@gmail.com>
parents:
diff
changeset
|
312 } |