1*e4b17023SJohn Marino // The template and inlines for the -*- C++ -*- slice_array class.
2*e4b17023SJohn Marino 
3*e4b17023SJohn Marino // Copyright (C) 1997, 1998, 1999, 2000, 2001, 2002, 2004, 2005, 2006, 2009,
4*e4b17023SJohn Marino // 2010, 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/slice_array.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 _SLICE_ARRAY_H
34*e4b17023SJohn Marino #define _SLICE_ARRAY_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 one-dimensional subset of an array.
49*e4b17023SJohn Marino    *
50*e4b17023SJohn Marino    *  The slice class represents a one-dimensional subset of an array,
51*e4b17023SJohn Marino    *  specified by three parameters: start offset, size, and stride.  The
52*e4b17023SJohn Marino    *  start offset is the index of the first element of the array that is part
53*e4b17023SJohn Marino    *  of the subset.  The size is the total number of elements in the subset.
54*e4b17023SJohn Marino    *  Stride is the distance between each successive array element to include
55*e4b17023SJohn Marino    *  in the subset.
56*e4b17023SJohn Marino    *
57*e4b17023SJohn Marino    *  For example, with an array of size 10, and a slice with offset 1, size 3
58*e4b17023SJohn Marino    *  and stride 2, the subset consists of array elements 1, 3, and 5.
59*e4b17023SJohn Marino    */
60*e4b17023SJohn Marino   class slice
61*e4b17023SJohn Marino   {
62*e4b17023SJohn Marino   public:
63*e4b17023SJohn Marino     ///  Construct an empty slice.
64*e4b17023SJohn Marino     slice();
65*e4b17023SJohn Marino 
66*e4b17023SJohn Marino     /**
67*e4b17023SJohn Marino      *  @brief  Construct a slice.
68*e4b17023SJohn Marino      *
69*e4b17023SJohn Marino      *  @param  __o  Offset in array of first element.
70*e4b17023SJohn Marino      *  @param  __d  Number of elements in slice.
71*e4b17023SJohn Marino      *  @param  __s  Stride between array elements.
72*e4b17023SJohn Marino      */
73*e4b17023SJohn Marino     slice(size_t __o, size_t __d, size_t __s);
74*e4b17023SJohn Marino 
75*e4b17023SJohn Marino     ///  Return array offset of first slice element.
76*e4b17023SJohn Marino     size_t start() const;
77*e4b17023SJohn Marino     ///  Return size of slice.
78*e4b17023SJohn Marino     size_t size() const;
79*e4b17023SJohn Marino     ///  Return array stride of slice.
80*e4b17023SJohn Marino     size_t stride() const;
81*e4b17023SJohn Marino 
82*e4b17023SJohn Marino   private:
83*e4b17023SJohn Marino     size_t _M_off;                      // offset
84*e4b17023SJohn Marino     size_t _M_sz;			// size
85*e4b17023SJohn Marino     size_t _M_st;			// stride unit
86*e4b17023SJohn Marino   };
87*e4b17023SJohn Marino 
88*e4b17023SJohn Marino   // _GLIBCXX_RESOLVE_LIB_DEFECTS
89*e4b17023SJohn Marino   // 543. valarray slice default constructor
90*e4b17023SJohn Marino   inline
91*e4b17023SJohn Marino   slice::slice()
92*e4b17023SJohn Marino   : _M_off(0), _M_sz(0), _M_st(0) {}
93*e4b17023SJohn Marino 
94*e4b17023SJohn Marino   inline
95*e4b17023SJohn Marino   slice::slice(size_t __o, size_t __d, size_t __s)
96*e4b17023SJohn Marino   : _M_off(__o), _M_sz(__d), _M_st(__s) {}
97*e4b17023SJohn Marino 
98*e4b17023SJohn Marino   inline size_t
99*e4b17023SJohn Marino   slice::start() const
100*e4b17023SJohn Marino   { return _M_off; }
101*e4b17023SJohn Marino 
102*e4b17023SJohn Marino   inline size_t
103*e4b17023SJohn Marino   slice::size() const
104*e4b17023SJohn Marino   { return _M_sz; }
105*e4b17023SJohn Marino 
106*e4b17023SJohn Marino   inline size_t
107*e4b17023SJohn Marino   slice::stride() const
108*e4b17023SJohn Marino   { return _M_st; }
109*e4b17023SJohn Marino 
110*e4b17023SJohn Marino   /**
111*e4b17023SJohn Marino    *  @brief  Reference to one-dimensional subset of an array.
112*e4b17023SJohn Marino    *
113*e4b17023SJohn Marino    *  A slice_array is a reference to the actual elements of an array
114*e4b17023SJohn Marino    *  specified by a slice.  The way to get a slice_array is to call
115*e4b17023SJohn Marino    *  operator[](slice) on a valarray.  The returned slice_array then permits
116*e4b17023SJohn Marino    *  carrying operations out on the referenced subset of elements in the
117*e4b17023SJohn Marino    *  original valarray.  For example, operator+=(valarray) will add values
118*e4b17023SJohn Marino    *  to the subset of elements in the underlying valarray this slice_array
119*e4b17023SJohn Marino    *  refers to.
120*e4b17023SJohn Marino    *
121*e4b17023SJohn Marino    *  @param  Tp  Element type.
122*e4b17023SJohn Marino    */
123*e4b17023SJohn Marino   template<typename _Tp>
124*e4b17023SJohn Marino     class slice_array
125*e4b17023SJohn Marino     {
126*e4b17023SJohn Marino     public:
127*e4b17023SJohn Marino       typedef _Tp value_type;
128*e4b17023SJohn Marino 
129*e4b17023SJohn Marino       // _GLIBCXX_RESOLVE_LIB_DEFECTS
130*e4b17023SJohn Marino       // 253. valarray helper functions are almost entirely useless
131*e4b17023SJohn Marino 
132*e4b17023SJohn Marino       ///  Copy constructor.  Both slices refer to the same underlying array.
133*e4b17023SJohn Marino       slice_array(const slice_array&);
134*e4b17023SJohn Marino 
135*e4b17023SJohn Marino       ///  Assignment operator.  Assigns slice elements to corresponding
136*e4b17023SJohn Marino       ///  elements of @a a.
137*e4b17023SJohn Marino       slice_array& operator=(const slice_array&);
138*e4b17023SJohn Marino 
139*e4b17023SJohn Marino       ///  Assign slice elements to corresponding elements of @a v.
140*e4b17023SJohn Marino       void operator=(const valarray<_Tp>&) const;
141*e4b17023SJohn Marino       ///  Multiply slice elements by corresponding elements of @a v.
142*e4b17023SJohn Marino       void operator*=(const valarray<_Tp>&) const;
143*e4b17023SJohn Marino       ///  Divide slice elements by corresponding elements of @a v.
144*e4b17023SJohn Marino       void operator/=(const valarray<_Tp>&) const;
145*e4b17023SJohn Marino       ///  Modulo slice elements by corresponding elements of @a v.
146*e4b17023SJohn Marino       void operator%=(const valarray<_Tp>&) const;
147*e4b17023SJohn Marino       ///  Add corresponding elements of @a v to slice elements.
148*e4b17023SJohn Marino       void operator+=(const valarray<_Tp>&) const;
149*e4b17023SJohn Marino       ///  Subtract corresponding elements of @a v from slice elements.
150*e4b17023SJohn Marino       void operator-=(const valarray<_Tp>&) const;
151*e4b17023SJohn Marino       ///  Logical xor slice elements with corresponding elements of @a v.
152*e4b17023SJohn Marino       void operator^=(const valarray<_Tp>&) const;
153*e4b17023SJohn Marino       ///  Logical and slice elements with corresponding elements of @a v.
154*e4b17023SJohn Marino       void operator&=(const valarray<_Tp>&) const;
155*e4b17023SJohn Marino       ///  Logical or slice elements with corresponding elements of @a v.
156*e4b17023SJohn Marino       void operator|=(const valarray<_Tp>&) const;
157*e4b17023SJohn Marino       ///  Left shift slice elements by corresponding elements of @a v.
158*e4b17023SJohn Marino       void operator<<=(const valarray<_Tp>&) const;
159*e4b17023SJohn Marino       ///  Right shift slice elements by corresponding elements of @a v.
160*e4b17023SJohn Marino       void operator>>=(const valarray<_Tp>&) const;
161*e4b17023SJohn Marino       ///  Assign all slice elements to @a t.
162*e4b17023SJohn Marino       void operator=(const _Tp &) const;
163*e4b17023SJohn Marino       //        ~slice_array ();
164*e4b17023SJohn Marino 
165*e4b17023SJohn Marino       template<class _Dom>
166*e4b17023SJohn Marino         void operator=(const _Expr<_Dom, _Tp>&) const;
167*e4b17023SJohn Marino       template<class _Dom>
168*e4b17023SJohn Marino 	void operator*=(const _Expr<_Dom, _Tp>&) const;
169*e4b17023SJohn Marino       template<class _Dom>
170*e4b17023SJohn Marino 	void operator/=(const _Expr<_Dom, _Tp>&) const;
171*e4b17023SJohn Marino       template<class _Dom>
172*e4b17023SJohn Marino 	void operator%=(const _Expr<_Dom, _Tp>&) const;
173*e4b17023SJohn Marino       template<class _Dom>
174*e4b17023SJohn Marino 	void operator+=(const _Expr<_Dom, _Tp>&) const;
175*e4b17023SJohn Marino       template<class _Dom>
176*e4b17023SJohn Marino 	void operator-=(const _Expr<_Dom, _Tp>&) const;
177*e4b17023SJohn Marino       template<class _Dom>
178*e4b17023SJohn Marino 	void operator^=(const _Expr<_Dom, _Tp>&) const;
179*e4b17023SJohn Marino       template<class _Dom>
180*e4b17023SJohn Marino 	void operator&=(const _Expr<_Dom, _Tp>&) const;
181*e4b17023SJohn Marino       template<class _Dom>
182*e4b17023SJohn Marino 	void operator|=(const _Expr<_Dom, _Tp>&) const;
183*e4b17023SJohn Marino       template<class _Dom>
184*e4b17023SJohn Marino 	void operator<<=(const _Expr<_Dom, _Tp>&) const;
185*e4b17023SJohn Marino       template<class _Dom>
186*e4b17023SJohn Marino 	void operator>>=(const _Expr<_Dom, _Tp>&) const;
187*e4b17023SJohn Marino 
188*e4b17023SJohn Marino     private:
189*e4b17023SJohn Marino       friend class valarray<_Tp>;
190*e4b17023SJohn Marino       slice_array(_Array<_Tp>, const slice&);
191*e4b17023SJohn Marino 
192*e4b17023SJohn Marino       const size_t      _M_sz;
193*e4b17023SJohn Marino       const size_t      _M_stride;
194*e4b17023SJohn Marino       const _Array<_Tp> _M_array;
195*e4b17023SJohn Marino 
196*e4b17023SJohn Marino       // not implemented
197*e4b17023SJohn Marino       slice_array();
198*e4b17023SJohn Marino     };
199*e4b17023SJohn Marino 
200*e4b17023SJohn Marino   template<typename _Tp>
201*e4b17023SJohn Marino     inline
202*e4b17023SJohn Marino     slice_array<_Tp>::slice_array(_Array<_Tp> __a, const slice& __s)
203*e4b17023SJohn Marino     : _M_sz(__s.size()), _M_stride(__s.stride()),
204*e4b17023SJohn Marino       _M_array(__a.begin() + __s.start()) {}
205*e4b17023SJohn Marino 
206*e4b17023SJohn Marino   template<typename _Tp>
207*e4b17023SJohn Marino     inline
208*e4b17023SJohn Marino     slice_array<_Tp>::slice_array(const slice_array<_Tp>& a)
209*e4b17023SJohn Marino     : _M_sz(a._M_sz), _M_stride(a._M_stride), _M_array(a._M_array) {}
210*e4b17023SJohn Marino 
211*e4b17023SJohn Marino   //    template<typename _Tp>
212*e4b17023SJohn Marino   //    inline slice_array<_Tp>::~slice_array () {}
213*e4b17023SJohn Marino 
214*e4b17023SJohn Marino   template<typename _Tp>
215*e4b17023SJohn Marino     inline slice_array<_Tp>&
216*e4b17023SJohn Marino     slice_array<_Tp>::operator=(const slice_array<_Tp>& __a)
217*e4b17023SJohn Marino     {
218*e4b17023SJohn Marino       std::__valarray_copy(__a._M_array, __a._M_sz, __a._M_stride,
219*e4b17023SJohn Marino 			   _M_array, _M_stride);
220*e4b17023SJohn Marino       return *this;
221*e4b17023SJohn Marino     }
222*e4b17023SJohn Marino 
223*e4b17023SJohn Marino   template<typename _Tp>
224*e4b17023SJohn Marino     inline void
225*e4b17023SJohn Marino     slice_array<_Tp>::operator=(const _Tp& __t) const
226*e4b17023SJohn Marino     { std::__valarray_fill(_M_array, _M_sz, _M_stride, __t); }
227*e4b17023SJohn Marino 
228*e4b17023SJohn Marino   template<typename _Tp>
229*e4b17023SJohn Marino     inline void
230*e4b17023SJohn Marino     slice_array<_Tp>::operator=(const valarray<_Tp>& __v) const
231*e4b17023SJohn Marino     { std::__valarray_copy(_Array<_Tp>(__v), _M_array, _M_sz, _M_stride); }
232*e4b17023SJohn Marino 
233*e4b17023SJohn Marino   template<typename _Tp>
234*e4b17023SJohn Marino   template<class _Dom>
235*e4b17023SJohn Marino     inline void
236*e4b17023SJohn Marino     slice_array<_Tp>::operator=(const _Expr<_Dom,_Tp>& __e) const
237*e4b17023SJohn Marino     { std::__valarray_copy(__e, _M_sz, _M_array, _M_stride); }
238*e4b17023SJohn Marino 
239*e4b17023SJohn Marino #undef _DEFINE_VALARRAY_OPERATOR
240*e4b17023SJohn Marino #define _DEFINE_VALARRAY_OPERATOR(_Op,_Name)				\
241*e4b17023SJohn Marino   template<typename _Tp>						\
242*e4b17023SJohn Marino     inline void								\
243*e4b17023SJohn Marino     slice_array<_Tp>::operator _Op##=(const valarray<_Tp>& __v) const	\
244*e4b17023SJohn Marino     {									\
245*e4b17023SJohn Marino       _Array_augmented_##_Name(_M_array, _M_sz, _M_stride, _Array<_Tp>(__v));\
246*e4b17023SJohn Marino     }									\
247*e4b17023SJohn Marino 									\
248*e4b17023SJohn Marino   template<typename _Tp>                                                \
249*e4b17023SJohn Marino     template<class _Dom>				                \
250*e4b17023SJohn Marino       inline void							\
251*e4b17023SJohn Marino       slice_array<_Tp>::operator _Op##=(const _Expr<_Dom,_Tp>& __e) const\
252*e4b17023SJohn Marino       {									\
253*e4b17023SJohn Marino 	  _Array_augmented_##_Name(_M_array, _M_stride, __e, _M_sz);	\
254*e4b17023SJohn Marino       }
255*e4b17023SJohn Marino 
256*e4b17023SJohn Marino 
257*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(*, __multiplies)
258*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(/, __divides)
259*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(%, __modulus)
260*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(+, __plus)
261*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(-, __minus)
262*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(^, __bitwise_xor)
263*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(&, __bitwise_and)
264*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(|, __bitwise_or)
265*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(<<, __shift_left)
266*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(>>, __shift_right)
267*e4b17023SJohn Marino 
268*e4b17023SJohn Marino #undef _DEFINE_VALARRAY_OPERATOR
269*e4b17023SJohn Marino 
270*e4b17023SJohn Marino   // @} group numeric_arrays
271*e4b17023SJohn Marino 
272*e4b17023SJohn Marino _GLIBCXX_END_NAMESPACE_VERSION
273*e4b17023SJohn Marino } // namespace
274*e4b17023SJohn Marino 
275*e4b17023SJohn Marino #endif /* _SLICE_ARRAY_H */
276