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