| 1 | // The template and inlines for the -*- C++ -*- gslice class. | 
 
 
 
 
 | 2 |  | 
 
 
 
 
 | 3 | // Copyright (C) 1997-2021 Free Software Foundation, Inc. | 
 
 
 
 
 | 4 | // | 
 
 
 
 
 | 5 | // This file is part of the GNU ISO C++ Library.  This library is free | 
 
 
 
 
 | 6 | // software; you can redistribute it and/or modify it under the | 
 
 
 
 
 | 7 | // terms of the GNU General Public License as published by the | 
 
 
 
 
 | 8 | // Free Software Foundation; either version 3, or (at your option) | 
 
 
 
 
 | 9 | // any later version. | 
 
 
 
 
 | 10 |  | 
 
 
 
 
 | 11 | // This library is distributed in the hope that it will be useful, | 
 
 
 
 
 | 12 | // but WITHOUT ANY WARRANTY; without even the implied warranty of | 
 
 
 
 
 | 13 | // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the | 
 
 
 
 
 | 14 | // GNU General Public License for more details. | 
 
 
 
 
 | 15 |  | 
 
 
 
 
 | 16 | // Under Section 7 of GPL version 3, you are granted additional | 
 
 
 
 
 | 17 | // permissions described in the GCC Runtime Library Exception, version | 
 
 
 
 
 | 18 | // 3.1, as published by the Free Software Foundation. | 
 
 
 
 
 | 19 |  | 
 
 
 
 
 | 20 | // You should have received a copy of the GNU General Public License and | 
 
 
 
 
 | 21 | // a copy of the GCC Runtime Library Exception along with this program; | 
 
 
 
 
 | 22 | // see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see | 
 
 
 
 
 | 23 | // <http://www.gnu.org/licenses/>. | 
 
 
 
 
 | 24 |  | 
 
 
 
 
 | 25 | /** @file bits/gslice.h | 
 
 
 
 
 | 26 | *  This is an internal header file, included by other library headers. | 
 
 
 
 
 | 27 | *  Do not attempt to use it directly. @headername{valarray} | 
 
 
 
 
 | 28 | */ | 
 
 
 
 
 | 29 |  | 
 
 
 
 
 | 30 | // Written by Gabriel Dos Reis <Gabriel.Dos-Reis@DPTMaths.ENS-Cachan.Fr> | 
 
 
 
 
 | 31 |  | 
 
 
 
 
 | 32 | #ifndef _GSLICE_H | 
 
 
 
 
 | 33 | #define _GSLICE_H 1 | 
 
 
 
 
 | 34 |  | 
 
 
 
 
 | 35 | #pragma GCC system_header | 
 
 
 
 
 | 36 |  | 
 
 
 
 
 | 37 | namespace std _GLIBCXX_VISIBILITY(default) | 
 
 
 
 
 | 38 | { | 
 
 
 
 
 | 39 | _GLIBCXX_BEGIN_NAMESPACE_VERSION | 
 
 
 
 
 | 40 |  | 
 
 
 
 
 | 41 | /** | 
 
 
 
 
 | 42 | * @addtogroup numeric_arrays | 
 
 
 
 
 | 43 | * @{ | 
 
 
 
 
 | 44 | */ | 
 
 
 
 
 | 45 |  | 
 
 
 
 
 | 46 | /** | 
 
 
 
 
 | 47 | *  @brief  Class defining multi-dimensional subset of an array. | 
 
 
 
 
 | 48 | * | 
 
 
 
 
 | 49 | *  The slice class represents a multi-dimensional subset of an array, | 
 
 
 
 
 | 50 | *  specified by three parameter sets: start offset, size array, and stride | 
 
 
 
 
 | 51 | *  array.  The start offset is the index of the first element of the array | 
 
 
 
 
 | 52 | *  that is part of the subset.  The size and stride array describe each | 
 
 
 
 
 | 53 | *  dimension of the slice.  Size is the number of elements in that | 
 
 
 
 
 | 54 | *  dimension, and stride is the distance in the array between successive | 
 
 
 
 
 | 55 | *  elements in that dimension.  Each dimension's size and stride is taken | 
 
 
 
 
 | 56 | *  to begin at an array element described by the previous dimension.  The | 
 
 
 
 
 | 57 | *  size array and stride array must be the same size. | 
 
 
 
 
 | 58 | * | 
 
 
 
 
 | 59 | *  For example, if you have offset==3, stride[0]==11, size[1]==3, | 
 
 
 
 
 | 60 | *  stride[1]==3, then slice[0,0]==array[3], slice[0,1]==array[6], | 
 
 
 
 
 | 61 | *  slice[0,2]==array[9], slice[1,0]==array[14], slice[1,1]==array[17], | 
 
 
 
 
 | 62 | *  slice[1,2]==array[20]. | 
 
 
 
 
 | 63 | */ | 
 
 
 
 
 | 64 | class gslice | 
 
 
 
 
 | 65 | { | 
 
 
 
 
 | 66 | public: | 
 
 
 
 
 | 67 | ///  Construct an empty slice. | 
 
 
 
 
 | 68 | gslice(); | 
 
 
 
 
 | 69 |  | 
 
 
 
 
 | 70 | /** | 
 
 
 
 
 | 71 | *  @brief  Construct a slice. | 
 
 
 
 
 | 72 | * | 
 
 
 
 
 | 73 | *  Constructs a slice with as many dimensions as the length of the @a l | 
 
 
 
 
 | 74 | *  and @a s arrays. | 
 
 
 
 
 | 75 | * | 
 
 
 
 
 | 76 | *  @param  __o  Offset in array of first element. | 
 
 
 
 
 | 77 | *  @param  __l  Array of dimension lengths. | 
 
 
 
 
 | 78 | *  @param  __s  Array of dimension strides between array elements. | 
 
 
 
 
 | 79 | */ | 
 
 
 
 
 | 80 | gslice(size_t __o, const valarray<size_t>& __l, | 
 
 
 
 
 | 81 | const valarray<size_t>& __s); | 
 
 
 
 
 | 82 |  | 
 
 
 
 
 | 83 | // XXX: the IS says the copy-ctor and copy-assignment operators are | 
 
 
 
 
 | 84 | //      synthesized by the compiler but they are just unsuitable | 
 
 
 
 
 | 85 | //      for a ref-counted semantic | 
 
 
 
 
 | 86 | ///  Copy constructor. | 
 
 
 
 
 | 87 | gslice(const gslice&); | 
 
 
 
 
 | 88 |  | 
 
 
 
 
 | 89 | ///  Destructor. | 
 
 
 
 
 | 90 | ~gslice(); | 
 
 
 
 
 | 91 |  | 
 
 
 
 
 | 92 | // XXX: See the note above. | 
 
 
 
 
 | 93 | ///  Assignment operator. | 
 
 
 
 
 | 94 | gslice& operator=(const gslice&); | 
 
 
 
 
 | 95 |  | 
 
 
 
 
 | 96 | ///  Return array offset of first slice element. | 
 
 
 
 
 | 97 | size_t           start() const; | 
 
 
 
 
 | 98 |  | 
 
 
 
 
 | 99 | ///  Return array of sizes of slice dimensions. | 
 
 
 
 
 | 100 | valarray<size_t> size() const; | 
 
 
 
 
 | 101 |  | 
 
 
 
 
 | 102 | ///  Return array of array strides for each dimension. | 
 
 
 
 
 | 103 | valarray<size_t> stride() const; | 
 
 
 
 
 | 104 |  | 
 
 
 
 
 | 105 | private: | 
 
 
 
 
 | 106 | struct _Indexer | 
 
 
 
 
 | 107 | { | 
 
 
 
 
 | 108 | size_t _M_count; | 
 
 
 
 
 | 109 | size_t _M_start; | 
 
 
 
 
 | 110 | valarray<size_t> _M_size; | 
 
 
 
 
 | 111 | valarray<size_t> _M_stride; | 
 
 
 
 
 | 112 | valarray<size_t> _M_index; // Linear array of referenced indices | 
 
 
 
 
 | 113 |  | 
 
 
 
 
 | 114 | _Indexer() | 
 
 
 
 
 | 115 | : _M_count(1), _M_start(0), _M_size(), _M_stride(), _M_index() {} | 
 
 
 
 
 | 116 |  | 
 
 
 
 
 | 117 | _Indexer(size_t, const valarray<size_t>&, | 
 
 
 
 
 | 118 | const valarray<size_t>&); | 
 
 
 
 
 | 119 |  | 
 
 
 
 
 | 120 | void | 
 
 
 
 
 | 121 | _M_increment_use() | 
 
 
 
 
 | 122 | { ++_M_count; } | 
 
 
 
 
 | 123 |  | 
 
 
 
 
 | 124 | size_t | 
 
 
 
 
 | 125 | _M_decrement_use() | 
 
 
 
 
 | 126 | { return --_M_count; } | 
 
 
 
 
 | 127 | }; | 
 
 
 
 
 | 128 |  | 
 
 
 
 
 | 129 | _Indexer* _M_index; | 
 
 
 
 
 | 130 |  | 
 
 
 
 
 | 131 | template<typename _Tp> friend class valarray; | 
 
 
 
 
 | 132 | }; | 
 
 
 
 
 | 133 |  | 
 
 
 
 
 | 134 | inline size_t | 
 
 
 
 
 | 135 | gslice::start() const | 
 
 
 
 
 | 136 | { return _M_index ? _M_index->_M_start : 0; } | 
 
 
 
 
 | 137 |  | 
 
 
 
 
 | 138 | inline valarray<size_t> | 
 
 
 
 
 | 139 | gslice::size() const | 
 
 
 
 
 | 140 | { return _M_index ? _M_index->_M_size : valarray<size_t>(); } | 
 
 
 
 
 | 141 |  | 
 
 
 
 
 | 142 | inline valarray<size_t> | 
 
 
 
 
 | 143 | gslice::stride() const | 
 
 
 
 
 | 144 | { return _M_index ? _M_index->_M_stride : valarray<size_t>(); } | 
 
 
 
 
 | 145 |  | 
 
 
 
 
 | 146 | // _GLIBCXX_RESOLVE_LIB_DEFECTS | 
 
 
 
 
 | 147 | // 543. valarray slice default constructor | 
 
 
 
 
 | 148 | inline | 
 
 
 
 
 | 149 | gslice::gslice() | 
 
 
 
 
 | 150 | : _M_index(new gslice::_Indexer()) {} | 
 
 
 
 
 | 151 |  | 
 
 
 
 
 | 152 | inline | 
 
 
 
 
 | 153 | gslice::gslice(size_t __o, const valarray<size_t>& __l, | 
 
 
 
 
 | 154 | const valarray<size_t>& __s) | 
 
 
 
 
 | 155 | : _M_index(new gslice::_Indexer(__o, __l, __s)) {} | 
 
 
 
 
 | 156 |  | 
 
 
 
 
 | 157 | inline | 
 
 
 
 
 | 158 | gslice::gslice(const gslice& __g) | 
 
 
 
 
 | 159 | : _M_index(__g._M_index) | 
 
 
 
 
 | 160 | { if (_M_index) _M_index->_M_increment_use(); } | 
 
 
 
 
 | 161 |  | 
 
 
 
 
 | 162 | inline | 
 
 
 
 
 | 163 | gslice::~gslice() | 
 
 
 
 
 | 164 | { | 
 
 
 
 
 | 165 | if (_M_index && _M_index->_M_decrement_use() == 0) | 
 
 
 
 
 | 166 | delete _M_index; | 
 
 
 
 
 | 167 | } | 
 
 
 
 
 | 168 |  | 
 
 
 
 
 | 169 | inline gslice& | 
 
 
 
 
 | 170 | gslice::operator=(const gslice& __g) | 
 
 
 
 
 | 171 | { | 
 
 
 
 
 | 172 | if (__g._M_index) | 
 
 
 
 
 | 173 | __g._M_index->_M_increment_use(); | 
 
 
 
 
 | 174 | if (_M_index && _M_index->_M_decrement_use() == 0) | 
 
 
 
 
 | 175 | delete _M_index; | 
 
 
 
 
 | 176 | _M_index = __g._M_index; | 
 
 
 
 
 | 177 | return *this; | 
 
 
 
 
 | 178 | } | 
 
 
 
 
 | 179 |  | 
 
 
 
 
 | 180 | /// @} group numeric_arrays | 
 
 
 
 
 | 181 |  | 
 
 
 
 
 | 182 | _GLIBCXX_END_NAMESPACE_VERSION | 
 
 
 
 
 | 183 | } // namespace | 
 
 
 
 
 | 184 |  | 
 
 
 
 
 | 185 | #endif /* _GSLICE_H */ |