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 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