100db7afdSDavid E. O'Brien // The template and inlines for the -*- C++ -*- gslice class.
200db7afdSDavid E. O'Brien 
3*f8a1b7d9SAlexander Kabaev // Copyright (C) 1997, 1998, 1999, 2000, 2001, 2004, 2005
4ffeaf689SAlexander Kabaev // Free Software Foundation, Inc.
500db7afdSDavid E. O'Brien //
600db7afdSDavid E. O'Brien // This file is part of the GNU ISO C++ Library.  This library is free
700db7afdSDavid E. O'Brien // software; you can redistribute it and/or modify it under the
800db7afdSDavid E. O'Brien // terms of the GNU General Public License as published by the
900db7afdSDavid E. O'Brien // Free Software Foundation; either version 2, or (at your option)
1000db7afdSDavid E. O'Brien // any later version.
1100db7afdSDavid E. O'Brien 
1200db7afdSDavid E. O'Brien // This library is distributed in the hope that it will be useful,
1300db7afdSDavid E. O'Brien // but WITHOUT ANY WARRANTY; without even the implied warranty of
1400db7afdSDavid E. O'Brien // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
1500db7afdSDavid E. O'Brien // GNU General Public License for more details.
1600db7afdSDavid E. O'Brien 
1700db7afdSDavid E. O'Brien // You should have received a copy of the GNU General Public License along
1800db7afdSDavid E. O'Brien // with this library; see the file COPYING.  If not, write to the Free
19*f8a1b7d9SAlexander Kabaev // Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301,
2000db7afdSDavid E. O'Brien // USA.
2100db7afdSDavid E. O'Brien 
2200db7afdSDavid E. O'Brien // As a special exception, you may use this file as part of a free software
2300db7afdSDavid E. O'Brien // library without restriction.  Specifically, if other files instantiate
2400db7afdSDavid E. O'Brien // templates or use macros or inline functions from this file, or you compile
2500db7afdSDavid E. O'Brien // this file and link it with other files to produce an executable, this
2600db7afdSDavid E. O'Brien // file does not by itself cause the resulting executable to be covered by
2700db7afdSDavid E. O'Brien // the GNU General Public License.  This exception does not however
2800db7afdSDavid E. O'Brien // invalidate any other reasons why the executable file might be covered by
2900db7afdSDavid E. O'Brien // the GNU General Public License.
3000db7afdSDavid E. O'Brien 
3100db7afdSDavid E. O'Brien /** @file gslice.h
3200db7afdSDavid E. O'Brien  *  This is an internal header file, included by other library headers.
3300db7afdSDavid E. O'Brien  *  You should not attempt to use it directly.
3400db7afdSDavid E. O'Brien  */
3500db7afdSDavid E. O'Brien 
36*f8a1b7d9SAlexander Kabaev // Written by Gabriel Dos Reis <[email protected]>
37*f8a1b7d9SAlexander Kabaev 
38ffeaf689SAlexander Kabaev #ifndef _GSLICE_H
39ffeaf689SAlexander Kabaev #define _GSLICE_H 1
4000db7afdSDavid E. O'Brien 
4100db7afdSDavid E. O'Brien #pragma GCC system_header
4200db7afdSDavid E. O'Brien 
_GLIBCXX_BEGIN_NAMESPACE(std)43*f8a1b7d9SAlexander Kabaev _GLIBCXX_BEGIN_NAMESPACE(std)
4400db7afdSDavid E. O'Brien 
45ffeaf689SAlexander Kabaev   /**
46ffeaf689SAlexander Kabaev    *  @brief  Class defining multi-dimensional subset of an array.
47ffeaf689SAlexander Kabaev    *
48ffeaf689SAlexander Kabaev    *  The slice class represents a multi-dimensional subset of an array,
49ffeaf689SAlexander Kabaev    *  specified by three parameter sets: start offset, size array, and stride
50ffeaf689SAlexander Kabaev    *  array.  The start offset is the index of the first element of the array
51ffeaf689SAlexander Kabaev    *  that is part of the subset.  The size and stride array describe each
52ffeaf689SAlexander Kabaev    *  dimension of the slice.  Size is the number of elements in that
53ffeaf689SAlexander Kabaev    *  dimension, and stride is the distance in the array between successive
54ffeaf689SAlexander Kabaev    *  elements in that dimension.  Each dimension's size and stride is taken
55ffeaf689SAlexander Kabaev    *  to begin at an array element described by the previous dimension.  The
56ffeaf689SAlexander Kabaev    *  size array and stride array must be the same size.
57ffeaf689SAlexander Kabaev    *
58ffeaf689SAlexander Kabaev    *  For example, if you have offset==3, stride[0]==11, size[1]==3,
59ffeaf689SAlexander Kabaev    *  stride[1]==3, then slice[0,0]==array[3], slice[0,1]==array[6],
60ffeaf689SAlexander Kabaev    *  slice[0,2]==array[9], slice[1,0]==array[14], slice[1,1]==array[17],
61ffeaf689SAlexander Kabaev    *  slice[1,2]==array[20].
62ffeaf689SAlexander Kabaev    */
6300db7afdSDavid E. O'Brien   class gslice
6400db7afdSDavid E. O'Brien   {
6500db7afdSDavid E. O'Brien   public:
66ffeaf689SAlexander Kabaev     ///  Construct an empty slice.
6700db7afdSDavid E. O'Brien     gslice ();
68ffeaf689SAlexander Kabaev 
69ffeaf689SAlexander Kabaev     /**
70ffeaf689SAlexander Kabaev      *  @brief  Construct a slice.
71ffeaf689SAlexander Kabaev      *
72ffeaf689SAlexander Kabaev      *  Constructs a slice with as many dimensions as the length of the @a l
73ffeaf689SAlexander Kabaev      *  and @a s arrays.
74ffeaf689SAlexander Kabaev      *
75ffeaf689SAlexander Kabaev      *  @param  o  Offset in array of first element.
76ffeaf689SAlexander Kabaev      *  @param  l  Array of dimension lengths.
77ffeaf689SAlexander Kabaev      *  @param  s  Array of dimension strides between array elements.
78ffeaf689SAlexander Kabaev      */
7900db7afdSDavid E. O'Brien     gslice(size_t, const valarray<size_t>&, const valarray<size_t>&);
80ffeaf689SAlexander Kabaev 
8100db7afdSDavid E. O'Brien     // XXX: the IS says the copy-ctor and copy-assignment operators are
8200db7afdSDavid E. O'Brien     //      synthetized by the compiler but they are just unsuitable
8300db7afdSDavid E. O'Brien     //      for a ref-counted semantic
84ffeaf689SAlexander Kabaev     ///  Copy constructor.
8500db7afdSDavid E. O'Brien     gslice(const gslice&);
86ffeaf689SAlexander Kabaev 
87ffeaf689SAlexander Kabaev     ///  Destructor.
8800db7afdSDavid E. O'Brien     ~gslice();
8900db7afdSDavid E. O'Brien 
9000db7afdSDavid E. O'Brien     // XXX: See the note above.
91ffeaf689SAlexander Kabaev     ///  Assignment operator.
9200db7afdSDavid E. O'Brien     gslice& operator=(const gslice&);
9300db7afdSDavid E. O'Brien 
94ffeaf689SAlexander Kabaev     ///  Return array offset of first slice element.
9500db7afdSDavid E. O'Brien     size_t           start() const;
96ffeaf689SAlexander Kabaev 
97ffeaf689SAlexander Kabaev     ///  Return array of sizes of slice dimensions.
9800db7afdSDavid E. O'Brien     valarray<size_t> size() const;
99ffeaf689SAlexander Kabaev 
100ffeaf689SAlexander Kabaev     ///  Return array of array strides for each dimension.
10100db7afdSDavid E. O'Brien     valarray<size_t> stride() const;
10200db7afdSDavid E. O'Brien 
10300db7afdSDavid E. O'Brien   private:
104*f8a1b7d9SAlexander Kabaev     struct _Indexer
105*f8a1b7d9SAlexander Kabaev     {
10600db7afdSDavid E. O'Brien       size_t _M_count;
10700db7afdSDavid E. O'Brien       size_t _M_start;
10800db7afdSDavid E. O'Brien       valarray<size_t> _M_size;
10900db7afdSDavid E. O'Brien       valarray<size_t> _M_stride;
110ffeaf689SAlexander Kabaev       valarray<size_t> _M_index; // Linear array of referenced indices
11100db7afdSDavid E. O'Brien       _Indexer(size_t, const valarray<size_t>&,
11200db7afdSDavid E. O'Brien 	       const valarray<size_t>&);
113*f8a1b7d9SAlexander Kabaev       void
114*f8a1b7d9SAlexander Kabaev       _M_increment_use()
115*f8a1b7d9SAlexander Kabaev       { ++_M_count; }
116*f8a1b7d9SAlexander Kabaev 
117*f8a1b7d9SAlexander Kabaev       size_t
118*f8a1b7d9SAlexander Kabaev       _M_decrement_use()
119*f8a1b7d9SAlexander Kabaev       { return --_M_count; }
12000db7afdSDavid E. O'Brien     };
12100db7afdSDavid E. O'Brien 
12200db7afdSDavid E. O'Brien     _Indexer* _M_index;
12300db7afdSDavid E. O'Brien 
12400db7afdSDavid E. O'Brien     template<typename _Tp> friend class valarray;
12500db7afdSDavid E. O'Brien   };
12600db7afdSDavid E. O'Brien 
12700db7afdSDavid E. O'Brien   inline size_t
start()12800db7afdSDavid E. O'Brien   gslice::start () const
12900db7afdSDavid E. O'Brien   { return _M_index ? _M_index->_M_start : 0; }
13000db7afdSDavid E. O'Brien 
13100db7afdSDavid E. O'Brien   inline valarray<size_t>
size()13200db7afdSDavid E. O'Brien   gslice::size () const
13300db7afdSDavid E. O'Brien   { return _M_index ? _M_index->_M_size : valarray<size_t>(); }
13400db7afdSDavid E. O'Brien 
13500db7afdSDavid E. O'Brien   inline valarray<size_t>
stride()13600db7afdSDavid E. O'Brien   gslice::stride () const
13700db7afdSDavid E. O'Brien   { return _M_index ? _M_index->_M_stride : valarray<size_t>(); }
13800db7afdSDavid E. O'Brien 
gslice()13900db7afdSDavid E. O'Brien   inline gslice::gslice () : _M_index(0) {}
14000db7afdSDavid E. O'Brien 
14100db7afdSDavid E. O'Brien   inline
gslice(size_t __o,const valarray<size_t> & __l,const valarray<size_t> & __s)14200db7afdSDavid E. O'Brien   gslice::gslice(size_t __o, const valarray<size_t>& __l,
14300db7afdSDavid E. O'Brien 		 const valarray<size_t>& __s)
14400db7afdSDavid E. O'Brien   : _M_index(new gslice::_Indexer(__o, __l, __s)) {}
14500db7afdSDavid E. O'Brien 
14600db7afdSDavid E. O'Brien   inline
gslice(const gslice & __g)14700db7afdSDavid E. O'Brien   gslice::gslice(const gslice& __g) : _M_index(__g._M_index)
14800db7afdSDavid E. O'Brien   { if (_M_index) _M_index->_M_increment_use(); }
14900db7afdSDavid E. O'Brien 
15000db7afdSDavid E. O'Brien   inline
~gslice()15100db7afdSDavid E. O'Brien   gslice::~gslice()
152*f8a1b7d9SAlexander Kabaev   {
153*f8a1b7d9SAlexander Kabaev     if (_M_index && _M_index->_M_decrement_use() == 0)
154*f8a1b7d9SAlexander Kabaev       delete _M_index;
155*f8a1b7d9SAlexander Kabaev   }
15600db7afdSDavid E. O'Brien 
15700db7afdSDavid E. O'Brien   inline gslice&
15800db7afdSDavid E. O'Brien   gslice::operator= (const gslice& __g)
15900db7afdSDavid E. O'Brien   {
160*f8a1b7d9SAlexander Kabaev     if (__g._M_index)
161*f8a1b7d9SAlexander Kabaev       __g._M_index->_M_increment_use();
162*f8a1b7d9SAlexander Kabaev     if (_M_index && _M_index->_M_decrement_use() == 0)
163*f8a1b7d9SAlexander Kabaev       delete _M_index;
16400db7afdSDavid E. O'Brien     _M_index = __g._M_index;
16500db7afdSDavid E. O'Brien     return *this;
16600db7afdSDavid E. O'Brien   }
16700db7afdSDavid E. O'Brien 
168*f8a1b7d9SAlexander Kabaev _GLIBCXX_END_NAMESPACE
16900db7afdSDavid E. O'Brien 
170ffeaf689SAlexander Kabaev #endif /* _GSLICE_H */
171