octave: liboctave/idx-vector.h annotate

annotate liboctave/idx-vector.h @ 13735:2cd903770748

Fix broken gcc compilation due to changeset d13b6086aea9. * idx-vector.cc: Don't include <memory>. * idx-vector.h: Move here.

author	Michael Goffioul <michael.goffioul@gmail.com>
date	Thu, 20 Oct 2011 14:08:25 +0100
parents	d13b6086aea9
children	72c96de7a403

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
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 11503 diff changeset	3 Copyright (C) 1993-2011 John W. Eaton
fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 11503 diff changeset	4 Copyright (C) 2008-2009 Jaroslav Hajek
9601 a9b37bae1802 add a couple of missing copyright statements Jaroslav Hajek <highegg@gmail.com> parents: 9479 diff changeset	5 Copyright (C) 2009 VZLU Prague
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
8950 d865363208d6 include <iosfwd> instead of <iostream> in header files John W. Eaton <jwe@octave.org> parents: 8934 diff changeset	28 #include <cassert>
d865363208d6 include <iosfwd> instead of <iostream> in header files John W. Eaton <jwe@octave.org> parents: 8934 diff changeset	29
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	30 #include <algorithm>
8950 d865363208d6 include <iosfwd> instead of <iostream> in header files John W. Eaton <jwe@octave.org> parents: 8934 diff changeset	31 #include <iosfwd>
13735 2cd903770748 Fix broken gcc compilation due to changeset d13b6086aea9. Michael Goffioul <michael.goffioul@gmail.com> parents: 13729 diff changeset	32 #include <memory>
3503 d14c483b3c12 [project @ 2000-02-01 04:06:07 by jwe] jwe parents: 3145 diff changeset	33
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	34 #include "dim-vector.h"
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	35 #include "oct-inttypes.h"
8680 a2bcd96b9344 use allocators for idx_vector_rep subclasses Jaroslav Hajek <highegg@gmail.com> parents: 8384 diff changeset	36 #include "oct-alloc.h"
9773 01f897d8a130 optimize memory manipulation by arrays & indexing Jaroslav Hajek <highegg@gmail.com> parents: 9725 diff changeset	37 #include "oct-mem.h"
12125 a21a3875ca83 implement a common class for reference counts Jaroslav Hajek <highegg@gmail.com> parents: 11586 diff changeset	38 #include "oct-refcount.h"
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	39
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	40 template<class T> class Array;
8384 a99b9113c58c optimize sparse bool indexing Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	41 template<class T> class Sparse;
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	42 class Range;
4504 f6a61399bc5c [project @ 2003-09-09 17:48:00 by jwe] jwe parents: 4461 diff changeset	43
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	44 // Design rationale:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	45 // 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	46 // 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	47 // 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	48 // 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	49 // type out of loops.
164 e2c950dd96d2 [project @ 1993-10-18 19:32:00 by jwe] jwe parents: 1 diff changeset	50
453 393e95f46b51 [project @ 1994-06-06 00:05:20 by jwe] jwe parents: 434 diff changeset	51 class
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5781 diff changeset	52 OCTAVE_API
453 393e95f46b51 [project @ 1994-06-06 00:05:20 by jwe] jwe parents: 434 diff changeset	53 idx_vector
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	54 {
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	55 public:
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	56
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	57 enum idx_class_type
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 class_invalid = -1,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	60 class_colon = 0,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	61 class_range,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	62 class_scalar,
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	63 class_vector,
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	64 class_mask
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	65 };
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	66
13729 d13b6086aea9 Make std::auto_ptr friend of idx_vector to allow access to private class idx_vector_rep. Michael Goffioul <michael.goffioul@gmail.com> parents: 12125 diff changeset	67 template<class T> friend class std::auto_ptr;
d13b6086aea9 Make std::auto_ptr friend of idx_vector to allow access to private class idx_vector_rep. Michael Goffioul <michael.goffioul@gmail.com> parents: 12125 diff changeset	68
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	69 private:
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	70
8319 c374691576f6 Fix for MSVC compilation Michael Goffioul <michael.goffioul@gmail.com> parents: 8298 diff changeset	71 class OCTAVE_API idx_base_rep
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	72 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	73 public:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	74 idx_base_rep (void) : count (1), err (false) { }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	75
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	76 virtual ~idx_base_rep (void) { }
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 // Non-range-checking element query.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	79 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	80
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	81 // Range-checking element query.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	82 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	83
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	84 // Length of the index vector.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	85 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	86
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	87 // The maximum index + 1. The actual dimension is passed in.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	88 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	89
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	90 // Index class.
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	91 virtual idx_class_type idx_class (void) const { return class_invalid; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	92
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	93 // Sorts, maybe uniqifies, and returns a clone object pointer.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	94 virtual idx_base_rep *sort_uniq_clone (bool uniq = false) = 0;
10273 3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	95 // Sorts, and returns a sorting permutation (aka Array::sort).
3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	96 virtual idx_base_rep *sort_idx (Array<octave_idx_type>&) = 0;
8290 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 // Checks whether the index is colon or a range equivalent to colon.
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	99 virtual bool is_colon_equiv (octave_idx_type) const { return false; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	100
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	101 // The original dimensions of this object (used when subscribing by matrices).
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	102 virtual dim_vector orig_dimensions (void) const { return dim_vector (); }
8290 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 // i/o
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	105 virtual std::ostream& print (std::ostream& os) const = 0;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	106
10339 de2d43bcb083 optimize some lazy index operations Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	107 virtual Array<octave_idx_type> as_array (void);
de2d43bcb083 optimize some lazy index operations Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	108
12125 a21a3875ca83 implement a common class for reference counts Jaroslav Hajek <highegg@gmail.com> parents: 11586 diff changeset	109 octave_refcount<int> count;
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	110
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	111 bool err;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	112
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	113 private:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	114
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	115 // No copying!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	116 idx_base_rep (const idx_base_rep&);
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	117 idx_base_rep& operator = (const idx_base_rep&);
8290 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
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	120 // The magic colon index.
8319 c374691576f6 Fix for MSVC compilation Michael Goffioul <michael.goffioul@gmail.com> parents: 8298 diff changeset	121 class OCTAVE_API idx_colon_rep : public idx_base_rep
1879 1b43d3c06c04 [project @ 1996-02-05 16:24:56 by jwe] jwe parents: 1560 diff changeset	122 {
1b43d3c06c04 [project @ 1996-02-05 16:24:56 by jwe] jwe parents: 1560 diff changeset	123 public:
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	124 idx_colon_rep (void) { }
1879 1b43d3c06c04 [project @ 1996-02-05 16:24:56 by jwe] jwe parents: 1560 diff changeset	125
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	126 idx_colon_rep (char c);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	127
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	128 octave_idx_type xelem (octave_idx_type i) const { return i; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	129
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	130 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	131
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	132 octave_idx_type length (octave_idx_type n) const { return n; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	133
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	134 octave_idx_type extent (octave_idx_type n) const { return n; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	135
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	136 idx_class_type idx_class (void) const { return class_colon; }
1879 1b43d3c06c04 [project @ 1996-02-05 16:24:56 by jwe] jwe parents: 1560 diff changeset	137
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	138 idx_base_rep *sort_uniq_clone (bool = false)
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	139 { count++; return this; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	140
10273 3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	141 idx_base_rep *sort_idx (Array<octave_idx_type>&);
3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	142
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	143 bool is_colon_equiv (octave_idx_type) const { return true; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	144
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	145 std::ostream& print (std::ostream& os) const;
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	146
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	147 private:
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	148
8680 a2bcd96b9344 use allocators for idx_vector_rep subclasses Jaroslav Hajek <highegg@gmail.com> parents: 8384 diff changeset	149 DECLARE_OCTAVE_ALLOCATOR
a2bcd96b9344 use allocators for idx_vector_rep subclasses Jaroslav Hajek <highegg@gmail.com> parents: 8384 diff changeset	150
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	151 // No copying!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	152 idx_colon_rep (const idx_colon_rep& idx);
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	153 idx_colon_rep& operator = (const idx_colon_rep& idx);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	154 };
2386 4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 2356 diff changeset	155
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	156 // To distinguish the "direct" constructors that blindly trust the data.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	157 enum direct { DIRECT };
2828 92826d6e8bd9 [project @ 1997-03-25 23:41:41 by jwe] jwe parents: 2802 diff changeset	158
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	159 // The integer range index.
8319 c374691576f6 Fix for MSVC compilation Michael Goffioul <michael.goffioul@gmail.com> parents: 8298 diff changeset	160 class OCTAVE_API idx_range_rep : public idx_base_rep
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	161 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	162 public:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	163 idx_range_rep (octave_idx_type _start, octave_idx_type _len,
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	164 octave_idx_type _step, direct)
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	165 : 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	166
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	167 idx_range_rep (void)
8298 7e87d3d76a56 fix extent query for empty ranges Jaroslav Hajek <highegg@gmail.com> parents: 8290 diff changeset	168 : start(0), len(0), step(1) { }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	169
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	170 // Zero-based constructor.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	171 idx_range_rep (octave_idx_type _start, octave_idx_type _limit,
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	172 octave_idx_type _step);
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_range_rep (const Range&);
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	175
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	176 octave_idx_type xelem (octave_idx_type i) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	177 { return start + i * step; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	178
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	179 octave_idx_type checkelem (octave_idx_type i) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	180
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	181 octave_idx_type length (octave_idx_type) const { return len; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	182
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	183 octave_idx_type extent (octave_idx_type n) const
8298 7e87d3d76a56 fix extent query for empty ranges Jaroslav Hajek <highegg@gmail.com> parents: 8290 diff changeset	184 { return len ? std::max (n, (start + 1 + (step < 0 ? 0 : step * (len - 1)))) : n; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	185
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	186 idx_class_type idx_class (void) const { return class_range; }
1879 1b43d3c06c04 [project @ 1996-02-05 16:24:56 by jwe] jwe parents: 1560 diff changeset	187
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	188 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	189
10273 3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	190 idx_base_rep *sort_idx (Array<octave_idx_type>&);
3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	191
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	192 bool is_colon_equiv (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	193 { return start == 0 && step == 1 && len == n; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	194
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	195 dim_vector orig_dimensions (void) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	196 { return dim_vector (1, len); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	197
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	198 octave_idx_type get_start (void) const { return start; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	199
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	200 octave_idx_type get_step (void) const { return step; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	201
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	202 std::ostream& print (std::ostream& os) const;
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	203
9894 83bd7f34f9da improve idx_vector->octave_value conversion Jaroslav Hajek <highegg@gmail.com> parents: 9886 diff changeset	204 Range unconvert (void) const;
83bd7f34f9da improve idx_vector->octave_value conversion Jaroslav Hajek <highegg@gmail.com> parents: 9886 diff changeset	205
10339 de2d43bcb083 optimize some lazy index operations Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	206 Array<octave_idx_type> as_array (void);
de2d43bcb083 optimize some lazy index operations Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	207
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	208 private:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	209
8680 a2bcd96b9344 use allocators for idx_vector_rep subclasses Jaroslav Hajek <highegg@gmail.com> parents: 8384 diff changeset	210 DECLARE_OCTAVE_ALLOCATOR
a2bcd96b9344 use allocators for idx_vector_rep subclasses Jaroslav Hajek <highegg@gmail.com> parents: 8384 diff changeset	211
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	212 // No copying!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	213 idx_range_rep (const idx_range_rep& idx);
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	214 idx_range_rep& operator = (const idx_range_rep& idx);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	215
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	216 octave_idx_type start, len, step;
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 };
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 // The integer scalar index.
8319 c374691576f6 Fix for MSVC compilation Michael Goffioul <michael.goffioul@gmail.com> parents: 8298 diff changeset	221 class OCTAVE_API idx_scalar_rep : public idx_base_rep
8290 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 public:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	224 idx_scalar_rep (octave_idx_type i, direct)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	225 : data (i) { }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	226
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	227 idx_scalar_rep (void)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	228 : data (0) { }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	229
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	230 // Zero-based constructor.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	231 idx_scalar_rep (octave_idx_type i);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	232
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	233 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	234 idx_scalar_rep (T x);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	235
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	236 octave_idx_type xelem (octave_idx_type) const { return data; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	237
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	238 octave_idx_type checkelem (octave_idx_type i) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	239
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	240 octave_idx_type length (octave_idx_type) const { return 1; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	241
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	242 octave_idx_type extent (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	243 { return std::max (n, data + 1); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	244
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	245 idx_class_type idx_class (void) const { return class_scalar; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	246
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	247 idx_base_rep *sort_uniq_clone (bool = false)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	248 { count++; return this; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	249
10273 3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	250 idx_base_rep *sort_idx (Array<octave_idx_type>&);
3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	251
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	252 bool is_colon_equiv (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	253 { return n == 1 && data == 0; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	254
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	255 dim_vector orig_dimensions (void) const { return dim_vector (1, 1); }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	256
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	257 octave_idx_type get_data (void) const { return data; }
8290 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 std::ostream& print (std::ostream& os) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	260
9894 83bd7f34f9da improve idx_vector->octave_value conversion Jaroslav Hajek <highegg@gmail.com> parents: 9886 diff changeset	261 double unconvert (void) const;
83bd7f34f9da improve idx_vector->octave_value conversion Jaroslav Hajek <highegg@gmail.com> parents: 9886 diff changeset	262
10339 de2d43bcb083 optimize some lazy index operations Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	263 Array<octave_idx_type> as_array (void);
de2d43bcb083 optimize some lazy index operations Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	264
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	265 private:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	266
8680 a2bcd96b9344 use allocators for idx_vector_rep subclasses Jaroslav Hajek <highegg@gmail.com> parents: 8384 diff changeset	267 DECLARE_OCTAVE_ALLOCATOR
a2bcd96b9344 use allocators for idx_vector_rep subclasses Jaroslav Hajek <highegg@gmail.com> parents: 8384 diff changeset	268
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	269 // No copying!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	270 idx_scalar_rep (const idx_scalar_rep& idx);
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	271 idx_scalar_rep& operator = (const idx_scalar_rep& idx);
8290 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 octave_idx_type data;
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 };
4504 f6a61399bc5c [project @ 2003-09-09 17:48:00 by jwe] jwe parents: 4461 diff changeset	276
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	277 // The integer vector index.
8319 c374691576f6 Fix for MSVC compilation Michael Goffioul <michael.goffioul@gmail.com> parents: 8298 diff changeset	278 class OCTAVE_API idx_vector_rep : public idx_base_rep
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	279 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	280 public:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	281 // Direct constructor.
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	282 idx_vector_rep (octave_idx_type *_data, octave_idx_type _len,
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	283 octave_idx_type _ext, const dim_vector& od, direct)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	284 : 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	285
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	286 idx_vector_rep (void)
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	287 : data (0), len (0), ext (0), aowner (0), orig_dims ()
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	288 { }
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	289
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	290 // Zero-based constructor.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	291 idx_vector_rep (const Array<octave_idx_type>& inda);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	292
10285 22a7913bbeb5 optimize return values of find and sort Jaroslav Hajek <highegg@gmail.com> parents: 10273 diff changeset	293 idx_vector_rep (const Array<octave_idx_type>& inda,
22a7913bbeb5 optimize return values of find and sort Jaroslav Hajek <highegg@gmail.com> parents: 10273 diff changeset	294 octave_idx_type _ext, direct);
22a7913bbeb5 optimize return values of find and sort Jaroslav Hajek <highegg@gmail.com> parents: 10273 diff changeset	295
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	296 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	297 idx_vector_rep (const Array<T>&);
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 idx_vector_rep (bool);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	300
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	301 idx_vector_rep (const Array<bool>&, octave_idx_type = -1);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	302
8384 a99b9113c58c optimize sparse bool indexing Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	303 idx_vector_rep (const Sparse<bool>&);
a99b9113c58c optimize sparse bool indexing Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	304
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	305 ~idx_vector_rep (void);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	306
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	307 octave_idx_type xelem (octave_idx_type i) const { return data[i]; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	308
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	309 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	310
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	311 octave_idx_type length (octave_idx_type) const { return len; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	312
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	313 octave_idx_type extent (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	314 { return std::max (n, ext); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	315
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	316 idx_class_type idx_class (void) const { return class_vector; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	317
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	318 idx_base_rep *sort_uniq_clone (bool uniq = false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	319
10273 3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	320 idx_base_rep *sort_idx (Array<octave_idx_type>&);
3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	321
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	322 dim_vector orig_dimensions (void) const { return orig_dims; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	323
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	324 const octave_idx_type *get_data (void) const { return data; }
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	325
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	326 std::ostream& print (std::ostream& os) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	327
9894 83bd7f34f9da improve idx_vector->octave_value conversion Jaroslav Hajek <highegg@gmail.com> parents: 9886 diff changeset	328 Array<double> unconvert (void) const;
83bd7f34f9da improve idx_vector->octave_value conversion Jaroslav Hajek <highegg@gmail.com> parents: 9886 diff changeset	329
10339 de2d43bcb083 optimize some lazy index operations Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	330 Array<octave_idx_type> as_array (void);
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	331
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	332 private:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	333
8680 a2bcd96b9344 use allocators for idx_vector_rep subclasses Jaroslav Hajek <highegg@gmail.com> parents: 8384 diff changeset	334 DECLARE_OCTAVE_ALLOCATOR
a2bcd96b9344 use allocators for idx_vector_rep subclasses Jaroslav Hajek <highegg@gmail.com> parents: 8384 diff changeset	335
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	336 // No copying!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	337 idx_vector_rep (const idx_vector_rep& idx);
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	338 idx_vector_rep& operator = (const idx_vector_rep& idx);
8290 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 const octave_idx_type *data;
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	341 octave_idx_type len;
b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	342 octave_idx_type ext;
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	343
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	344 // This is a trick to allow user-given zero-based arrays to be used
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	345 // as indices without copying. If the following pointer is nonzero,
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	346 // we do not own the data, but rather have an Array<octave_idx_type>
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	347 // object that provides us the data. Note that we need a pointer
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	348 // because we deferred the Array<T> declaration and we do not want
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	349 // it yet to be defined.
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	350
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	351 Array<octave_idx_type> *aowner;
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	352
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	353 dim_vector orig_dims;
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	354 };
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	355
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	356 // The logical mask index.
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	357 class OCTAVE_API idx_mask_rep : public idx_base_rep
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	358 {
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	359 public:
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	360 // Direct constructor.
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	361 idx_mask_rep (bool *_data, octave_idx_type _len,
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	362 octave_idx_type _ext, const dim_vector& od, direct)
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	363 : data (_data), len (_len), ext (_ext), lsti (-1), lste (-1),
b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	364 aowner (0), orig_dims (od) { }
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	365
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	366 idx_mask_rep (void)
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	367 : data (0), len (0), ext (0), lsti (-1), lste (-1), aowner (0),
b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	368 orig_dims ()
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	369 { }
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	370
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	371 idx_mask_rep (bool);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	372
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	373 idx_mask_rep (const Array<bool>&, octave_idx_type = -1);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	374
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	375 ~idx_mask_rep (void);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	376
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	377 octave_idx_type xelem (octave_idx_type i) const;
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	378
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	379 octave_idx_type checkelem (octave_idx_type i) const;
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	380
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	381 octave_idx_type length (octave_idx_type) const { return len; }
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	382
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	383 octave_idx_type extent (octave_idx_type n) const
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	384 { return std::max (n, ext); }
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	385
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	386 idx_class_type idx_class (void) const { return class_mask; }
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	387
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	388 idx_base_rep *sort_uniq_clone (bool = false)
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	389 { count++; return this; }
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	390
10273 3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	391 idx_base_rep *sort_idx (Array<octave_idx_type>&);
3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	392
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	393 dim_vector orig_dimensions (void) const { return orig_dims; }
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	394
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	395 bool is_colon_equiv (octave_idx_type n) const
9886 cddd5c3d5f04 fix & extend special-case optimizations for indexed assignment Jaroslav Hajek <highegg@gmail.com> parents: 9879 diff changeset	396 { return len == n && ext == n; }
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	397
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	398 const bool *get_data (void) const { return data; }
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	399
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	400 std::ostream& print (std::ostream& os) const;
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	401
9894 83bd7f34f9da improve idx_vector->octave_value conversion Jaroslav Hajek <highegg@gmail.com> parents: 9886 diff changeset	402 Array<bool> unconvert (void) const;
83bd7f34f9da improve idx_vector->octave_value conversion Jaroslav Hajek <highegg@gmail.com> parents: 9886 diff changeset	403
10425 0677c5d80b77 rewrite 1D sparse indexing Jaroslav Hajek <highegg@gmail.com> parents: 10339 diff changeset	404 Array<octave_idx_type> as_array (void);
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	405
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	406 private:
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	407
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	408 DECLARE_OCTAVE_ALLOCATOR
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	409
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	410 // No copying!
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	411 idx_mask_rep (const idx_mask_rep& idx);
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	412 idx_mask_rep& operator = (const idx_mask_rep& idx);
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	413
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	414 const bool *data;
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	415 octave_idx_type len;
b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	416 octave_idx_type ext;
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	417
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	418 // FIXME: I'm not sure if this is a good design. Maybe it would be
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	419 // better to employ some sort of generalized iteration scheme.
11503 b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	420 mutable octave_idx_type lsti;
b3f511c48b14 idx-vector.h, idx-vector.cc: data member initialization and constructor fixes John W. Eaton <jwe@octave.org> parents: 10516 diff changeset	421 mutable octave_idx_type lste;
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	422
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	423 // This is a trick to allow user-given mask arrays to be used as
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	424 // indices without copying. If the following pointer is nonzero, we
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	425 // do not own the data, but rather have an Array<bool> object that
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	426 // provides us the data. Note that we need a pointer because we
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	427 // deferred the Array<T> declaration and we do not want it yet to be
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	428 // defined.
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	429
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	430 Array<bool> *aowner;
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	431
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	432 dim_vector orig_dims;
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	433 };
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	434
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	435 idx_vector (idx_base_rep *r) : rep (r) { }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	436
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	437 // The shared empty vector representation (for fast default
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	438 // constructor).
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	439 static idx_vector_rep *nil_rep (void)
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	440 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	441 static idx_vector_rep ivr;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	442 return &ivr;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	443 }
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	444
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	445 // The shared empty vector representation with the error flag set.
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	446 static idx_vector_rep *err_rep (void)
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	447 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	448 static idx_vector_rep ivr;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	449 ivr.err = true;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	450 return &ivr;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	451 }
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	452
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	453 // If there was an error in constructing the rep, replace it with
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	454 // empty vector for safety.
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	455 void chkerr (void)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	456 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	457 if (rep->err)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	458 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	459 if (--rep->count == 0)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	460 delete rep;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	461 rep = err_rep ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	462 rep->count++;
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 }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	465
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	466 public:
1551 8f63ee44726c [project @ 1995-10-08 23:19:55 by jwe] jwe parents: 1315 diff changeset	467
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	468 // Fast empty constructor.
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	469 idx_vector (void) : rep (nil_rep ()) { rep->count++; }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	470
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	471 // Zero-based constructors (for use from C++).
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	472 idx_vector (octave_idx_type i) : rep (new idx_scalar_rep (i))
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	473 { chkerr (); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	474
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	475 idx_vector (octave_idx_type start, octave_idx_type limit,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	476 octave_idx_type step = 1)
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	477 : rep (new idx_range_rep (start, limit, step))
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	478 { chkerr (); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	479
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	480 static idx_vector
9479 d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	481 make_range (octave_idx_type start, octave_idx_type step,
d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	482 octave_idx_type len)
d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	483 {
d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	484 return idx_vector (new idx_range_rep (start, len, step, DIRECT));
d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	485 }
d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	486
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	487 idx_vector (const Array<octave_idx_type>& inda)
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	488 : rep (new idx_vector_rep (inda))
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	489 { chkerr (); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	490
10285 22a7913bbeb5 optimize return values of find and sort Jaroslav Hajek <highegg@gmail.com> parents: 10273 diff changeset	491 // Directly pass extent, no checking.
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	492 idx_vector (const Array<octave_idx_type>& inda, octave_idx_type ext)
10285 22a7913bbeb5 optimize return values of find and sort Jaroslav Hajek <highegg@gmail.com> parents: 10273 diff changeset	493 : rep (new idx_vector_rep (inda, ext, DIRECT))
22a7913bbeb5 optimize return values of find and sort Jaroslav Hajek <highegg@gmail.com> parents: 10273 diff changeset	494 { }
22a7913bbeb5 optimize return values of find and sort Jaroslav Hajek <highegg@gmail.com> parents: 10273 diff changeset	495
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	496 // 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	497 static const idx_vector colon;
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	498
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	499 // or passing ':' here
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	500 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	501
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	502 // Conversion constructors (used by interpreter).
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	503
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	504 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	505 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	506
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	507 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	508
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	509 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	510
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	511 // A scalar bool does not necessarily map to scalar index.
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	512 idx_vector (bool x) : rep (new idx_mask_rep (x)) { chkerr (); }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	513
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	514 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	515 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	516 { chkerr (); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	517
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	518 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	519 { chkerr (); }
2828 92826d6e8bd9 [project @ 1997-03-25 23:41:41 by jwe] jwe parents: 2802 diff changeset	520
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	521 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	522 { chkerr (); }
4938 703d97b89507 [project @ 2004-08-09 20:21:59 by jwe] jwe parents: 4653 diff changeset	523
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	524 idx_vector (const Array<bool>& nda);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	525
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	526 idx_vector (const Range& r)
10489 d47802f0e557 Back out changeset 4e64fbbd5c58 Jaroslav Hajek <highegg@gmail.com> parents: 10487 diff changeset	527 : rep (new idx_range_rep (r))
d47802f0e557 Back out changeset 4e64fbbd5c58 Jaroslav Hajek <highegg@gmail.com> parents: 10487 diff changeset	528 { chkerr (); }
2828 92826d6e8bd9 [project @ 1997-03-25 23:41:41 by jwe] jwe parents: 2802 diff changeset	529
8384 a99b9113c58c optimize sparse bool indexing Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	530 idx_vector (const Sparse<bool>& nda) : rep (new idx_vector_rep (nda))
a99b9113c58c optimize sparse bool indexing Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	531 { chkerr (); }
a99b9113c58c optimize sparse bool indexing Jaroslav Hajek <highegg@gmail.com> parents: 8367 diff changeset	532
4979 ce4e3d39d05b [project @ 2004-09-08 23:25:05 by jwe] jwe parents: 4938 diff changeset	533 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	534
2802 bca5fbab2e52 [project @ 1997-03-09 11:02:27 by jwe] jwe parents: 2800 diff changeset	535 ~idx_vector (void)
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	536 {
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	537 if (--rep->count == 0)
10312 cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10285 diff changeset	538 delete rep;
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	539 }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	540
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	541 idx_vector& operator = (const idx_vector& a)
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	542 {
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	543 if (this != &a)
10312 cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10285 diff changeset	544 {
cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10285 diff changeset	545 if (--rep->count == 0)
cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10285 diff changeset	546 delete rep;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	547
10312 cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10285 diff changeset	548 rep = a.rep;
cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10285 diff changeset	549 rep->count++;
cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10285 diff changeset	550 }
1560 27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	551 return *this;
1551 8f63ee44726c [project @ 1995-10-08 23:19:55 by jwe] jwe parents: 1315 diff changeset	552 }
8f63ee44726c [project @ 1995-10-08 23:19:55 by jwe] jwe parents: 1315 diff changeset	553
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	554 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	555
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	556 octave_idx_type length (octave_idx_type n = 0) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	557 { return rep->length (n); }
191 b6b4d8c513fe [project @ 1993-10-29 23:41:39 by jwe] jwe parents: 164 diff changeset	558
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	559 octave_idx_type extent (octave_idx_type n) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	560 { return rep->extent (n); }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	561
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	562 octave_idx_type xelem (octave_idx_type n) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	563 { return rep->xelem (n); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	564
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	565 octave_idx_type checkelem (octave_idx_type n) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	566 { return rep->checkelem (n); }
1551 8f63ee44726c [project @ 1995-10-08 23:19:55 by jwe] jwe parents: 1315 diff changeset	567
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	568 octave_idx_type operator () (octave_idx_type n) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	569 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	570 #if defined (BOUNDS_CHECKING)
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	571 return rep->checkelem (n);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	572 #else
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	573 return rep->xelem (n);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	574 #endif
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	575 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	576
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	577 operator bool (void) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	578 { return ! rep->err; }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	579
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	580 bool is_colon (void) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	581 { return rep->idx_class () == class_colon; }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	582
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	583 bool is_scalar (void) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	584 { return rep->idx_class () == class_scalar; }
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	585
9479 d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	586 bool is_range (void) const
d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	587 { return rep->idx_class () == class_range; }
d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	588
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	589 bool is_colon_equiv (octave_idx_type n) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	590 { return rep->is_colon_equiv (n); }
3079 eb827eb9b8ff [project @ 1997-08-03 20:44:07 by jwe] jwe parents: 2847 diff changeset	591
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	592 idx_vector sorted (bool uniq = false) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	593 { return idx_vector (rep->sort_uniq_clone (uniq)); }
208 380cb9ccc889 [project @ 1993-11-12 10:35:05 by jwe] jwe parents: 191 diff changeset	594
10273 3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	595 idx_vector sorted (Array<octave_idx_type>& sidx) const
3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	596 { return idx_vector (rep->sort_idx (sidx)); }
3a8c13b71612 implement special-case optimization for sort of index vectors Jaroslav Hajek <highegg@gmail.com> parents: 10158 diff changeset	597
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	598 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	599
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	600 octave_idx_type orig_rows (void) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	601 { return orig_dimensions () (0); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	602
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	603 octave_idx_type orig_columns (void) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	604 { return orig_dimensions () (1); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	605
2663 00b2eff19bf5 [project @ 1997-02-08 21:58:49 by jwe] jwe parents: 2386 diff changeset	606 int orig_empty (void) const
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	607 { return (! is_colon () && orig_dimensions().any_zero ()); }
2663 00b2eff19bf5 [project @ 1997-02-08 21:58:49 by jwe] jwe parents: 2386 diff changeset	608
4653 14ab7b05a572 [project @ 2003-11-24 05:02:46 by jwe] jwe parents: 4650 diff changeset	609 // i/o
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	610
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	611 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	612
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 3503 diff changeset	613 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	614 { return a.print (os); }
27a03373de41 [project @ 1995-10-12 07:22:26 by jwe] jwe parents: 1553 diff changeset	615
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	616 // Slice with specializations. No checking of bounds!
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 // This is equivalent to the following loop (but much faster):
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	619 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	620 // 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	621 // dest[i] = src[idx(i)];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	622 // return i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	623 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	624 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	625 octave_idx_type
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	626 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	627 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	628 octave_idx_type len = rep->length (n);
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	629
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	630 switch (rep->idx_class ())
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	631 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	632 case class_colon:
9780 6dafc60dde31 rename oct-mem functions Jaroslav Hajek <highegg@gmail.com> parents: 9773 diff changeset	633 copy_or_memcpy (len, src, dest);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	634 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	635
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	636 case class_range:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	637 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	638 idx_range_rep * r = dynamic_cast<idx_range_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	639 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	640 const T *ssrc = src + start;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	641 if (step == 1)
9780 6dafc60dde31 rename oct-mem functions Jaroslav Hajek <highegg@gmail.com> parents: 9773 diff changeset	642 copy_or_memcpy (len, ssrc, dest);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	643 else if (step == -1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	644 std::reverse_copy (ssrc - len + 1, ssrc + 1, dest);
8972 5fa53d1b6247 more indexing optimizations Jaroslav Hajek <highegg@gmail.com> parents: 8950 diff changeset	645 else if (step == 0)
5fa53d1b6247 more indexing optimizations Jaroslav Hajek <highegg@gmail.com> parents: 8950 diff changeset	646 std::fill_n (dest, len, *ssrc);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	647 else
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	648 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	649 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	650 dest[i] = ssrc[j];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	651 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	652 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	653 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	654
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	655 case class_scalar:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	656 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	657 idx_scalar_rep * r = dynamic_cast<idx_scalar_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	658 dest[0] = src[r->get_data ()];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	659 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	660 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	661
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	662 case class_vector:
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 idx_vector_rep * r = dynamic_cast<idx_vector_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	665 const octave_idx_type *data = r->get_data ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	666 for (octave_idx_type i = 0; i < len; i++)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	667 dest[i] = src[data[i]];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	668 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	669 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	670
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	671 case class_mask:
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	672 {
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	673 idx_mask_rep * r = dynamic_cast<idx_mask_rep *> (rep);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	674 const bool *data = r->get_data ();
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	675 octave_idx_type ext = r->extent (0);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	676 for (octave_idx_type i = 0; i < ext; i++)
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	677 if (data[i]) *dest++ = src[i];
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	678 }
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	679 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	680
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	681 default:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	682 assert (false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	683 break;
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
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	686 return len;
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
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	689 // Slice assignment with specializations. No checking of bounds!
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 // This is equivalent to the following loop (but much faster):
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	692 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	693 // 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	694 // dest[idx(i)] = src[i];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	695 // return 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 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	698 octave_idx_type
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	699 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	700 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	701 octave_idx_type len = rep->length (n);
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	702
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	703 switch (rep->idx_class ())
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 case class_colon:
9780 6dafc60dde31 rename oct-mem functions Jaroslav Hajek <highegg@gmail.com> parents: 9773 diff changeset	706 copy_or_memcpy (len, src, dest);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	707 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	708
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	709 case class_range:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	710 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	711 idx_range_rep * r = dynamic_cast<idx_range_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	712 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	713 T *sdest = dest + start;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	714 if (step == 1)
9780 6dafc60dde31 rename oct-mem functions Jaroslav Hajek <highegg@gmail.com> parents: 9773 diff changeset	715 copy_or_memcpy (len, src, sdest);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	716 else if (step == -1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	717 std::reverse_copy (src, src + len, sdest - len + 1);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	718 else
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	719 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	720 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	721 sdest[j] = src[i];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	722 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	723 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	724 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	725
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	726 case class_scalar:
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 idx_scalar_rep * r = dynamic_cast<idx_scalar_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	729 dest[r->get_data ()] = src[0];
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 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	732
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	733 case class_vector:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	734 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	735 idx_vector_rep * r = dynamic_cast<idx_vector_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	736 const octave_idx_type *data = r->get_data ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	737 for (octave_idx_type i = 0; i < len; i++)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	738 dest[data[i]] = src[i];
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	739 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	740 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	741
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	742 case class_mask:
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	743 {
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	744 idx_mask_rep * r = dynamic_cast<idx_mask_rep *> (rep);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	745 const bool *data = r->get_data ();
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	746 octave_idx_type ext = r->extent (0);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	747 for (octave_idx_type i = 0; i < ext; i++)
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	748 if (data[i]) dest[i] = *src++;
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	749 }
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	750 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	751
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	752 default:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	753 assert (false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	754 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	755 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	756
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	757 return len;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	758 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	759
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	760 // Slice fill with specializations. No checking of bounds!
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	761 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	762 // This is equivalent to the following loop (but much faster):
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	763 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	764 // 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	765 // dest[idx(i)] = val;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	766 // return i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	767 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	768 template <class T>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	769 octave_idx_type
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	770 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	771 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	772 octave_idx_type len = rep->length (n);
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	773
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	774 switch (rep->idx_class ())
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	775 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	776 case class_colon:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	777 std::fill (dest, dest + len, val);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	778 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	779
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	780 case class_range:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	781 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	782 idx_range_rep * r = dynamic_cast<idx_range_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	783 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	784 T *sdest = dest + start;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	785 if (step == 1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	786 std::fill (sdest, sdest + len, val);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	787 else if (step == -1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	788 std::fill (sdest - len + 1, sdest + 1, val);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	789 else
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	790 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	791 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	792 sdest[j] = val;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	793 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	794 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	795 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	796
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	797 case class_scalar:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	798 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	799 idx_scalar_rep * r = dynamic_cast<idx_scalar_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	800 dest[r->get_data ()] = val;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	801 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	802 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	803
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	804 case class_vector:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	805 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	806 idx_vector_rep * r = dynamic_cast<idx_vector_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	807 const octave_idx_type *data = r->get_data ();
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	808 for (octave_idx_type i = 0; i < len; i++)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	809 dest[data[i]] = val;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	810 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	811 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	812
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	813 case class_mask:
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	814 {
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	815 idx_mask_rep * r = dynamic_cast<idx_mask_rep *> (rep);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	816 const bool *data = r->get_data ();
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	817 octave_idx_type ext = r->extent (0);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	818 for (octave_idx_type i = 0; i < ext; i++)
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	819 if (data[i]) dest[i] = val;
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	820 }
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	821 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	822
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	823 default:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	824 assert (false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	825 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	826 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	827
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	828 return len;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	829 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	830
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	831 // Generic non-breakable indexed loop. The loop body should be
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	832 // encapsulated in a single functor body. This is equivalent to the
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	833 // following loop (but faster, at least for simple inlined bodies):
8934 c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	834 //
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	835 // for (octave_idx_type i = 0; i < idx->length (n); i++) body (idx(i));
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	836
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	837 template <class Functor>
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	838 void
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	839 loop (octave_idx_type n, Functor body) const
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	840 {
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	841 octave_idx_type len = rep->length (n);
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	842
8934 c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	843 switch (rep->idx_class ())
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	844 {
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	845 case class_colon:
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	846 for (octave_idx_type i = 0; i < len; i++) body (i);
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	847 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	848
8934 c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	849 case class_range:
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	850 {
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	851 idx_range_rep * r = dynamic_cast<idx_range_rep *> (rep);
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	852 octave_idx_type start = r->get_start (), step = r->get_step ();
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	853 octave_idx_type i, j;
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	854 if (step == 1)
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	855 for (i = start, j = start + len; i < j; i++) body (i);
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	856 else if (step == -1)
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	857 for (i = start, j = start - len; i > j; i--) body (i);
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	858 else
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	859 for (i = 0, j = start; i < len; i++, j += step) body (j);
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	860 }
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	861 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	862
8934 c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	863 case class_scalar:
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	864 {
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	865 idx_scalar_rep * r = dynamic_cast<idx_scalar_rep *> (rep);
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	866 body (r->get_data ());
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	867 }
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	868 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	869
8934 c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	870 case class_vector:
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	871 {
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	872 idx_vector_rep * r = dynamic_cast<idx_vector_rep *> (rep);
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	873 const octave_idx_type *data = r->get_data ();
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	874 for (octave_idx_type i = 0; i < len; i++) body (data[i]);
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	875 }
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	876 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	877
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	878 case class_mask:
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	879 {
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	880 idx_mask_rep * r = dynamic_cast<idx_mask_rep *> (rep);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	881 const bool *data = r->get_data ();
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	882 octave_idx_type ext = r->extent (0);
9985 43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	883 for (octave_idx_type i = 0; i < ext; i++)
43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	884 if (data[i]) body (i);
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	885 }
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	886 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	887
8934 c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	888 default:
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	889 assert (false);
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	890 break;
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	891 }
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	892
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	893 }
c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	894
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	895 // Generic breakable indexed loop. The loop body should be
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	896 // encapsulated in a single functor body. This is equivalent to the
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	897 // following loop (but faster, at least for simple inlined bodies):
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	898 //
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	899 // for (octave_idx_type i = 0; i < idx->length (n); i++)
8934 c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	900 // if (body (idx(i))) break;
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	901 // return i;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	902 //
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	903
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	904 template <class Functor>
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	905 octave_idx_type
8934 c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	906 bloop (octave_idx_type n, Functor body) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	907 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	908 octave_idx_type len = rep->length (n), ret;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	909
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	910 switch (rep->idx_class ())
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	911 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	912 case class_colon:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	913 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	914 octave_idx_type i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	915 for (i = 0; i < len && body (i); i++) ;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	916 ret = i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	917 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	918 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	919
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	920 case class_range:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	921 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	922 idx_range_rep * r = dynamic_cast<idx_range_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	923 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	924 octave_idx_type i, j;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	925 if (step == 1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	926 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	927 else if (step == -1)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	928 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	929 else
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	930 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	931 ret = i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	932 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	933 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	934
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	935 case class_scalar:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	936 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	937 idx_scalar_rep * r = dynamic_cast<idx_scalar_rep *> (rep);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	938 ret = body (r->get_data ()) ? 1 : 0;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	939 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	940 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	941
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	942 case class_vector:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	943 {
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	944 idx_vector_rep * r = dynamic_cast<idx_vector_rep *> (rep);
8934 c2099a4d12ea partially optimize accumarray Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	945 const octave_idx_type *data = r->get_data ();
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	946 octave_idx_type i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	947 for (i = 0; i < len && body (data[i]); i++) ;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	948 ret = i;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	949 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	950 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	951
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	952 case class_mask:
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	953 {
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	954 idx_mask_rep * r = dynamic_cast<idx_mask_rep *> (rep);
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	955 const bool *data = r->get_data ();
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	956 octave_idx_type ext = r->extent (0), j = 0;
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	957 for (octave_idx_type i = 0; i < ext; i++)
9985 43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	958 {
43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	959 if (data[i])
43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	960 {
43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	961 if (body (i))
43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	962 break;
43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	963 else
43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	964 j++;
43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	965 }
43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	966 }
43a29eeda994 fix idx_vector::loop for masks Jaroslav Hajek <highegg@gmail.com> parents: 9894 diff changeset	967
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	968 ret = j;
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	969 }
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	970 break;
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	971
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	972 default:
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	973 assert (false);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	974 break;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	975 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	976
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	977 return ret;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	978 }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	979
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	980 // Rationale:
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	981 // This method is the key to "smart indexing". When indexing cartesian
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	982 // arrays, sometimes consecutive index vectors can be reduced into a
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	983 // single index. If rows (A) = k and i.maybe_reduce (j) gives k, then
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	984 // A(i,j)(:) is equal to A(k)(:).
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	985
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	986 // 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	987 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	988 octave_idx_type nj);
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	989
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	990 bool is_cont_range (octave_idx_type n,
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	991 octave_idx_type& l, octave_idx_type& u) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	992
9725 aea3a3a950e1 implement nth_element Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	993 // Returns the increment for ranges and colon, 0 for scalars and empty
aea3a3a950e1 implement nth_element Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	994 // vectors, 1st difference otherwise.
aea3a3a950e1 implement nth_element Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	995 octave_idx_type increment (void) const;
aea3a3a950e1 implement nth_element Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	996
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	997 idx_vector
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	998 complement (octave_idx_type n) const;
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	999
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8333 diff changeset	1000 bool is_permutation (octave_idx_type n) const;
445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 8333 diff changeset	1001
10512 aac9f4265048 rewrite sparse indexed assignment Jaroslav Hajek <highegg@gmail.com> parents: 10489 diff changeset	1002 // Returns the inverse permutation. If this is not a permutation on 1:n, the
aac9f4265048 rewrite sparse indexed assignment Jaroslav Hajek <highegg@gmail.com> parents: 10489 diff changeset	1003 // result is undefined (but no error unless extent () != n).
aac9f4265048 rewrite sparse indexed assignment Jaroslav Hajek <highegg@gmail.com> parents: 10489 diff changeset	1004 idx_vector inverse_permutation (octave_idx_type n) const;
aac9f4265048 rewrite sparse indexed assignment Jaroslav Hajek <highegg@gmail.com> parents: 10489 diff changeset	1005
9479 d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	1006 // Copies all the indices to a given array. Not allowed for colons.
d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	1007 void copy_data (octave_idx_type *data) const;
d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	1008
9879 034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	1009 // If the index is a mask, convert it to index vector.
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	1010 idx_vector unmask (void) const;
034677ab6865 smarter treatment of mask indexing Jaroslav Hajek <highegg@gmail.com> parents: 9780 diff changeset	1011
9894 83bd7f34f9da improve idx_vector->octave_value conversion Jaroslav Hajek <highegg@gmail.com> parents: 9886 diff changeset	1012 // Unconverts the index to a scalar, Range, double array or a mask.
9479 d9716e3ee0dd supply optimized compiled sub2ind & ind2sub Jaroslav Hajek <highegg@gmail.com> parents: 8972 diff changeset	1013 void unconvert (idx_class_type& iclass,
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	1014 double& scalar, Range& range,
9894 83bd7f34f9da improve idx_vector->octave_value conversion Jaroslav Hajek <highegg@gmail.com> parents: 9886 diff changeset	1015 Array<double>& array, Array<bool>& mask) const;
10339 de2d43bcb083 optimize some lazy index operations Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	1016
de2d43bcb083 optimize some lazy index operations Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	1017 Array<octave_idx_type> as_array (void) const;
10475 739109ef88ee helper idx_vector methods Jaroslav Hajek <highegg@gmail.com> parents: 10425 diff changeset	1018
10487 942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	1019 // Raw pointer to index array. This is non-const because it may be
942386d6d1a5 style fixes John W. Eaton <jwe@octave.org> parents: 10486 diff changeset	1020 // necessary to mutate the index.
10475 739109ef88ee helper idx_vector methods Jaroslav Hajek <highegg@gmail.com> parents: 10425 diff changeset	1021 const octave_idx_type *raw (void);
10516 f0266ee4aabe optimize some special indexing & assignment cases Jaroslav Hajek <highegg@gmail.com> parents: 10512 diff changeset	1022
f0266ee4aabe optimize some special indexing & assignment cases Jaroslav Hajek <highegg@gmail.com> parents: 10512 diff changeset	1023 bool is_vector (void) const;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	1024
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	1025 // FIXME -- these are here for compatibility. They should be removed
9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	1026 // when no longer in use.
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1027
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	1028 octave_idx_type elem (octave_idx_type n) const
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1029 { return (*this) (n); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1030
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1031 bool is_colon_equiv (octave_idx_type n, int) const
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1032 { return is_colon_equiv (n); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1033
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	1034 octave_idx_type
9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	1035 freeze (octave_idx_type z_len, const char *tag, bool resize_ok = false);
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1036
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1037 void sort (bool uniq = false)
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1038 { *this = sorted (uniq); }
7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1039
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	1040 octave_idx_type ones_count (void) const;
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1041
8333 9238637cb81c style fixes John W. Eaton <jwe@octave.org> parents: 8319 diff changeset	1042 octave_idx_type max (void) const { return extent (1) - 1; }
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	1043
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	1044 private:
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	1045
8290 7cbe01c21986 improve dense array indexing Jaroslav Hajek <highegg@gmail.com> parents: 7573 diff changeset	1046 idx_base_rep *rep;
1 78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	1047
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	1048 };
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	1049
78fd87e624cb [project @ 1993-08-08 01:13:40 by jwe] jwe parents: diff changeset	1050 #endif

Mercurial > octave

annotate liboctave/idx-vector.h @ 13735:2cd903770748