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