octave-lyh: liboctave/dim-vector.h annotate

annotate liboctave/dim-vector.h @ 10419:afe44ee90cbd

implement generic macro magic for repeated decls

author	Jaroslav Hajek <highegg@gmail.com>
date	Wed, 17 Mar 2010 10:30:36 +0100
parents	6c19d6fcd7e5
children	3373fdc0b14a

rev	line source
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	1 /*
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	2
8920 eb63fbe60fab update copyright notices John W. Eaton <jwe@octave.org> parents: 8721 diff changeset	3 Copyright (C) 2003, 2004, 2005, 2006, 2007, 2008, 2009 John W. Eaton
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	4 Copyirght (C) 2009 VZLU Prague
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	5
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	6 This file is part of Octave.
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	7
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	8 Octave is free software; you can redistribute it and/or modify it
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	9 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6216 diff changeset	10 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6216 diff changeset	11 option) any later version.
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	12
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	13 Octave is distributed in the hope that it will be useful, but WITHOUT
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	16 for more details.
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	17
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	18 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6216 diff changeset	19 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6216 diff changeset	20 <http://www.gnu.org/licenses/>.
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	21
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	22 */
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	23
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	24 #if !defined (octave_dim_vector_h)
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	25 #define octave_dim_vector_h 1
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	26
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	27 #include <cassert>
9840 c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	28 #include <limits>
8950 d865363208d6 include <iosfwd> instead of <iostream> in header files John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	29
5765 7ba9ad1fec11 [project @ 2006-04-17 05:05:15 by jwe] jwe parents: 5307 diff changeset	30 #include <sstream>
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	31 #include <string>
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	32
6216 74ca56e4414c [project @ 2006-12-22 17:45:47 by jwe] jwe parents: 5765 diff changeset	33 #include "lo-error.h"
10419 afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	34 #include "lo-macros.h"
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	35
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	36 // Rationale: This implementation is more tricky than Array, but the
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	37 // big plus is that dim_vector requires only one allocation instead of
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	38 // two. It is (slightly) patterned after GCC's basic_string
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	39 // implementation. rep is a pointer to an array of memory, comprising
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	40 // count, length, and the data:
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	41 //
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	42 // <count>
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	43 // <ndims>
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	44 // rep --> <dims[0]>
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	45 // <dims[1]>
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	46 // ...
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	47 //
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	48 // The inlines count(), ndims() recover this data from the rep. Note
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	49 // that rep points to the beginning of dims to grant faster access
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	50 // (reinterpret_cast is assumed to be an inexpensive operation).
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	51
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	52 class
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	53 dim_vector
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	54 {
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	55 private:
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	56
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	57 octave_idx_type *rep;
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	58
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	59 octave_idx_type& ndims (void) const { return rep[-1]; }
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	60
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	61 octave_idx_type& count (void) const { return rep[-2]; }
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	62
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	63 // Construct a new rep with count = 1 and ndims given.
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	64
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	65 static octave_idx_type *newrep (int ndims)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	66 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	67 octave_idx_type *r = new octave_idx_type[ndims + 2];
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	68
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	69 *r++ = 1;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	70 *r++ = ndims;
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	71
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	72 return r;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	73 }
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	74
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	75 // Clone this->rep.
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	76
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	77 octave_idx_type *clonerep (void)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	78 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	79 int l = ndims ();
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	80
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	81 octave_idx_type *r = new octave_idx_type[l + 2];
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	82
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	83 *r++ = 1;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	84 *r++ = l;
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	85
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	86 for (int i = 0; i < l; i++)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	87 r[i] = rep[i];
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	88
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	89 return r;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	90 }
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	91
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	92 // Clone and resize this->rep to length n, filling by given value.
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	93
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	94 octave_idx_type *resizerep (int n, octave_idx_type fill_value)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	95 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	96 int l = ndims ();
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	97
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	98 if (n < 2)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	99 n = 2;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	100
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	101 octave_idx_type *r = new octave_idx_type[n + 2];
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	102
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	103 *r++ = 1;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	104 *r++ = n;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	105
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	106 if (l > n)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	107 l = n;
4673 6639c2574023 [project @ 2003-12-10 16:11:54 by jwe] jwe parents: 4655 diff changeset	108
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	109 int j;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	110 for (j = 0; j < l; j++)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	111 r[j] = rep[j];
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	112 for (; j < n; j++)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	113 r[j] = fill_value;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	114
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	115 return r;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	116 }
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	117
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	118 // Free the rep.
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	119
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	120 void freerep (void)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	121 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	122 assert (count () == 0);
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	123 delete [] (rep - 2);
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	124 }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	125
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	126 void make_unique (void)
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	127 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	128 if (count () > 1)
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	129 {
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	130 --count();
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	131 rep = clonerep ();
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	132 }
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	133 }
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	134
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	135 public:
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	136
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	137 explicit dim_vector (octave_idx_type n) GCC_ATTR_DEPRECATED
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	138 : rep (newrep (2))
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	139 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	140 rep[0] = n;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	141 rep[1] = 1;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	142 }
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	143
10418 6c19d6fcd7e5 up to 7th-order dim_vector constructors Jaroslav Hajek <highegg@gmail.com> parents: 10403 diff changeset	144 dim_vector (octave_idx_type r, octave_idx_type c)
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	145 : rep (newrep (2))
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	146 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	147 rep[0] = r;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	148 rep[1] = c;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	149 }
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	150
10418 6c19d6fcd7e5 up to 7th-order dim_vector constructors Jaroslav Hajek <highegg@gmail.com> parents: 10403 diff changeset	151 dim_vector (octave_idx_type r, octave_idx_type c, octave_idx_type p)
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	152 : rep (newrep (3))
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	153 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	154 rep[0] = r;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	155 rep[1] = c;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	156 rep[2] = p;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	157 }
10419 afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	158
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	159 #define ASSIGN_REP(i) rep[i] = d ## i;
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	160 #define DIM_VECTOR_CTOR(N) \
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	161 dim_vector (OCT_MAKE_DECL_LIST(octave_idx_type, d, N)) \
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	162 : rep (newrep (N)) \
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	163 { \
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	164 OCT_ITERATE_MACRO(ASSIGN_REP, N) \
10418 6c19d6fcd7e5 up to 7th-order dim_vector constructors Jaroslav Hajek <highegg@gmail.com> parents: 10403 diff changeset	165 }
10419 afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	166
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	167 // Add more if needed.
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	168 DIM_VECTOR_CTOR(4)
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	169 DIM_VECTOR_CTOR(5)
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	170 DIM_VECTOR_CTOR(6)
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	171 DIM_VECTOR_CTOR(7)
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	172
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	173 #undef ASSIGN_REP
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	174 #undef DIM_VECTOR_CTOR
afe44ee90cbd implement generic macro magic for repeated decls Jaroslav Hajek <highegg@gmail.com> parents: 10418 diff changeset	175
10403 69ecfbffcf4f liboctave/dim-vector.h: New constructor accepting a C array of dimensions. Soren Hauberg <hauberg@gmail.com> parents: 10366 diff changeset	176 dim_vector (const octave_idx_type *vec, size_t vec_size)
69ecfbffcf4f liboctave/dim-vector.h: New constructor accepting a C array of dimensions. Soren Hauberg <hauberg@gmail.com> parents: 10366 diff changeset	177 : rep (newrep (vec_size))
69ecfbffcf4f liboctave/dim-vector.h: New constructor accepting a C array of dimensions. Soren Hauberg <hauberg@gmail.com> parents: 10366 diff changeset	178 {
69ecfbffcf4f liboctave/dim-vector.h: New constructor accepting a C array of dimensions. Soren Hauberg <hauberg@gmail.com> parents: 10366 diff changeset	179 for (size_t k = 0; k < vec_size; k++)
69ecfbffcf4f liboctave/dim-vector.h: New constructor accepting a C array of dimensions. Soren Hauberg <hauberg@gmail.com> parents: 10366 diff changeset	180 rep[k] = vec[k];
69ecfbffcf4f liboctave/dim-vector.h: New constructor accepting a C array of dimensions. Soren Hauberg <hauberg@gmail.com> parents: 10366 diff changeset	181 }
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	182
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	183 octave_idx_type& elem (int i)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	184 {
10366 e5ae13b8b2c2 improve Array indexing error messages Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	185 #ifdef BOUNDS_CHECKING
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	186 assert (i >= 0 && i < ndims ());
10366 e5ae13b8b2c2 improve Array indexing error messages Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	187 #endif
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	188 make_unique ();
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	189 return rep[i];
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	190 }
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	191
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	192 octave_idx_type elem (int i) const
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	193 {
10366 e5ae13b8b2c2 improve Array indexing error messages Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	194 #ifdef BOUNDS_CHECKING
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	195 assert (i >= 0 && i < ndims ());
10366 e5ae13b8b2c2 improve Array indexing error messages Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	196 #endif
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	197 return rep[i];
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	198 }
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	199
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	200 void chop_trailing_singletons (void)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	201 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	202 int l = ndims ();
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	203 if (l > 2 && rep[l-1] == 1)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	204 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	205 make_unique ();
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	206 do
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	207 l--;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	208 while (l > 2 && rep[l-1] == 1);
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	209 ndims () = l;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	210 }
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	211 }
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	212
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	213 void chop_all_singletons (void)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	214 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	215 make_unique ();
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	216
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	217 int j = 0;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	218 int l = ndims();
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	219
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	220 for (int i = 0; i < l; i++)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	221 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	222 if (rep[i] != 1)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	223 rep[j++] = rep[i];
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	224 }
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	225
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	226 if (j == 1)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	227 rep[1] = 1;
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	228
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	229 ndims () = j > 2 ? j : 2;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	230 }
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	231
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	232 private:
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	233
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	234 static octave_idx_type *nil_rep (void)
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	235 {
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	236 static dim_vector zv (0, 0);
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	237 return zv.rep;
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	238 }
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	239
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	240 explicit dim_vector (octave_idx_type *r)
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	241 : rep (r) { }
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	242
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	243 public:
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	244
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	245 explicit dim_vector (void) : rep (nil_rep ()) { count()++; }
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	246
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	247 dim_vector (const dim_vector& dv) : rep (dv.rep) { count()++; }
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	248
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	249 static dim_vector alloc (int n)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	250 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	251 return dim_vector (newrep (n < 2 ? 2 : n));
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	252 }
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	253
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	254 dim_vector& operator = (const dim_vector& dv)
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	255 {
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	256 if (&dv != this)
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	257 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	258 if (--count() <= 0)
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	259 freerep ();
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	260
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	261 rep = dv.rep;
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	262 count()++;
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	263 }
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	264
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	265 return *this;
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	266 }
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	267
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	268 ~dim_vector (void)
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	269 {
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	270 if (--count() <= 0)
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	271 freerep ();
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	272 }
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	273
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	274 int length (void) const { return ndims (); }
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	275
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	276 octave_idx_type& operator () (int i) { return elem (i); }
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	277
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	278 octave_idx_type operator () (int i) const { return elem (i); }
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	279
4887 bbddd4339cf2 [project @ 2004-04-30 17:08:55 by jwe] jwe parents: 4735 diff changeset	280 void resize (int n, int fill_value = 0)
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	281 {
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	282 int len = length ();
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	283
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	284 if (n != len)
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	285 {
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	286 octave_idx_type *r = resizerep (n, fill_value);
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	287
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	288 if (--count() <= 0)
b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	289 freerep ();
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	290
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	291 rep = r;
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	292 }
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	293 }
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	294
4559 802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	295 std::string str (char sep = 'x') const
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	296 {
5765 7ba9ad1fec11 [project @ 2006-04-17 05:05:15 by jwe] jwe parents: 5307 diff changeset	297 std::ostringstream buf;
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	298
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	299 for (int i = 0; i < length (); i++)
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	300 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	301 buf << elem (i);
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	302
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	303 if (i < length () - 1)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	304 buf << sep;
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	305 }
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	306
5765 7ba9ad1fec11 [project @ 2006-04-17 05:05:15 by jwe] jwe parents: 5307 diff changeset	307 std::string retval = buf.str ();
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	308
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	309 return retval;
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	310 }
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	311
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	312 bool all_zero (void) const
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	313 {
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	314 bool retval = true;
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	315
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	316 for (int i = 0; i < length (); i++)
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	317 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	318 if (elem (i) != 0)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	319 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	320 retval = false;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	321 break;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	322 }
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	323 }
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	324
4548 820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	325 return retval;
820323598f4f [project @ 2003-10-26 03:16:42 by jwe] jwe parents: 4543 diff changeset	326 }
4559 802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	327
9886 cddd5c3d5f04 fix & extend special-case optimizations for indexed assignment Jaroslav Hajek <highegg@gmail.com> parents: 9840 diff changeset	328 bool zero_by_zero (void) const
cddd5c3d5f04 fix & extend special-case optimizations for indexed assignment Jaroslav Hajek <highegg@gmail.com> parents: 9840 diff changeset	329 {
cddd5c3d5f04 fix & extend special-case optimizations for indexed assignment Jaroslav Hajek <highegg@gmail.com> parents: 9840 diff changeset	330 return length () == 2 && elem (0) == 0 && elem (1) == 0;
cddd5c3d5f04 fix & extend special-case optimizations for indexed assignment Jaroslav Hajek <highegg@gmail.com> parents: 9840 diff changeset	331 }
cddd5c3d5f04 fix & extend special-case optimizations for indexed assignment Jaroslav Hajek <highegg@gmail.com> parents: 9840 diff changeset	332
4559 802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	333 bool any_zero (void) const
802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	334 {
802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	335 bool retval = false;
802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	336
802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	337 for (int i = 0; i < length (); i++)
802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	338 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	339 if (elem (i) == 0)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	340 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	341 retval = true;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	342 break;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	343 }
4559 802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	344 }
802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	345
802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	346 return retval;
802818bfac91 [project @ 2003-10-28 19:13:43 by jwe] jwe parents: 4548 diff changeset	347 }
4567 fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	348
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	349 int num_ones (void) const
4635 f2cd320cbf6e [project @ 2003-11-20 01:48:23 by jwe] jwe parents: 4567 diff changeset	350 {
f2cd320cbf6e [project @ 2003-11-20 01:48:23 by jwe] jwe parents: 4567 diff changeset	351 int retval = 0;
f2cd320cbf6e [project @ 2003-11-20 01:48:23 by jwe] jwe parents: 4567 diff changeset	352
f2cd320cbf6e [project @ 2003-11-20 01:48:23 by jwe] jwe parents: 4567 diff changeset	353 for (int i = 0; i < length (); i++)
f2cd320cbf6e [project @ 2003-11-20 01:48:23 by jwe] jwe parents: 4567 diff changeset	354 if (elem (i) == 1)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	355 retval++;
4635 f2cd320cbf6e [project @ 2003-11-20 01:48:23 by jwe] jwe parents: 4567 diff changeset	356
f2cd320cbf6e [project @ 2003-11-20 01:48:23 by jwe] jwe parents: 4567 diff changeset	357 return retval;
f2cd320cbf6e [project @ 2003-11-20 01:48:23 by jwe] jwe parents: 4567 diff changeset	358 }
f2cd320cbf6e [project @ 2003-11-20 01:48:23 by jwe] jwe parents: 4567 diff changeset	359
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	360 bool all_ones (void) const
4655 c8829691db47 [project @ 2003-11-24 21:24:37 by jwe] jwe parents: 4635 diff changeset	361 {
c8829691db47 [project @ 2003-11-24 21:24:37 by jwe] jwe parents: 4635 diff changeset	362 return (num_ones () == length ());
c8829691db47 [project @ 2003-11-24 21:24:37 by jwe] jwe parents: 4635 diff changeset	363 }
c8829691db47 [project @ 2003-11-24 21:24:37 by jwe] jwe parents: 4635 diff changeset	364
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	365 // Return the number of elements that a matrix with this dimension
4567 fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	366 // vector would have, NOT the number of dimensions (elements in the
fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	367 // dimension vector).
fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	368
9027 9a46ba093db4 generalize dim_vector::numel Jaroslav Hajek <highegg@gmail.com> parents: 8950 diff changeset	369 octave_idx_type numel (int n = 0) const
4567 fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	370 {
fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	371 int n_dims = length ();
fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	372
9027 9a46ba093db4 generalize dim_vector::numel Jaroslav Hajek <highegg@gmail.com> parents: 8950 diff changeset	373 octave_idx_type retval = 1;
4567 fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	374
9027 9a46ba093db4 generalize dim_vector::numel Jaroslav Hajek <highegg@gmail.com> parents: 8950 diff changeset	375 for (int i = n; i < n_dims; i++)
4567 fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	376 retval *= elem (i);
fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	377
fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	378 return retval;
fc30061d01da [project @ 2003-10-29 20:11:15 by jwe] jwe parents: 4565 diff changeset	379 }
4673 6639c2574023 [project @ 2003-12-10 16:11:54 by jwe] jwe parents: 4655 diff changeset	380
9840 c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	381 // The following function will throw a std::bad_alloc ()
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	382 // exception if the requested size is larger than can be indexed by
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	383 // octave_idx_type. This may be smaller than the actual amount of
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	384 // memory that can be safely allocated on a system. However, if we
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	385 // don't fail here, we can end up with a mysterious crash inside a
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	386 // function that is iterating over an array using octave_idx_type
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	387 // indices.
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	388
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	389 octave_idx_type safe_numel (void) const
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	390 {
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	391 octave_idx_type idx_max = std::numeric_limits<octave_idx_type>::max () - 1;
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	392 octave_idx_type n = 1;
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	393 int n_dims = length ();
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	394
9840 c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	395 for (int i = 0; i < n_dims; i++)
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	396 {
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	397 n *= rep[i];
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	398 if (rep[i] != 0)
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	399 idx_max /= rep[i];
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	400 if (idx_max <= 0)
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	401 throw std::bad_alloc ();
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	402 }
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	403
9840 c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	404 return n;
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	405 }
c0b54271904b improve safe numel() calculation for arrays Jaroslav Hajek <highegg@gmail.com> parents: 9743 diff changeset	406
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	407 bool any_neg (void) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	408 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	409 int n_dims = length ();
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	410 int i;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	411
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	412 for (i = 0; i < n_dims; i++)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	413 if (elem (i) < 0)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	414 break;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	415
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	416 return i < n_dims;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	417 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	418
4735 24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	419 dim_vector squeeze (void) const
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	420 {
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	421 dim_vector new_dims = *this;
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	422
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	423 bool dims_changed = 1;
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	424
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	425 int k = 0;
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	426
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	427 for (int i = 0; i < length (); i++)
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	428 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	429 if (elem (i) == 1)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	430 dims_changed = true;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	431 else
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	432 new_dims(k++) = elem (i);
4735 24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	433 }
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	434
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	435 if (dims_changed)
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	436 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	437 if (k == 0)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	438 new_dims = dim_vector (1, 1);
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	439 else if (k == 1)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	440 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	441 // There is one non-singleton dimension, so we need
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	442 // to decide the correct orientation.
4735 24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	443
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	444 if (elem (0) == 1)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	445 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	446 // The original dimension vector had a leading
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	447 // singleton dimension.
4735 24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	448
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	449 octave_idx_type tmp = new_dims(0);
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	450
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	451 new_dims.resize (2);
4735 24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	452
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	453 new_dims(0) = 1;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	454 new_dims(1) = tmp;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	455 }
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	456 else
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	457 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	458 // The first element of the original dimension vector
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	459 // was not a singleton dimension.
4735 24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	460
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	461 new_dims.resize (2);
4735 24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	462
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	463 new_dims(1) = 1;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	464 }
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	465 }
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	466 else
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	467 new_dims.resize(k);
4735 24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	468 }
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	469
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	470 return new_dims;
24c7bc6354ba [project @ 2004-02-05 17:26:28 by jwe] jwe parents: 4673 diff changeset	471 }
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	472
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	473 bool concat (const dim_vector& dvb, int dim = 0)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	474 {
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	475 if (all_zero ())
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	476 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	477 *this = dvb;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	478 return true;
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	479 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	480
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	481 if (dvb.all_zero ())
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	482 return true;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	483
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	484 int na = length ();
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	485 int nb = dvb.length ();
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	486
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	487 // Find the max and min value of na and nb
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	488 int n_max = na > nb ? na : nb;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	489 int n_min = na < nb ? na : nb;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	490
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	491 // The elements of the dimension vectors can only differ
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	492 // if the dim variable differs from the actual dimension
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	493 // they differ.
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	494
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	495 for (int i = 0; i < n_min; i++)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	496 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	497 if (elem(i) != dvb(i) && dim != i)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	498 return false;
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	499 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	500
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	501 // Ditto.
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	502 for (int i = n_min; i < n_max; i++)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	503 {
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	504 if (na > n_min)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	505 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	506 if (elem(i) != 1 && dim != i)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	507 return false;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	508 }
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	509 else
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	510 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	511 if (dvb(i) != 1 && dim != i)
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	512 return false;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	513 }
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	514 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	515
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	516 // If we want to add the dimension vectors at a dimension
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	517 // larger than both, then we need to set n_max to this number
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	518 // so that we resize *this to the right dimension.
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	519
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	520 n_max = n_max > (dim + 1) ? n_max : (dim + 1);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	521
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	522 // Resize *this to the appropriate dimensions.
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	523
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	524 if (n_max > na)
9507 b096d11237be dim_vector improvements Jaroslav Hajek <highegg@gmail.com> parents: 9420 diff changeset	525 resize (n_max, 1);
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	526
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	527 // Larger or equal since dim has been decremented by one.
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	528
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	529 if (dim >= nb)
4940 954cc2ba6a49 [project @ 2004-08-24 19:30:31 by jwe] jwe parents: 4915 diff changeset	530 elem (dim)++;
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	531 else
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	532 elem (dim) += dvb(dim);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	533
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	534 return true;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4887 diff changeset	535 }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	536
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	537 // Force certain dimensionality, preserving numel (). Missing
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	538 // dimensions are set to 1, redundant are folded into the trailing
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	539 // one. If n = 1, the result is 2d and the second dim is 1
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	540 // (dim_vectors are always at least 2D). If the original dimensions
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	541 // were all zero, the padding value is zero.
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	542
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	543 dim_vector redim (int n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	544 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	545 int n_dims = length ();
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	546
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	547 if (n_dims == n)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	548 return *this;
9666 a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	549 else if (n_dims < n)
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	550 {
9666 a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	551 dim_vector retval = alloc (n);
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	552
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	553 int pad = 0;
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	554 for (int i = 0; i < n_dims; i++)
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	555 {
9666 a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	556 retval.rep[i] = rep[i];
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	557 if (rep[i] != 0)
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	558 pad = 1;
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	559 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	560
9666 a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	561 for (int i = n_dims; i < n; i++)
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	562 retval.rep[i] = pad;
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	563
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	564 return retval;
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	565 }
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	566 else
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	567 {
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	568 if (n < 1) n = 1;
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	569
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	570 dim_vector retval = alloc (n);
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	571
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	572 retval.rep[1] = 1;
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	573
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	574 for (int i = 0; i < n-1; i++)
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	575 retval.rep[i] = rep[i];
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	576
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	577 int k = rep[n-1];
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	578 for (int i = n; i < n_dims; i++)
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	579 k *= rep[i];
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	580
a531dec450c4 allow 1D case for sub2ind and ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 9648 diff changeset	581 retval.rep[n-1] = k;
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	582
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	583 return retval;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	584 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	585 }
8721 e9cb742df9eb imported patch sort3.diff Jaroslav Hajek <highegg@gmail.com> parents: 8290 diff changeset	586
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	587 dim_vector as_column (void) const
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	588 {
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	589 if (length () == 2 && elem (1) == 1)
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	590 return *this;
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	591 else
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	592 return dim_vector (numel (), 1);
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	593 }
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	594
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	595 dim_vector as_row (void) const
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	596 {
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	597 if (length () == 2 && elem (0) == 1)
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	598 return *this;
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	599 else
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	600 return dim_vector (1, numel ());
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	601 }
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10150 diff changeset	602
8721 e9cb742df9eb imported patch sort3.diff Jaroslav Hajek <highegg@gmail.com> parents: 8290 diff changeset	603 bool is_vector (void) const
e9cb742df9eb imported patch sort3.diff Jaroslav Hajek <highegg@gmail.com> parents: 8290 diff changeset	604 {
e9cb742df9eb imported patch sort3.diff Jaroslav Hajek <highegg@gmail.com> parents: 8290 diff changeset	605 return (length () == 2 && (elem (0) == 1 \|\| elem (1) == 1));
e9cb742df9eb imported patch sort3.diff Jaroslav Hajek <highegg@gmail.com> parents: 8290 diff changeset	606 }
e9cb742df9eb imported patch sort3.diff Jaroslav Hajek <highegg@gmail.com> parents: 8290 diff changeset	607
9511 cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	608 int first_non_singleton (int def = 0) const
cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	609 {
cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	610 for (int i = 0; i < length (); i++)
cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	611 {
cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	612 if (elem (i) != 1)
cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	613 return i;
cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	614 }
cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	615
cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	616 return def;
cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	617 }
cc1fd3084cb2 implement dim_vector::first_non_singleton Jaroslav Hajek <highegg@gmail.com> parents: 9507 diff changeset	618
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	619 // Compute a linear index from an index tuple.
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	620
9739 13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	621 octave_idx_type compute_index (const octave_idx_type *idx)
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	622 {
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	623 octave_idx_type k = 0;
9743 26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9739 diff changeset	624 for (int i = length () - 1; i >= 0; i--)
9739 13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	625 k = k * rep[i] + idx[i];
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	626
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	627 return k;
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	628 }
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	629
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	630 // Increment a multi-dimensional index tuple, optionally starting
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	631 // from an offset position and return the index of the last index
9739 13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	632 // position that was changed, or length () if just cycled over.
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	633
9739 13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	634 int increment_index (octave_idx_type *idx, int start = 0)
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	635 {
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	636 int i;
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	637 for (i = start; i < length (); i++)
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	638 {
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	639 if (++(*idx) == rep[i])
9743 26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9739 diff changeset	640 *idx++ = 0;
9739 13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	641 else
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	642 break;
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	643 }
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	644 return i;
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	645 }
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	646
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	647 // Return cumulative dimensions.
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	648
9739 13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	649 dim_vector cumulative (void) const
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	650 {
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	651 int nd = length ();
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	652 dim_vector retval = alloc (nd);
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	653
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	654 octave_idx_type k = 1;
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	655 for (int i = 0; i < nd; i++)
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	656 retval.rep[i] = k *= rep[i];
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	657
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	658 return retval;
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	659 }
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	660
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	661 // Compute a linear index from an index tuple. Dimensions are
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	662 // required to be cumulative.
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	663
9739 13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	664 octave_idx_type cum_compute_index (const octave_idx_type *idx)
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	665 {
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	666 octave_idx_type k = idx[0];
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	667
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	668 for (int i = 1; i < length (); i++)
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	669 k += rep[i-1] * idx[i];
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	670
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	671 return k;
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	672 }
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	673
13b57eec9440 a few handy methods for dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	674
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 9666 diff changeset	675 friend bool operator == (const dim_vector& a, const dim_vector& b);
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	676 };
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	677
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 9666 diff changeset	678 inline bool
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	679 operator == (const dim_vector& a, const dim_vector& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	680 {
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 9666 diff changeset	681 // Fast case.
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 9666 diff changeset	682 if (a.rep == b.rep)
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 9666 diff changeset	683 return true;
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 9666 diff changeset	684
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	685 bool retval = true;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	686
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	687 int a_len = a.length ();
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	688 int b_len = b.length ();
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	689
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	690 if (a_len != b_len)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	691 retval = false;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	692 else
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	693 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	694 for (int i = 0; i < a_len; i++)
10150 3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	695 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	696 if (a(i) != b(i))
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	697 {
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	698 retval = false;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	699 break;
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	700 }
3778ea95adea style fixes John W. Eaton <jwe@octave.org> parents: 10095 diff changeset	701 }
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	702 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	703
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	704 return retval;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	705 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	706
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 9666 diff changeset	707 inline bool
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	708 operator != (const dim_vector& a, const dim_vector& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	709 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	710 return ! operator == (a, b);
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	711 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	712
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: diff changeset	713 #endif

Mercurial > hg > octave-lyh

annotate liboctave/dim-vector.h @ 10419:afe44ee90cbd