100db7afdSDavid E. O'Brien // The template and inlines for the -*- C++ -*- slice_array class.
200db7afdSDavid E. O'Brien 
3*f8a1b7d9SAlexander Kabaev // Copyright (C) 1997, 1998, 1999, 2000, 2001, 2002, 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 slice_array.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 _SLICE_ARRAY_H
39ffeaf689SAlexander Kabaev #define _SLICE_ARRAY_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)
44*f8a1b7d9SAlexander Kabaev 
45ffeaf689SAlexander Kabaev   /**
46ffeaf689SAlexander Kabaev    *  @brief  Class defining one-dimensional subset of an array.
47ffeaf689SAlexander Kabaev    *
48ffeaf689SAlexander Kabaev    *  The slice class represents a one-dimensional subset of an array,
49ffeaf689SAlexander Kabaev    *  specified by three parameters: start offset, size, and stride.  The
50ffeaf689SAlexander Kabaev    *  start offset is the index of the first element of the array that is part
51ffeaf689SAlexander Kabaev    *  of the subset.  The size is the total number of elements in the subset.
52ffeaf689SAlexander Kabaev    *  Stride is the distance between each successive array element to include
53ffeaf689SAlexander Kabaev    *  in the subset.
54ffeaf689SAlexander Kabaev    *
55ffeaf689SAlexander Kabaev    *  For example, with an array of size 10, and a slice with offset 1, size 3
56ffeaf689SAlexander Kabaev    *  and stride 2, the subset consists of array elements 1, 3, and 5.
57ffeaf689SAlexander Kabaev    */
581b86b14eSAlexander Kabaev   class slice
591b86b14eSAlexander Kabaev   {
601b86b14eSAlexander Kabaev   public:
61ffeaf689SAlexander Kabaev     ///  Construct an empty slice.
621b86b14eSAlexander Kabaev     slice();
63ffeaf689SAlexander Kabaev 
64ffeaf689SAlexander Kabaev     /**
65ffeaf689SAlexander Kabaev      *  @brief  Construct a slice.
66ffeaf689SAlexander Kabaev      *
67ffeaf689SAlexander Kabaev      *  @param  o  Offset in array of first element.
68ffeaf689SAlexander Kabaev      *  @param  d  Number of elements in slice.
69ffeaf689SAlexander Kabaev      *  @param  s  Stride between array elements.
70ffeaf689SAlexander Kabaev      */
711b86b14eSAlexander Kabaev     slice(size_t, size_t, size_t);
721b86b14eSAlexander Kabaev 
73ffeaf689SAlexander Kabaev     ///  Return array offset of first slice element.
741b86b14eSAlexander Kabaev     size_t start() const;
75ffeaf689SAlexander Kabaev     ///  Return size of slice.
761b86b14eSAlexander Kabaev     size_t size() const;
77ffeaf689SAlexander Kabaev     ///  Return array stride of slice.
781b86b14eSAlexander Kabaev     size_t stride() const;
791b86b14eSAlexander Kabaev 
801b86b14eSAlexander Kabaev   private:
811b86b14eSAlexander Kabaev     size_t _M_off;                      // offset
821b86b14eSAlexander Kabaev     size_t _M_sz;			// size
831b86b14eSAlexander Kabaev     size_t _M_st;			// stride unit
841b86b14eSAlexander Kabaev   };
851b86b14eSAlexander Kabaev 
861b86b14eSAlexander Kabaev   // The default constructor constructor is not required to initialize
871b86b14eSAlexander Kabaev   // data members with any meaningful values, so we choose to do nothing.
881b86b14eSAlexander Kabaev   inline
slice()891b86b14eSAlexander Kabaev   slice::slice() {}
901b86b14eSAlexander Kabaev 
911b86b14eSAlexander Kabaev   inline
slice(size_t __o,size_t __d,size_t __s)921b86b14eSAlexander Kabaev   slice::slice(size_t __o, size_t __d, size_t __s)
931b86b14eSAlexander Kabaev   : _M_off(__o), _M_sz(__d), _M_st(__s) {}
941b86b14eSAlexander Kabaev 
951b86b14eSAlexander Kabaev   inline size_t
start()961b86b14eSAlexander Kabaev   slice::start() const
971b86b14eSAlexander Kabaev   { return _M_off; }
981b86b14eSAlexander Kabaev 
991b86b14eSAlexander Kabaev   inline size_t
size()1001b86b14eSAlexander Kabaev   slice::size() const
1011b86b14eSAlexander Kabaev   { return _M_sz; }
1021b86b14eSAlexander Kabaev 
1031b86b14eSAlexander Kabaev   inline size_t
stride()1041b86b14eSAlexander Kabaev   slice::stride() const
1051b86b14eSAlexander Kabaev   { return _M_st; }
10600db7afdSDavid E. O'Brien 
107ffeaf689SAlexander Kabaev   /**
108ffeaf689SAlexander Kabaev    *  @brief  Reference to one-dimensional subset of an array.
109ffeaf689SAlexander Kabaev    *
110ffeaf689SAlexander Kabaev    *  A slice_array is a reference to the actual elements of an array
111ffeaf689SAlexander Kabaev    *  specified by a slice.  The way to get a slice_array is to call
112ffeaf689SAlexander Kabaev    *  operator[](slice) on a valarray.  The returned slice_array then permits
113ffeaf689SAlexander Kabaev    *  carrying operations out on the referenced subset of elements in the
114ffeaf689SAlexander Kabaev    *  original valarray.  For example, operator+=(valarray) will add values
115ffeaf689SAlexander Kabaev    *  to the subset of elements in the underlying valarray this slice_array
116ffeaf689SAlexander Kabaev    *  refers to.
117ffeaf689SAlexander Kabaev    *
118ffeaf689SAlexander Kabaev    *  @param  Tp  Element type.
119ffeaf689SAlexander Kabaev    */
12000db7afdSDavid E. O'Brien   template<typename _Tp>
12100db7afdSDavid E. O'Brien     class slice_array
12200db7afdSDavid E. O'Brien     {
12300db7afdSDavid E. O'Brien     public:
12400db7afdSDavid E. O'Brien       typedef _Tp value_type;
12500db7afdSDavid E. O'Brien 
126ffeaf689SAlexander Kabaev       // _GLIBCXX_RESOLVE_LIB_DEFECTS
127ffeaf689SAlexander Kabaev       // 253. valarray helper functions are almost entirely useless
128ffeaf689SAlexander Kabaev 
129ffeaf689SAlexander Kabaev       ///  Copy constructor.  Both slices refer to the same underlying array.
13000db7afdSDavid E. O'Brien       slice_array(const slice_array&);
13100db7afdSDavid E. O'Brien 
132ffeaf689SAlexander Kabaev       ///  Assignment operator.  Assigns slice elements to corresponding
133ffeaf689SAlexander Kabaev       ///  elements of @a a.
13400db7afdSDavid E. O'Brien       slice_array& operator=(const slice_array&);
13500db7afdSDavid E. O'Brien 
136ffeaf689SAlexander Kabaev       ///  Assign slice elements to corresponding elements of @a v.
13700db7afdSDavid E. O'Brien       void operator=(const valarray<_Tp>&) const;
138ffeaf689SAlexander Kabaev       ///  Multiply slice elements by corresponding elements of @a v.
13900db7afdSDavid E. O'Brien       void operator*=(const valarray<_Tp>&) const;
140ffeaf689SAlexander Kabaev       ///  Divide slice elements by corresponding elements of @a v.
14100db7afdSDavid E. O'Brien       void operator/=(const valarray<_Tp>&) const;
142ffeaf689SAlexander Kabaev       ///  Modulo slice elements by corresponding elements of @a v.
14300db7afdSDavid E. O'Brien       void operator%=(const valarray<_Tp>&) const;
144ffeaf689SAlexander Kabaev       ///  Add corresponding elements of @a v to slice elements.
14500db7afdSDavid E. O'Brien       void operator+=(const valarray<_Tp>&) const;
146ffeaf689SAlexander Kabaev       ///  Subtract corresponding elements of @a v from slice elements.
14700db7afdSDavid E. O'Brien       void operator-=(const valarray<_Tp>&) const;
148ffeaf689SAlexander Kabaev       ///  Logical xor slice elements with corresponding elements of @a v.
14900db7afdSDavid E. O'Brien       void operator^=(const valarray<_Tp>&) const;
150ffeaf689SAlexander Kabaev       ///  Logical and slice elements with corresponding elements of @a v.
15100db7afdSDavid E. O'Brien       void operator&=(const valarray<_Tp>&) const;
152ffeaf689SAlexander Kabaev       ///  Logical or slice elements with corresponding elements of @a v.
15300db7afdSDavid E. O'Brien       void operator|=(const valarray<_Tp>&) const;
154ffeaf689SAlexander Kabaev       ///  Left shift slice elements by corresponding elements of @a v.
15500db7afdSDavid E. O'Brien       void operator<<=(const valarray<_Tp>&) const;
156ffeaf689SAlexander Kabaev       ///  Right shift slice elements by corresponding elements of @a v.
15700db7afdSDavid E. O'Brien       void operator>>=(const valarray<_Tp>&) const;
158ffeaf689SAlexander Kabaev       ///  Assign all slice elements to @a t.
1591b86b14eSAlexander Kabaev       void operator=(const _Tp &) const;
16000db7afdSDavid E. O'Brien       //        ~slice_array ();
16100db7afdSDavid E. O'Brien 
16200db7afdSDavid E. O'Brien       template<class _Dom>
16300db7afdSDavid E. O'Brien         void operator=(const _Expr<_Dom, _Tp>&) const;
16400db7afdSDavid E. O'Brien       template<class _Dom>
16500db7afdSDavid E. O'Brien 	void operator*=(const _Expr<_Dom, _Tp>&) const;
16600db7afdSDavid E. O'Brien       template<class _Dom>
16700db7afdSDavid E. O'Brien 	void operator/=(const _Expr<_Dom, _Tp>&) const;
16800db7afdSDavid E. O'Brien       template<class _Dom>
16900db7afdSDavid E. O'Brien 	void operator%=(const _Expr<_Dom, _Tp>&) const;
17000db7afdSDavid E. O'Brien       template<class _Dom>
17100db7afdSDavid E. O'Brien 	void operator+=(const _Expr<_Dom, _Tp>&) const;
17200db7afdSDavid E. O'Brien       template<class _Dom>
17300db7afdSDavid E. O'Brien 	void operator-=(const _Expr<_Dom, _Tp>&) const;
17400db7afdSDavid E. O'Brien       template<class _Dom>
17500db7afdSDavid E. O'Brien 	void operator^=(const _Expr<_Dom, _Tp>&) const;
17600db7afdSDavid E. O'Brien       template<class _Dom>
17700db7afdSDavid E. O'Brien 	void operator&=(const _Expr<_Dom, _Tp>&) const;
17800db7afdSDavid E. O'Brien       template<class _Dom>
17900db7afdSDavid E. O'Brien 	void operator|=(const _Expr<_Dom, _Tp>&) const;
18000db7afdSDavid E. O'Brien       template<class _Dom>
18100db7afdSDavid E. O'Brien 	void operator<<=(const _Expr<_Dom, _Tp>&) const;
18200db7afdSDavid E. O'Brien       template<class _Dom>
18300db7afdSDavid E. O'Brien 	void operator>>=(const _Expr<_Dom, _Tp>&) const;
18400db7afdSDavid E. O'Brien 
18500db7afdSDavid E. O'Brien     private:
18600db7afdSDavid E. O'Brien       friend class valarray<_Tp>;
18700db7afdSDavid E. O'Brien       slice_array(_Array<_Tp>, const slice&);
18800db7afdSDavid E. O'Brien 
18900db7afdSDavid E. O'Brien       const size_t      _M_sz;
19000db7afdSDavid E. O'Brien       const size_t      _M_stride;
19100db7afdSDavid E. O'Brien       const _Array<_Tp> _M_array;
19200db7afdSDavid E. O'Brien 
19300db7afdSDavid E. O'Brien       // not implemented
19400db7afdSDavid E. O'Brien       slice_array();
19500db7afdSDavid E. O'Brien     };
19600db7afdSDavid E. O'Brien 
19700db7afdSDavid E. O'Brien   template<typename _Tp>
1981b86b14eSAlexander Kabaev     inline
slice_array(_Array<_Tp> __a,const slice & __s)1991b86b14eSAlexander Kabaev     slice_array<_Tp>::slice_array(_Array<_Tp> __a, const slice& __s)
20000db7afdSDavid E. O'Brien     : _M_sz(__s.size()), _M_stride(__s.stride()),
20100db7afdSDavid E. O'Brien       _M_array(__a.begin() + __s.start()) {}
20200db7afdSDavid E. O'Brien 
20300db7afdSDavid E. O'Brien   template<typename _Tp>
2041b86b14eSAlexander Kabaev     inline
slice_array(const slice_array<_Tp> & a)2051b86b14eSAlexander Kabaev     slice_array<_Tp>::slice_array(const slice_array<_Tp>& a)
20600db7afdSDavid E. O'Brien     : _M_sz(a._M_sz), _M_stride(a._M_stride), _M_array(a._M_array) {}
20700db7afdSDavid E. O'Brien 
20800db7afdSDavid E. O'Brien   //    template<typename _Tp>
20900db7afdSDavid E. O'Brien   //    inline slice_array<_Tp>::~slice_array () {}
21000db7afdSDavid E. O'Brien 
21100db7afdSDavid E. O'Brien   template<typename _Tp>
21200db7afdSDavid E. O'Brien     inline slice_array<_Tp>&
21300db7afdSDavid E. O'Brien     slice_array<_Tp>::operator=(const slice_array<_Tp>& __a)
21400db7afdSDavid E. O'Brien     {
215ffeaf689SAlexander Kabaev       std::__valarray_copy(__a._M_array, __a._M_sz, __a._M_stride,
21661e268fbSAlexander Kabaev 			   _M_array, _M_stride);
21700db7afdSDavid E. O'Brien       return *this;
21800db7afdSDavid E. O'Brien     }
21900db7afdSDavid E. O'Brien 
22000db7afdSDavid E. O'Brien   template<typename _Tp>
22100db7afdSDavid E. O'Brien     inline void
2221b86b14eSAlexander Kabaev     slice_array<_Tp>::operator=(const _Tp& __t) const
223ffeaf689SAlexander Kabaev     { std::__valarray_fill(_M_array, _M_sz, _M_stride, __t); }
22400db7afdSDavid E. O'Brien 
22500db7afdSDavid E. O'Brien   template<typename _Tp>
22600db7afdSDavid E. O'Brien     inline void
22700db7afdSDavid E. O'Brien     slice_array<_Tp>::operator=(const valarray<_Tp>& __v) const
228ffeaf689SAlexander Kabaev     { std::__valarray_copy(_Array<_Tp>(__v), _M_array, _M_sz, _M_stride); }
22900db7afdSDavid E. O'Brien 
23000db7afdSDavid E. O'Brien   template<typename _Tp>
23100db7afdSDavid E. O'Brien   template<class _Dom>
23200db7afdSDavid E. O'Brien     inline void
23300db7afdSDavid E. O'Brien     slice_array<_Tp>::operator=(const _Expr<_Dom,_Tp>& __e) const
234ffeaf689SAlexander Kabaev     { std::__valarray_copy(__e, _M_sz, _M_array, _M_stride); }
23500db7afdSDavid E. O'Brien 
23600db7afdSDavid E. O'Brien #undef _DEFINE_VALARRAY_OPERATOR
2371b86b14eSAlexander Kabaev #define _DEFINE_VALARRAY_OPERATOR(_Op,_Name)				\
23800db7afdSDavid E. O'Brien   template<typename _Tp>						\
23900db7afdSDavid E. O'Brien     inline void								\
2401b86b14eSAlexander Kabaev     slice_array<_Tp>::operator _Op##=(const valarray<_Tp>& __v) const	\
24100db7afdSDavid E. O'Brien     {									\
2421b86b14eSAlexander Kabaev       _Array_augmented_##_Name(_M_array, _M_sz, _M_stride, _Array<_Tp>(__v));\
24300db7afdSDavid E. O'Brien     }									\
24400db7afdSDavid E. O'Brien 									\
2451b86b14eSAlexander Kabaev   template<typename _Tp>                                                \
2461b86b14eSAlexander Kabaev     template<class _Dom>				                \
24700db7afdSDavid E. O'Brien       inline void							\
2481b86b14eSAlexander Kabaev       slice_array<_Tp>::operator _Op##=(const _Expr<_Dom,_Tp>& __e) const\
24900db7afdSDavid E. O'Brien       {									\
2501b86b14eSAlexander Kabaev 	  _Array_augmented_##_Name(_M_array, _M_stride, __e, _M_sz);	\
25100db7afdSDavid E. O'Brien       }
25200db7afdSDavid E. O'Brien 
25300db7afdSDavid E. O'Brien 
2541b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(*, __multiplies)
2551b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(/, __divides)
2561b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(%, __modulus)
2571b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(+, __plus)
2581b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(-, __minus)
2591b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(^, __bitwise_xor)
2601b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(&, __bitwise_and)
2611b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(|, __bitwise_or)
2621b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(<<, __shift_left)
2631b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(>>, __shift_right)
26400db7afdSDavid E. O'Brien 
26500db7afdSDavid E. O'Brien #undef _DEFINE_VALARRAY_OPERATOR
26600db7afdSDavid E. O'Brien 
267*f8a1b7d9SAlexander Kabaev _GLIBCXX_END_NAMESPACE
26800db7afdSDavid E. O'Brien 
269ffeaf689SAlexander Kabaev #endif /* _SLICE_ARRAY_H */
270