100db7afdSDavid E. O'Brien // The template and inlines for the -*- C++ -*- gslice_array 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_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 _GSLICE_ARRAY_H
39ffeaf689SAlexander Kabaev #define _GSLICE_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)
4400db7afdSDavid E. O'Brien 
45ffeaf689SAlexander Kabaev   /**
46ffeaf689SAlexander Kabaev    *  @brief  Reference to multi-dimensional subset of an array.
47ffeaf689SAlexander Kabaev    *
48ffeaf689SAlexander Kabaev    *  A gslice_array is a reference to the actual elements of an array
49ffeaf689SAlexander Kabaev    *  specified by a gslice.  The way to get a gslice_array is to call
50ffeaf689SAlexander Kabaev    *  operator[](gslice) on a valarray.  The returned gslice_array then
51ffeaf689SAlexander Kabaev    *  permits carrying operations out on the referenced subset of elements in
52ffeaf689SAlexander Kabaev    *  the original valarray.  For example, operator+=(valarray) will add
53ffeaf689SAlexander Kabaev    *  values to the subset of elements in the underlying valarray this
54ffeaf689SAlexander Kabaev    *  gslice_array refers to.
55ffeaf689SAlexander Kabaev    *
56ffeaf689SAlexander Kabaev    *  @param  Tp  Element type.
57ffeaf689SAlexander Kabaev    */
581b86b14eSAlexander Kabaev   template<typename _Tp>
591b86b14eSAlexander Kabaev     class gslice_array
6000db7afdSDavid E. O'Brien     {
6100db7afdSDavid E. O'Brien     public:
6200db7afdSDavid E. O'Brien       typedef _Tp value_type;
6300db7afdSDavid E. O'Brien 
64ffeaf689SAlexander Kabaev       // _GLIBCXX_RESOLVE_LIB_DEFECTS
65ffeaf689SAlexander Kabaev       // 253. valarray helper functions are almost entirely useless
66ffeaf689SAlexander Kabaev 
67ffeaf689SAlexander Kabaev       ///  Copy constructor.  Both slices refer to the same underlying array.
68ffeaf689SAlexander Kabaev       gslice_array(const gslice_array&);
69ffeaf689SAlexander Kabaev 
70ffeaf689SAlexander Kabaev       ///  Assignment operator.  Assigns slice elements to corresponding
71ffeaf689SAlexander Kabaev       ///  elements of @a a.
72ffeaf689SAlexander Kabaev       gslice_array& operator=(const gslice_array&);
73ffeaf689SAlexander Kabaev 
74ffeaf689SAlexander Kabaev       ///  Assign slice elements to corresponding elements of @a v.
7500db7afdSDavid E. O'Brien       void operator=(const valarray<_Tp>&) const;
76ffeaf689SAlexander Kabaev       ///  Multiply slice elements by corresponding elements of @a v.
7700db7afdSDavid E. O'Brien       void operator*=(const valarray<_Tp>&) const;
78ffeaf689SAlexander Kabaev       ///  Divide slice elements by corresponding elements of @a v.
7900db7afdSDavid E. O'Brien       void operator/=(const valarray<_Tp>&) const;
80ffeaf689SAlexander Kabaev       ///  Modulo slice elements by corresponding elements of @a v.
8100db7afdSDavid E. O'Brien       void operator%=(const valarray<_Tp>&) const;
82ffeaf689SAlexander Kabaev       ///  Add corresponding elements of @a v to slice elements.
8300db7afdSDavid E. O'Brien       void operator+=(const valarray<_Tp>&) const;
84ffeaf689SAlexander Kabaev       ///  Subtract corresponding elements of @a v from slice elements.
8500db7afdSDavid E. O'Brien       void operator-=(const valarray<_Tp>&) const;
86ffeaf689SAlexander Kabaev       ///  Logical xor slice elements with corresponding elements of @a v.
8700db7afdSDavid E. O'Brien       void operator^=(const valarray<_Tp>&) const;
88ffeaf689SAlexander Kabaev       ///  Logical and slice elements with corresponding elements of @a v.
8900db7afdSDavid E. O'Brien       void operator&=(const valarray<_Tp>&) const;
90ffeaf689SAlexander Kabaev       ///  Logical or slice elements with corresponding elements of @a v.
9100db7afdSDavid E. O'Brien       void operator|=(const valarray<_Tp>&) const;
92ffeaf689SAlexander Kabaev       ///  Left shift slice elements by corresponding elements of @a v.
9300db7afdSDavid E. O'Brien       void operator<<=(const valarray<_Tp>&) const;
94ffeaf689SAlexander Kabaev       ///  Right shift slice elements by corresponding elements of @a v.
9500db7afdSDavid E. O'Brien       void operator>>=(const valarray<_Tp>&) const;
96ffeaf689SAlexander Kabaev       ///  Assign all slice elements to @a t.
971b86b14eSAlexander Kabaev       void operator=(const _Tp&) const;
9800db7afdSDavid E. O'Brien 
9900db7afdSDavid E. O'Brien       template<class _Dom>
10000db7afdSDavid E. O'Brien         void operator=(const _Expr<_Dom, _Tp>&) const;
10100db7afdSDavid E. O'Brien       template<class _Dom>
10200db7afdSDavid E. O'Brien         void operator*=(const _Expr<_Dom, _Tp>&) const;
10300db7afdSDavid E. O'Brien       template<class _Dom>
10400db7afdSDavid E. O'Brien         void operator/=(const _Expr<_Dom, _Tp>&) const;
10500db7afdSDavid E. O'Brien       template<class _Dom>
10600db7afdSDavid E. O'Brien         void operator%=(const _Expr<_Dom, _Tp>&) const;
10700db7afdSDavid E. O'Brien       template<class _Dom>
10800db7afdSDavid E. O'Brien         void operator+=(const _Expr<_Dom, _Tp>&) const;
10900db7afdSDavid E. O'Brien       template<class _Dom>
11000db7afdSDavid E. O'Brien         void operator-=(const _Expr<_Dom, _Tp>&) const;
11100db7afdSDavid E. O'Brien       template<class _Dom>
11200db7afdSDavid E. O'Brien         void operator^=(const _Expr<_Dom, _Tp>&) const;
11300db7afdSDavid E. O'Brien       template<class _Dom>
11400db7afdSDavid E. O'Brien         void operator&=(const _Expr<_Dom, _Tp>&) const;
11500db7afdSDavid E. O'Brien       template<class _Dom>
11600db7afdSDavid E. O'Brien         void operator|=(const _Expr<_Dom, _Tp>&) const;
11700db7afdSDavid E. O'Brien       template<class _Dom>
11800db7afdSDavid E. O'Brien         void operator<<=(const _Expr<_Dom, _Tp>&) const;
11900db7afdSDavid E. O'Brien       template<class _Dom>
12000db7afdSDavid E. O'Brien         void operator>>=(const _Expr<_Dom, _Tp>&) const;
12100db7afdSDavid E. O'Brien 
12200db7afdSDavid E. O'Brien     private:
12300db7afdSDavid E. O'Brien       _Array<_Tp>    _M_array;
12400db7afdSDavid E. O'Brien       const valarray<size_t>& _M_index;
12500db7afdSDavid E. O'Brien 
12600db7afdSDavid E. O'Brien       friend class valarray<_Tp>;
12700db7afdSDavid E. O'Brien 
12800db7afdSDavid E. O'Brien       gslice_array(_Array<_Tp>, const valarray<size_t>&);
12900db7afdSDavid E. O'Brien 
13000db7afdSDavid E. O'Brien       // not implemented
13100db7afdSDavid E. O'Brien       gslice_array();
13200db7afdSDavid E. O'Brien     };
13300db7afdSDavid E. O'Brien 
13400db7afdSDavid E. O'Brien   template<typename _Tp>
13500db7afdSDavid E. O'Brien     inline
gslice_array(_Array<_Tp> __a,const valarray<size_t> & __i)13600db7afdSDavid E. O'Brien     gslice_array<_Tp>::gslice_array(_Array<_Tp> __a,
13700db7afdSDavid E. O'Brien 				    const valarray<size_t>& __i)
13800db7afdSDavid E. O'Brien     : _M_array(__a), _M_index(__i) {}
13900db7afdSDavid E. O'Brien 
14000db7afdSDavid E. O'Brien   template<typename _Tp>
14100db7afdSDavid E. O'Brien     inline
gslice_array(const gslice_array<_Tp> & __a)14200db7afdSDavid E. O'Brien     gslice_array<_Tp>::gslice_array(const gslice_array<_Tp>& __a)
14300db7afdSDavid E. O'Brien     : _M_array(__a._M_array), _M_index(__a._M_index) {}
14400db7afdSDavid E. O'Brien 
14500db7afdSDavid E. O'Brien   template<typename _Tp>
146ffeaf689SAlexander Kabaev     inline gslice_array<_Tp>&
147ffeaf689SAlexander Kabaev     gslice_array<_Tp>::operator=(const gslice_array<_Tp>& __a)
148ffeaf689SAlexander Kabaev     {
149ffeaf689SAlexander Kabaev       std::__valarray_copy(_Array<_Tp>(__a._M_array),
150ffeaf689SAlexander Kabaev 			   _Array<size_t>(__a._M_index), _M_index.size(),
151ffeaf689SAlexander Kabaev 			   _M_array, _Array<size_t>(_M_index));
152ffeaf689SAlexander Kabaev       return *this;
153ffeaf689SAlexander Kabaev     }
154ffeaf689SAlexander Kabaev 
155ffeaf689SAlexander Kabaev   template<typename _Tp>
15600db7afdSDavid E. O'Brien     inline void
1571b86b14eSAlexander Kabaev     gslice_array<_Tp>::operator=(const _Tp& __t) const
15800db7afdSDavid E. O'Brien     {
159ffeaf689SAlexander Kabaev       std::__valarray_fill(_M_array, _Array<size_t>(_M_index),
16000db7afdSDavid E. O'Brien 			   _M_index.size(), __t);
16100db7afdSDavid E. O'Brien     }
16200db7afdSDavid E. O'Brien 
16300db7afdSDavid E. O'Brien   template<typename _Tp>
16400db7afdSDavid E. O'Brien     inline void
16500db7afdSDavid E. O'Brien     gslice_array<_Tp>::operator=(const valarray<_Tp>& __v) const
16600db7afdSDavid E. O'Brien     {
167ffeaf689SAlexander Kabaev       std::__valarray_copy(_Array<_Tp>(__v), __v.size(),
16800db7afdSDavid E. O'Brien 			   _M_array, _Array<size_t>(_M_index));
16900db7afdSDavid E. O'Brien     }
17000db7afdSDavid E. O'Brien 
17100db7afdSDavid E. O'Brien   template<typename _Tp>
1721b86b14eSAlexander Kabaev     template<class _Dom>
17300db7afdSDavid E. O'Brien       inline void
1741b86b14eSAlexander Kabaev       gslice_array<_Tp>::operator=(const _Expr<_Dom, _Tp>& __e) const
17500db7afdSDavid E. O'Brien       {
176ffeaf689SAlexander Kabaev 	std::__valarray_copy (__e, _M_index.size(), _M_array,
17700db7afdSDavid E. O'Brien 			      _Array<size_t>(_M_index));
17800db7afdSDavid E. O'Brien       }
17900db7afdSDavid E. O'Brien 
18000db7afdSDavid E. O'Brien #undef _DEFINE_VALARRAY_OPERATOR
1811b86b14eSAlexander Kabaev #define _DEFINE_VALARRAY_OPERATOR(_Op, _Name)				\
18200db7afdSDavid E. O'Brien   template<typename _Tp>						\
18300db7afdSDavid E. O'Brien     inline void								\
1841b86b14eSAlexander Kabaev     gslice_array<_Tp>::operator _Op##=(const valarray<_Tp>& __v) const	\
18500db7afdSDavid E. O'Brien     {									\
1861b86b14eSAlexander Kabaev       _Array_augmented_##_Name(_M_array, _Array<size_t>(_M_index),	\
18700db7afdSDavid E. O'Brien 			       _Array<_Tp>(__v), __v.size());		\
18800db7afdSDavid E. O'Brien     }									\
18900db7afdSDavid E. O'Brien 									\
1901b86b14eSAlexander Kabaev   template<typename _Tp>                                                \
1911b86b14eSAlexander Kabaev     template<class _Dom>				                \
19200db7afdSDavid E. O'Brien       inline void							\
1931b86b14eSAlexander Kabaev       gslice_array<_Tp>::operator _Op##= (const _Expr<_Dom, _Tp>& __e) const\
19400db7afdSDavid E. O'Brien       {									\
1951b86b14eSAlexander Kabaev 	_Array_augmented_##_Name(_M_array, _Array<size_t>(_M_index), __e,\
19600db7afdSDavid E. O'Brien 				 _M_index.size());			\
19700db7afdSDavid E. O'Brien       }
19800db7afdSDavid E. O'Brien 
1991b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(*, __multiplies)
2001b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(/, __divides)
2011b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(%, __modulus)
2021b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(+, __plus)
2031b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(-, __minus)
2041b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(^, __bitwise_xor)
2051b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(&, __bitwise_and)
2061b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(|, __bitwise_or)
2071b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(<<, __shift_left)
2081b86b14eSAlexander Kabaev _DEFINE_VALARRAY_OPERATOR(>>, __shift_right)
20900db7afdSDavid E. O'Brien 
21000db7afdSDavid E. O'Brien #undef _DEFINE_VALARRAY_OPERATOR
21100db7afdSDavid E. O'Brien 
212*f8a1b7d9SAlexander Kabaev _GLIBCXX_END_NAMESPACE
21300db7afdSDavid E. O'Brien 
214ffeaf689SAlexander Kabaev #endif /* _GSLICE_ARRAY_H */
215