octave: liboctave/idx-vector.h annotate

annotate liboctave/idx-vector.h @ 8290:7cbe01c21986

improve dense array indexing

author	Jaroslav Hajek <highegg@gmail.com>
date	Mon, 20 Oct 2008 16:54:28 +0200
parents	755bf7ecc29b
children	7e87d3d76a56

rev	line source
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	1 /*
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	2
7017 a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe] jwe parents: 7016 diff changeset	3 Copyright (C) 1993, 1994, 1995, 1996, 1997, 1998, 2000, 2002, 2003,
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe] jwe parents: 7016 diff changeset	4 2004, 2005, 2006, 2007 John W. Eaton
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	5 Copyright (C) 2008 Jaroslav Hajek
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	6
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	7 This file is part of Octave.
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	8
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	9 Octave is free software; you can redistribute it and/or modify it
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	10 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: 6108 diff changeset	11 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	12 option) any later version.
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	13
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	14 Octave is distributed in the hope that it will be useful, but WITHOUT
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	15 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	16 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	17 for more details.
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	18
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	19 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: 6108 diff changeset	20 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	21 <http://www.gnu.org/licenses/>.
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	22
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	23 */
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	24
383 0b52c68ec81f [project @ 1994-03-09 21:19:44 by jwe] jwe parents: 240 diff changeset	25 #if !defined (octave_idx_vector_h)
0b52c68ec81f [project @ 1994-03-09 21:19:44 by jwe] jwe parents: 240 diff changeset	26 #define octave_idx_vector_h 1
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	27
3503 d14c483b3c12 [project @ 2000-02-01 04:06:07 by jwe] jwe parents: 3145 diff changeset	28 #include <iostream>
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	29 #include <algorithm>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	30 #include <cassert>
3503 d14c483b3c12 [project @ 2000-02-01 04:06:07 by jwe] jwe parents: 3145 diff changeset	31
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	32 #include "dim-vector.h"
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	33 #include "oct-inttypes.h"
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	34
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	35 template<class T> class Array;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	36 class Range;
4504 f6a61399bc5c [project @ 2003-09-09 17:48:00 by jwe] jwe parents: 4461 diff changeset	37
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	38 // Design rationale:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	39 // idx_vector is a reference-counting, polymorphic pointer, that can contain
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	40 // 4 types of index objects: a magic colon, a range, a scalar, or an index vector.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	41 // Polymorphic methods for single element access are provided, as well as
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	42 // templates implementing "early dispatch", i.e. hoisting the checks for index
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	43 // type out of loops.
164 e2c950dd96d2 [project @ 1993-10-18 19:32:00 by jwe] jwe parents: 1 diff changeset	44
453 393e95f46b51 [project @ 1994-06-06 00:05:20 by jwe] jwe parents: 434 diff changeset	45 class
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5781 diff changeset	46 OCTAVE_API
453 393e95f46b51 [project @ 1994-06-06 00:05:20 by jwe] jwe parents: 434 diff changeset	47 idx_vector
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	48 {
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	49 public:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	50
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	51 enum idx_class_type
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	52 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	53 class_invalid = -1,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	54 class_colon = 0,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	55 class_range,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	56 class_scalar,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	57 class_vector
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	58 };
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	59
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	60 private:
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	61
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	62 class idx_base_rep
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	63 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	64 public:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	65 idx_base_rep (void) : count (1), err (false) { }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	66
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	67 virtual ~idx_base_rep (void) { }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	68
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	69 // Non-range-checking element query.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	70 virtual octave_idx_type xelem (octave_idx_type i) const = 0;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	71
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	72 // Range-checking element query.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	73 virtual octave_idx_type checkelem (octave_idx_type i) const = 0;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	74
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	75 // Length of the index vector.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	76 virtual octave_idx_type length (octave_idx_type n) const = 0;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	77
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	78 // The maximum index + 1. The actual dimension is passed in.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	79 virtual octave_idx_type extent (octave_idx_type n) const = 0;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	80
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	81 // Index class.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	82 virtual idx_class_type idx_class () const { return class_invalid; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	83
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	84 // Sorts, maybe uniqifies, and returns a clone object pointer.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	85 virtual idx_base_rep *sort_uniq_clone (bool uniq = false) = 0;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	86
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	87 // Checks whether the index is colon or a range equivalent to colon.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	88 virtual bool is_colon_equiv (octave_idx_type) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	89 { return false; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	90
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	91 // The original dimensions of this object (used when subscribing by matrices).
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	92 virtual dim_vector orig_dimensions () const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	93 { return dim_vector (); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	94
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	95 // i/o
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	96 virtual std::ostream& print (std::ostream& os) const = 0;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	97
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	98 int count;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	99
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	100 bool err;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	101
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	102 private:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	103
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	104 // No copying!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	105 idx_base_rep (const idx_base_rep&);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	106 };
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	107
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	108 // The magic colon index.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	109 class idx_colon_rep : public idx_base_rep
1879 1b43d3c06c04 [project @ 1996-02-05 16:24:56 by jwe] jwe parents: 1560 diff changeset	110 {
1b43d3c06c04 [project @ 1996-02-05 16:24:56 by jwe] jwe parents: 1560 diff changeset	111 public:
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	112 idx_colon_rep (void) { }
1879 1b43d3c06c04 [project @ 1996-02-05 16:24:56 by jwe] jwe parents: 1560 diff changeset	113
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	114 idx_colon_rep (char c);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	115
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	116 octave_idx_type xelem (octave_idx_type i) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	117 { return i; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	118
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	119 octave_idx_type checkelem (octave_idx_type i) const;
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	120
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	121 octave_idx_type length (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	122 { return n; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	123
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	124 octave_idx_type extent (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	125 { return n; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	126
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	127 idx_class_type idx_class () const { return class_colon; }
1879 1b43d3c06c04 [project @ 1996-02-05 16:24:56 by jwe] jwe parents: 1560 diff changeset	128
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	129 idx_base_rep *sort_uniq_clone (bool = false)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	130 { count++; return this; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	131
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	132 bool is_colon_equiv (octave_idx_type) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	133 { return true; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	134
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	135 std::ostream& print (std::ostream& os) const;
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	136
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	137 private:
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	138
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	139 // No copying!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	140 idx_colon_rep (const idx_colon_rep& idx);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	141 };
2386 4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 2356 diff changeset	142
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	143 // To distinguish the "direct" constructors that blindly trust the data.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	144 enum direct { DIRECT };
2828 92826d6e8bd9 [project @ 1997-03-25 23:41:41 by jwe] jwe parents: 2802 diff changeset	145
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	146 // The integer range index.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	147 class idx_range_rep : public idx_base_rep
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	148 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	149 public:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	150 idx_range_rep (octave_idx_type _start, octave_idx_type _len,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	151 octave_idx_type _step, direct)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	152 : idx_base_rep (), start(_start), len(_len), step(_step) { }
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	153
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	154 idx_range_rep (void)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	155 : start(0), len(1), step(1) { }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	156
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	157 // Zero-based constructor.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	158 idx_range_rep (octave_idx_type _start, octave_idx_type _limit,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	159 octave_idx_type _step);
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	160
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	161 idx_range_rep (const Range&);
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	162
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	163 octave_idx_type xelem (octave_idx_type i) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	164 { return start + i * step; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	165
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	166 octave_idx_type checkelem (octave_idx_type i) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	167
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	168 octave_idx_type length (octave_idx_type) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	169 { return len; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	170
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	171 octave_idx_type extent (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	172 { return std::max (n, (start + 1 + (step < 0 ? 0 : step * (len - 1)))); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	173
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	174 idx_class_type idx_class () const { return class_range; }
1879 1b43d3c06c04 [project @ 1996-02-05 16:24:56 by jwe] jwe parents: 1560 diff changeset	175
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	176 idx_base_rep *sort_uniq_clone (bool uniq = false);
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	177
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	178 bool is_colon_equiv (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	179 { return start == 0 && step == 1 && len == n; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	180
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	181 dim_vector orig_dimensions () const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	182 { return dim_vector (1, len); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	183
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	184 octave_idx_type get_start () const { return start; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	185
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	186 octave_idx_type get_step () const { return step; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	187
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	188 std::ostream& print (std::ostream& os) const;
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	189
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	190 private:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	191
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	192 // No copying!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	193 idx_range_rep (const idx_range_rep& idx);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	194
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	195 octave_idx_type start, len, step;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	196
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	197 };
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	198
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	199 // The integer scalar index.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	200 class idx_scalar_rep : public idx_base_rep
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	201 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	202 public:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	203 idx_scalar_rep (octave_idx_type i, direct)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	204 : data (i) { }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	205
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	206 idx_scalar_rep (void)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	207 : data (0) { }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	208
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	209 // Zero-based constructor.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	210 idx_scalar_rep (octave_idx_type i);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	211
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	212 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	213 idx_scalar_rep (T x);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	214
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	215 octave_idx_type xelem (octave_idx_type) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	216 { return data; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	217
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	218 octave_idx_type checkelem (octave_idx_type i) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	219
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	220 octave_idx_type length (octave_idx_type) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	221 { return 1; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	222
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	223 octave_idx_type extent (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	224 { return std::max (n, data + 1); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	225
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	226 idx_class_type idx_class () const { return class_scalar; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	227
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	228 idx_base_rep *sort_uniq_clone (bool = false)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	229 { count++; return this; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	230
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	231 bool is_colon_equiv (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	232 { return n == 1 && data == 0; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	233
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	234 dim_vector orig_dimensions () const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	235 { return dim_vector (1, 1); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	236
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	237 octave_idx_type get_data () const { return data; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	238
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	239 std::ostream& print (std::ostream& os) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	240
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	241 private:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	242
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	243 // No copying!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	244 idx_scalar_rep (const idx_scalar_rep& idx);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	245
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	246 octave_idx_type data;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	247
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	248 };
4504 f6a61399bc5c [project @ 2003-09-09 17:48:00 by jwe] jwe parents: 4461 diff changeset	249
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	250 // The integer vector index.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	251 class idx_vector_rep : public idx_base_rep
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	252 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	253 public:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	254 // Direct constructor.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	255 idx_vector_rep (octave_idx_type *_data, octave_idx_type _len,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	256 octave_idx_type _ext, const dim_vector& od, direct)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	257 : data (_data), len (_len), ext (_ext), aowner (0), orig_dims (od) { }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	258
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	259 idx_vector_rep (void)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	260 : data (0), len (0), aowner (0)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	261 { }
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	262
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	263 // Zero-based constructor.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	264 idx_vector_rep (const Array<octave_idx_type>& inda);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	265
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	266 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	267 idx_vector_rep (const Array<T>&);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	268
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	269 idx_vector_rep (bool);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	270
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	271 idx_vector_rep (const Array<bool>&);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	272
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	273 ~idx_vector_rep (void);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	274
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	275 octave_idx_type xelem (octave_idx_type i) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	276 { return data[i]; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	277
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	278 octave_idx_type checkelem (octave_idx_type i) const;
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	279
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	280 octave_idx_type length (octave_idx_type) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	281 { return len; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	282
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	283 octave_idx_type extent (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	284 { return std::max (n, ext); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	285
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	286 idx_class_type idx_class () const { return class_vector; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	287
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	288 idx_base_rep *sort_uniq_clone (bool uniq = false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	289
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	290 dim_vector orig_dimensions () const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	291 { return orig_dims; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	292
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	293 const octave_idx_type *get_data () const { return data; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	294
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	295 std::ostream& print (std::ostream& os) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	296
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	297 private:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	298
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	299 // No copying!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	300 idx_vector_rep (const idx_vector_rep& idx);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	301
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	302 const octave_idx_type *data;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	303 octave_idx_type len, ext;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	304
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	305 // This is a trick to allow user-given zero-based arrays to be used as indices
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	306 // without copying. If the following pointer is nonzero, we do not own the data,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	307 // but rather have an Array<octave_idx_type> object that provides us the data.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	308 // Note that we need a pointer because we deferred the Array<T> declaration and
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	309 // we do not want it yet to be defined.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	310
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	311 Array<octave_idx_type> *aowner;
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	312
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	313 dim_vector orig_dims;
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	314 };
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	315
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	316
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	317 idx_vector (idx_base_rep *r) : rep (r) { }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	318
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	319 // The shared empty vector representation (for fast default constructor)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	320 static idx_vector_rep *nil_rep ()
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	321 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	322 static idx_vector_rep ivr;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	323 return &ivr;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	324 }
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	325
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	326 // The shared empty vector representation with the error flag set.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	327 static idx_vector_rep *err_rep ()
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	328 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	329 static idx_vector_rep ivr;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	330 ivr.err = true;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	331 return &ivr;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	332 }
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	333
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	334 // If there was an error in constructing the rep, replace it with empty vector
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	335 // for safety.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	336 void chkerr (void)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	337 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	338 if (rep->err)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	339 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	340 if (--rep->count == 0)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	341 delete rep;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	342 rep = err_rep ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	343 rep->count++;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	344 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	345 }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	346
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	347 public:
1551 8f63ee44726c [project @ 1995-10-08 23:19:55 by jwe] jwe parents: 1315 diff changeset	348
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	349 // Fast empty constructor.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	350 idx_vector (void) : rep (nil_rep ()) { rep->count++; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	351
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	352 // Zero-based constructors (for use from C++).
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	353 idx_vector (octave_idx_type i) : rep (new idx_scalar_rep (i))
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	354 { chkerr (); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	355
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	356 idx_vector (octave_idx_type start, octave_idx_type limit,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	357 octave_idx_type step = 1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	358 : rep (new idx_range_rep (start, limit, step))
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	359 { chkerr (); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	360
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	361 idx_vector (const Array<octave_idx_type>& inda)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	362 : rep (new idx_vector_rep (inda))
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	363 { chkerr (); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	364
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	365 // Colon is best constructed by simply copying (or referencing) this member.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	366 static const idx_vector colon;
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	367
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	368 // or passing ':' here
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	369 idx_vector (char c) : rep (new idx_colon_rep (c)) { chkerr (); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	370
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	371 // Conversion constructors (used by interpreter).
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	372
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	373 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	374 idx_vector (octave_int<T> x) : rep (new idx_scalar_rep (x)) { chkerr (); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	375
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	376 idx_vector (double x) : rep (new idx_scalar_rep (x)) { chkerr (); }
2386 4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 2356 diff changeset	377
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	378 idx_vector (float x) : rep (new idx_scalar_rep (x)) { chkerr (); }
3928 e8627dc4bdf2 [project @ 2002-05-03 19:56:01 by jwe] jwe parents: 3504 diff changeset	379
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	380 // A scalar bool does not necessarily map to scalar index.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	381 idx_vector (bool x) : rep (new idx_vector_rep (x)) { chkerr (); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	382
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	383 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	384 idx_vector (const Array<octave_int<T> >& nda) : rep (new idx_vector_rep (nda))
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	385 { chkerr (); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	386
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	387 idx_vector (const Array<double>& nda) : rep (new idx_vector_rep (nda))
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	388 { chkerr (); }
2828 92826d6e8bd9 [project @ 1997-03-25 23:41:41 by jwe] jwe parents: 2802 diff changeset	389
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	390 idx_vector (const Array<float>& nda) : rep (new idx_vector_rep (nda))
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	391 { chkerr (); }
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	392
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	393 idx_vector (const Array<bool>& nda) : rep (new idx_vector_rep (nda))
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	394 { chkerr (); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	395
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	396 idx_vector (const Range& r)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	397 : rep (new idx_range_rep (r))
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	398 { chkerr (); }
2828 92826d6e8bd9 [project @ 1997-03-25 23:41:41 by jwe] jwe parents: 2802 diff changeset	399
4979 ce4e3d39d05b [project @ 2004-09-08 23:25:05 by jwe] jwe parents: 4938 diff changeset	400 idx_vector (const idx_vector& a) : rep (a.rep) { rep->count++; }
1551 8f63ee44726c [project @ 1995-10-08 23:19:55 by jwe] jwe parents: 1315 diff changeset	401
2802 bca5fbab2e52 [project @ 1997-03-09 11:02:27 by jwe] jwe parents: 2800 diff changeset	402 ~idx_vector (void)
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	403 {
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	404 if (--rep->count == 0)
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	405 delete rep;
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	406 }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	407
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	408 idx_vector& operator = (const idx_vector& a)
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	409 {
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	410 if (this != &a)
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	411 {
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	412 if (--rep->count == 0)
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	413 delete rep;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	414
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	415 rep = a.rep;
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	416 rep->count++;
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	417 }
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	418 return *this;
1551 8f63ee44726c [project @ 1995-10-08 23:19:55 by jwe] jwe parents: 1315 diff changeset	419 }
8f63ee44726c [project @ 1995-10-08 23:19:55 by jwe] jwe parents: 1315 diff changeset	420
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	421 idx_class_type idx_class (void) const { return rep->idx_class (); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	422
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	423 octave_idx_type length (octave_idx_type n = 0) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	424 { return rep->length (n); }
191 b6b4d8c513fe [project @ 1993-10-29 23:41:39 by jwe] jwe parents: 164 diff changeset	425
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	426 octave_idx_type extent (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	427 { return rep->extent (n); }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	428
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	429 octave_idx_type xelem (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	430 { return rep->xelem (n); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	431
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	432 octave_idx_type checkelem (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	433 { return rep->checkelem (n); }
1551 8f63ee44726c [project @ 1995-10-08 23:19:55 by jwe] jwe parents: 1315 diff changeset	434
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	435 octave_idx_type operator () (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	436 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	437 #if defined (BOUNDS_CHECKING)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	438 return rep->checkelem (n);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	439 #else
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	440 return rep->xelem (n);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	441 #endif
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	442 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	443
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	444 operator bool () const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	445 { return ! rep->err; }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	446
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	447 bool is_colon (void) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	448 { return rep->idx_class () == class_colon; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	449
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	450 bool is_scalar (void) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	451 { return rep->idx_class () == class_scalar; }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	452
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	453 bool is_colon_equiv (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	454 { return rep->is_colon_equiv (n); }
3079 eb827eb9b8ff [project @ 1997-08-03 20:44:07 by jwe] jwe parents: 2847 diff changeset	455
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	456 idx_vector sorted (bool uniq = false) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	457 { return idx_vector (rep->sort_uniq_clone (uniq)); }
208 380cb9ccc889 [project @ 1993-11-12 10:35:05 by jwe] jwe parents: 191 diff changeset	458
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	459 dim_vector orig_dimensions (void) const { return rep->orig_dimensions (); }
4504 f6a61399bc5c [project @ 2003-09-09 17:48:00 by jwe] jwe parents: 4461 diff changeset	460
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	461 octave_idx_type orig_rows () const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	462 { return orig_dimensions () (0); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	463
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	464 octave_idx_type orig_columns () const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	465 { return orig_dimensions () (1); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	466
2663 00b2eff19bf5 [project @ 1997-02-08 21:58:49 by jwe] jwe parents: 2386 diff changeset	467 int orig_empty (void) const
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	468 { return (! is_colon () && orig_dimensions().any_zero ()); }
2663 00b2eff19bf5 [project @ 1997-02-08 21:58:49 by jwe] jwe parents: 2386 diff changeset	469
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	470 // i/o
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	471
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	472 std::ostream& print (std::ostream& os) const { return rep->print (os); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	473
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	474 friend std::ostream& operator << (std::ostream& os, const idx_vector& a)
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	475 { return a.print (os); }
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	476
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	477 // Slice with specializations. No checking of bounds!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	478 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	479 // This is equivalent to the following loop (but much faster):
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	480 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	481 // for (octave_idx_type i = 0; i < idx->length (n); i++)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	482 // dest[i] = src[idx(i)];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	483 // return i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	484 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	485 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	486 octave_idx_type
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	487 index (const T src, octave_idx_type n, T dest) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	488 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	489 octave_idx_type len = rep->length (n);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	490 switch (rep->idx_class ())
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	491 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	492 case class_colon:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	493 std::copy (src, src + len, dest);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	494 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	495 case class_range:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	496 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	497 idx_range_rep * r = dynamic_cast<idx_range_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	498 octave_idx_type start = r->get_start (), step = r->get_step ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	499 const T *ssrc = src + start;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	500 if (step == 1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	501 std::copy (ssrc, ssrc + len, dest);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	502 else if (step == -1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	503 std::reverse_copy (ssrc - len + 1, ssrc + 1, dest);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	504 else
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	505 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	506 for (octave_idx_type i = 0, j = 0; i < len; i++, j += step)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	507 dest[i] = ssrc[j];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	508 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	509 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	510 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	511 case class_scalar:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	512 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	513 idx_scalar_rep * r = dynamic_cast<idx_scalar_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	514 dest[0] = src[r->get_data ()];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	515 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	516 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	517 case class_vector:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	518 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	519 idx_vector_rep * r = dynamic_cast<idx_vector_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	520 const octave_idx_type *data = r->get_data ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	521 for (octave_idx_type i = 0; i < len; i++)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	522 dest[i] = src[data[i]];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	523 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	524 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	525 default:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	526 assert (false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	527 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	528 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	529
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	530 return len;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	531 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	532
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	533 // Slice assignment with specializations. No checking of bounds!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	534 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	535 // This is equivalent to the following loop (but much faster):
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	536 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	537 // for (octave_idx_type i = 0; i < idx->length (n); i++)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	538 // dest[idx(i)] = src[i];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	539 // return i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	540 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	541 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	542 octave_idx_type
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	543 assign (const T src, octave_idx_type n, T dest) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	544 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	545 octave_idx_type len = rep->length (n);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	546 switch (rep->idx_class ())
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	547 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	548 case class_colon:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	549 std::copy (src, src + len, dest);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	550 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	551 case class_range:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	552 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	553 idx_range_rep * r = dynamic_cast<idx_range_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	554 octave_idx_type start = r->get_start (), step = r->get_step ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	555 T *sdest = dest + start;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	556 if (step == 1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	557 std::copy (src, src + len, sdest);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	558 else if (step == -1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	559 std::reverse_copy (src, src + len, sdest - len + 1);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	560 else
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	561 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	562 for (octave_idx_type i = 0, j = 0; i < len; i++, j += step)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	563 sdest[j] = src[i];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	564 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	565 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	566 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	567 case class_scalar:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	568 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	569 idx_scalar_rep * r = dynamic_cast<idx_scalar_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	570 dest[r->get_data ()] = src[0];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	571 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	572 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	573 case class_vector:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	574 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	575 idx_vector_rep * r = dynamic_cast<idx_vector_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	576 const octave_idx_type *data = r->get_data ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	577 for (octave_idx_type i = 0; i < len; i++)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	578 dest[data[i]] = src[i];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	579 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	580 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	581 default:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	582 assert (false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	583 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	584 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	585
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	586 return len;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	587 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	588
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	589 // Slice fill with specializations. No checking of bounds!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	590 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	591 // This is equivalent to the following loop (but much faster):
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	592 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	593 // for (octave_idx_type i = 0; i < idx->length (n); i++)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	594 // dest[idx(i)] = val;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	595 // return i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	596 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	597 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	598 octave_idx_type
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	599 fill (const T& val, octave_idx_type n, T *dest) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	600 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	601 octave_idx_type len = rep->length (n);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	602 switch (rep->idx_class ())
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	603 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	604 case class_colon:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	605 std::fill (dest, dest + len, val);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	606 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	607 case class_range:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	608 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	609 idx_range_rep * r = dynamic_cast<idx_range_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	610 octave_idx_type start = r->get_start (), step = r->get_step ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	611 T *sdest = dest + start;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	612 if (step == 1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	613 std::fill (sdest, sdest + len, val);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	614 else if (step == -1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	615 std::fill (sdest - len + 1, sdest + 1, val);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	616 else
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	617 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	618 for (octave_idx_type i = 0, j = 0; i < len; i++, j += step)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	619 sdest[j] = val;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	620 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	621 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	622 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	623 case class_scalar:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	624 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	625 idx_scalar_rep * r = dynamic_cast<idx_scalar_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	626 dest[r->get_data ()] = val;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	627 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	628 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	629 case class_vector:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	630 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	631 idx_vector_rep * r = dynamic_cast<idx_vector_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	632 const octave_idx_type *data = r->get_data ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	633 for (octave_idx_type i = 0; i < len; i++)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	634 dest[data[i]] = val;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	635 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	636 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	637 default:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	638 assert (false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	639 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	640 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	641
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	642 return len;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	643 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	644
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	645 // Generic breakable indexed loop. The loop body should be encapsulated in a
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	646 // single functor body.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	647 // This is equivalent to the following loop (but faster, at least for simple
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	648 // inlined bodies):
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	649 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	650 // for (octave_idx_type i = 0; i < idx->length (n); i++)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	651 // if (body (i, idx(i))) break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	652 // return i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	653 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	654
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	655 template <class Functor>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	656 octave_idx_type
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	657 loop (octave_idx_type n, const Functor& body) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	658 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	659 octave_idx_type len = rep->length (n), ret;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	660 switch (rep->idx_class ())
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	661 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	662 case class_colon:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	663 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	664 octave_idx_type i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	665 for (i = 0; i < len && body (i); i++) ;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	666 ret = i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	667 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	668 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	669 case class_range:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	670 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	671 idx_range_rep * r = dynamic_cast<idx_range_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	672 octave_idx_type start = r->get_start (), step = r->get_step ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	673 octave_idx_type i, j;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	674 if (step == 1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	675 for (i = start, j = start + len; i < j && body (i); i++) ;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	676 else if (step == -1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	677 for (i = start, j = start - len; i > j && body (i); i--) ;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	678 else
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	679 for (i = 0, j = start; i < len && body (j); i++, j += step) ;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	680 ret = i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	681 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	682 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	683 case class_scalar:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	684 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	685 idx_scalar_rep * r = dynamic_cast<idx_scalar_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	686 ret = body (r->get_data ()) ? 1 : 0;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	687 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	688 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	689 case class_vector:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	690 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	691 idx_vector_rep * r = dynamic_cast<idx_vector_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	692 octave_idx_type *data = r->get_data ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	693 octave_idx_type i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	694 for (i = 0; i < len && body (data[i]); i++) ;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	695 ret = i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	696 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	697 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	698 default:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	699 assert (false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	700 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	701 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	702
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	703 return ret;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	704 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	705
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	706 // Rationale:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	707 // This method is the key to "smart indexing". When indexing cartesian
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	708 // arrays, sometimes consecutive index vectors can be reduced into a single
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	709 // index. If rows (A) = k and i.maybe_reduce (j) gives k, then A(i,j)(:) is
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	710 // equal to A(k)(:).
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	711
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	712 // If the next index can be reduced, returns true and updates this.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	713 bool maybe_reduce (octave_idx_type n, const idx_vector& j,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	714 octave_idx_type nj);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	715
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	716 bool is_cont_range (octave_idx_type n,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	717 octave_idx_type& l, octave_idx_type& u) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	718
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	719 idx_vector
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	720 complement (octave_idx_type n) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	721
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	722 // FIXME: These are here for compatibility. They should be removed when no
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	723 // longer in use.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	724
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	725 octave_idx_type elem (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	726 { return (*this) (n); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	727
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	728 bool is_colon_equiv (octave_idx_type n, int) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	729 { return is_colon_equiv (n); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	730
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	731 octave_idx_type freeze (octave_idx_type z_len, const char *tag, bool resize_ok = false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	732
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	733 void sort (bool uniq = false)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	734 { *this = sorted (uniq); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	735
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	736 octave_idx_type ones_count () const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	737
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	738 octave_idx_type max () const { return extent (1) - 1; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	739
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	740 private:
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	741
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	742 idx_base_rep *rep;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	743
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	744 };
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	745
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	746 #endif
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	747
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	748 /*
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	749 ;;; Local Variables: ***
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	750 ;;; mode: C++ ***
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	751 ;;; End: ***
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	752 */

Mercurial > octave

annotate liboctave/idx-vector.h @ 8290:7cbe01c21986