| 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 */
|
|---|