1*e4b17023SJohn Marino // The template and inlines for the -*- C++ -*- gslice_array class.
2*e4b17023SJohn Marino
3*e4b17023SJohn Marino // Copyright (C) 1997, 1998, 1999, 2000, 2001, 2004, 2005, 2009, 2010
4*e4b17023SJohn Marino // 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_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 _GSLICE_ARRAY_H
34*e4b17023SJohn Marino #define _GSLICE_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 Reference to multi-dimensional subset of an array.
49*e4b17023SJohn Marino *
50*e4b17023SJohn Marino * A gslice_array is a reference to the actual elements of an array
51*e4b17023SJohn Marino * specified by a gslice. The way to get a gslice_array is to call
52*e4b17023SJohn Marino * operator[](gslice) on a valarray. The returned gslice_array then
53*e4b17023SJohn Marino * permits carrying operations out on the referenced subset of elements in
54*e4b17023SJohn Marino * the original valarray. For example, operator+=(valarray) will add
55*e4b17023SJohn Marino * values to the subset of elements in the underlying valarray this
56*e4b17023SJohn Marino * gslice_array refers to.
57*e4b17023SJohn Marino *
58*e4b17023SJohn Marino * @param Tp Element type.
59*e4b17023SJohn Marino */
60*e4b17023SJohn Marino template<typename _Tp>
61*e4b17023SJohn Marino class gslice_array
62*e4b17023SJohn Marino {
63*e4b17023SJohn Marino public:
64*e4b17023SJohn Marino typedef _Tp value_type;
65*e4b17023SJohn Marino
66*e4b17023SJohn Marino // _GLIBCXX_RESOLVE_LIB_DEFECTS
67*e4b17023SJohn Marino // 253. valarray helper functions are almost entirely useless
68*e4b17023SJohn Marino
69*e4b17023SJohn Marino /// Copy constructor. Both slices refer to the same underlying array.
70*e4b17023SJohn Marino gslice_array(const gslice_array&);
71*e4b17023SJohn Marino
72*e4b17023SJohn Marino /// Assignment operator. Assigns slice elements to corresponding
73*e4b17023SJohn Marino /// elements of @a a.
74*e4b17023SJohn Marino gslice_array& operator=(const gslice_array&);
75*e4b17023SJohn Marino
76*e4b17023SJohn Marino /// Assign slice elements to corresponding elements of @a v.
77*e4b17023SJohn Marino void operator=(const valarray<_Tp>&) const;
78*e4b17023SJohn Marino /// Multiply slice elements by corresponding elements of @a v.
79*e4b17023SJohn Marino void operator*=(const valarray<_Tp>&) const;
80*e4b17023SJohn Marino /// Divide slice elements by corresponding elements of @a v.
81*e4b17023SJohn Marino void operator/=(const valarray<_Tp>&) const;
82*e4b17023SJohn Marino /// Modulo slice elements by corresponding elements of @a v.
83*e4b17023SJohn Marino void operator%=(const valarray<_Tp>&) const;
84*e4b17023SJohn Marino /// Add corresponding elements of @a v to slice elements.
85*e4b17023SJohn Marino void operator+=(const valarray<_Tp>&) const;
86*e4b17023SJohn Marino /// Subtract corresponding elements of @a v from slice elements.
87*e4b17023SJohn Marino void operator-=(const valarray<_Tp>&) const;
88*e4b17023SJohn Marino /// Logical xor slice elements with corresponding elements of @a v.
89*e4b17023SJohn Marino void operator^=(const valarray<_Tp>&) const;
90*e4b17023SJohn Marino /// Logical and slice elements with corresponding elements of @a v.
91*e4b17023SJohn Marino void operator&=(const valarray<_Tp>&) const;
92*e4b17023SJohn Marino /// Logical or slice elements with corresponding elements of @a v.
93*e4b17023SJohn Marino void operator|=(const valarray<_Tp>&) const;
94*e4b17023SJohn Marino /// Left shift slice elements by corresponding elements of @a v.
95*e4b17023SJohn Marino void operator<<=(const valarray<_Tp>&) const;
96*e4b17023SJohn Marino /// Right shift slice elements by corresponding elements of @a v.
97*e4b17023SJohn Marino void operator>>=(const valarray<_Tp>&) const;
98*e4b17023SJohn Marino /// Assign all slice elements to @a t.
99*e4b17023SJohn Marino void operator=(const _Tp&) const;
100*e4b17023SJohn Marino
101*e4b17023SJohn Marino template<class _Dom>
102*e4b17023SJohn Marino void operator=(const _Expr<_Dom, _Tp>&) const;
103*e4b17023SJohn Marino template<class _Dom>
104*e4b17023SJohn Marino void operator*=(const _Expr<_Dom, _Tp>&) const;
105*e4b17023SJohn Marino template<class _Dom>
106*e4b17023SJohn Marino void operator/=(const _Expr<_Dom, _Tp>&) const;
107*e4b17023SJohn Marino template<class _Dom>
108*e4b17023SJohn Marino void operator%=(const _Expr<_Dom, _Tp>&) const;
109*e4b17023SJohn Marino template<class _Dom>
110*e4b17023SJohn Marino void operator+=(const _Expr<_Dom, _Tp>&) const;
111*e4b17023SJohn Marino template<class _Dom>
112*e4b17023SJohn Marino void operator-=(const _Expr<_Dom, _Tp>&) const;
113*e4b17023SJohn Marino template<class _Dom>
114*e4b17023SJohn Marino void operator^=(const _Expr<_Dom, _Tp>&) const;
115*e4b17023SJohn Marino template<class _Dom>
116*e4b17023SJohn Marino void operator&=(const _Expr<_Dom, _Tp>&) const;
117*e4b17023SJohn Marino template<class _Dom>
118*e4b17023SJohn Marino void operator|=(const _Expr<_Dom, _Tp>&) const;
119*e4b17023SJohn Marino template<class _Dom>
120*e4b17023SJohn Marino void operator<<=(const _Expr<_Dom, _Tp>&) const;
121*e4b17023SJohn Marino template<class _Dom>
122*e4b17023SJohn Marino void operator>>=(const _Expr<_Dom, _Tp>&) const;
123*e4b17023SJohn Marino
124*e4b17023SJohn Marino private:
125*e4b17023SJohn Marino _Array<_Tp> _M_array;
126*e4b17023SJohn Marino const valarray<size_t>& _M_index;
127*e4b17023SJohn Marino
128*e4b17023SJohn Marino friend class valarray<_Tp>;
129*e4b17023SJohn Marino
130*e4b17023SJohn Marino gslice_array(_Array<_Tp>, const valarray<size_t>&);
131*e4b17023SJohn Marino
132*e4b17023SJohn Marino // not implemented
133*e4b17023SJohn Marino gslice_array();
134*e4b17023SJohn Marino };
135*e4b17023SJohn Marino
136*e4b17023SJohn Marino template<typename _Tp>
137*e4b17023SJohn Marino inline
138*e4b17023SJohn Marino gslice_array<_Tp>::gslice_array(_Array<_Tp> __a,
139*e4b17023SJohn Marino const valarray<size_t>& __i)
140*e4b17023SJohn Marino : _M_array(__a), _M_index(__i) {}
141*e4b17023SJohn Marino
142*e4b17023SJohn Marino template<typename _Tp>
143*e4b17023SJohn Marino inline
144*e4b17023SJohn Marino gslice_array<_Tp>::gslice_array(const gslice_array<_Tp>& __a)
145*e4b17023SJohn Marino : _M_array(__a._M_array), _M_index(__a._M_index) {}
146*e4b17023SJohn Marino
147*e4b17023SJohn Marino template<typename _Tp>
148*e4b17023SJohn Marino inline gslice_array<_Tp>&
149*e4b17023SJohn Marino gslice_array<_Tp>::operator=(const gslice_array<_Tp>& __a)
150*e4b17023SJohn Marino {
151*e4b17023SJohn Marino std::__valarray_copy(_Array<_Tp>(__a._M_array),
152*e4b17023SJohn Marino _Array<size_t>(__a._M_index), _M_index.size(),
153*e4b17023SJohn Marino _M_array, _Array<size_t>(_M_index));
154*e4b17023SJohn Marino return *this;
155*e4b17023SJohn Marino }
156*e4b17023SJohn Marino
157*e4b17023SJohn Marino template<typename _Tp>
158*e4b17023SJohn Marino inline void
159*e4b17023SJohn Marino gslice_array<_Tp>::operator=(const _Tp& __t) const
160*e4b17023SJohn Marino {
161*e4b17023SJohn Marino std::__valarray_fill(_M_array, _Array<size_t>(_M_index),
162*e4b17023SJohn Marino _M_index.size(), __t);
163*e4b17023SJohn Marino }
164*e4b17023SJohn Marino
165*e4b17023SJohn Marino template<typename _Tp>
166*e4b17023SJohn Marino inline void
167*e4b17023SJohn Marino gslice_array<_Tp>::operator=(const valarray<_Tp>& __v) const
168*e4b17023SJohn Marino {
169*e4b17023SJohn Marino std::__valarray_copy(_Array<_Tp>(__v), __v.size(),
170*e4b17023SJohn Marino _M_array, _Array<size_t>(_M_index));
171*e4b17023SJohn Marino }
172*e4b17023SJohn Marino
173*e4b17023SJohn Marino template<typename _Tp>
174*e4b17023SJohn Marino template<class _Dom>
175*e4b17023SJohn Marino inline void
176*e4b17023SJohn Marino gslice_array<_Tp>::operator=(const _Expr<_Dom, _Tp>& __e) const
177*e4b17023SJohn Marino {
178*e4b17023SJohn Marino std::__valarray_copy (__e, _M_index.size(), _M_array,
179*e4b17023SJohn Marino _Array<size_t>(_M_index));
180*e4b17023SJohn Marino }
181*e4b17023SJohn Marino
182*e4b17023SJohn Marino #undef _DEFINE_VALARRAY_OPERATOR
183*e4b17023SJohn Marino #define _DEFINE_VALARRAY_OPERATOR(_Op, _Name) \
184*e4b17023SJohn Marino template<typename _Tp> \
185*e4b17023SJohn Marino inline void \
186*e4b17023SJohn Marino gslice_array<_Tp>::operator _Op##=(const valarray<_Tp>& __v) const \
187*e4b17023SJohn Marino { \
188*e4b17023SJohn Marino _Array_augmented_##_Name(_M_array, _Array<size_t>(_M_index), \
189*e4b17023SJohn Marino _Array<_Tp>(__v), __v.size()); \
190*e4b17023SJohn Marino } \
191*e4b17023SJohn Marino \
192*e4b17023SJohn Marino template<typename _Tp> \
193*e4b17023SJohn Marino template<class _Dom> \
194*e4b17023SJohn Marino inline void \
195*e4b17023SJohn Marino gslice_array<_Tp>::operator _Op##= (const _Expr<_Dom, _Tp>& __e) const\
196*e4b17023SJohn Marino { \
197*e4b17023SJohn Marino _Array_augmented_##_Name(_M_array, _Array<size_t>(_M_index), __e,\
198*e4b17023SJohn Marino _M_index.size()); \
199*e4b17023SJohn Marino }
200*e4b17023SJohn Marino
201*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(*, __multiplies)
202*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(/, __divides)
203*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(%, __modulus)
204*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(+, __plus)
205*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(-, __minus)
206*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(^, __bitwise_xor)
207*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(&, __bitwise_and)
208*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(|, __bitwise_or)
209*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(<<, __shift_left)
210*e4b17023SJohn Marino _DEFINE_VALARRAY_OPERATOR(>>, __shift_right)
211*e4b17023SJohn Marino
212*e4b17023SJohn Marino #undef _DEFINE_VALARRAY_OPERATOR
213*e4b17023SJohn Marino
214*e4b17023SJohn Marino // @} group numeric_arrays
215*e4b17023SJohn Marino
216*e4b17023SJohn Marino _GLIBCXX_END_NAMESPACE_VERSION
217*e4b17023SJohn Marino } // namespace
218*e4b17023SJohn Marino
219*e4b17023SJohn Marino #endif /* _GSLICE_ARRAY_H */
220