00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038 #ifndef _CPP_BITS_GSLICE_ARRAY
00039 #define _CPP_BITS_GSLICE_ARRAY 1
00040
00041 #pragma GCC system_header
00042
00043 namespace std {
00044
00045 template<typename _Tp>
00046 class gslice_array
00047 {
00048 public:
00049 typedef _Tp value_type;
00050
00051 gslice_array(const gslice_array&);
00052 gslice_array& operator=(const gslice_array&);
00053
00054 void operator=(const valarray<_Tp>&) const;
00055 void operator*=(const valarray<_Tp>&) const;
00056 void operator/=(const valarray<_Tp>&) const;
00057 void operator%=(const valarray<_Tp>&) const;
00058 void operator+=(const valarray<_Tp>&) const;
00059 void operator-=(const valarray<_Tp>&) const;
00060 void operator^=(const valarray<_Tp>&) const;
00061 void operator&=(const valarray<_Tp>&) const;
00062 void operator|=(const valarray<_Tp>&) const;
00063 void operator<<=(const valarray<_Tp>&) const;
00064 void operator>>=(const valarray<_Tp>&) const;
00065 void operator=(const _Tp&) const;
00066
00067 template<class _Dom>
00068 void operator=(const _Expr<_Dom,_Tp>&) const;
00069 template<class _Dom>
00070 void operator*=(const _Expr<_Dom,_Tp>&) const;
00071 template<class _Dom>
00072 void operator/=(const _Expr<_Dom,_Tp>&) const;
00073 template<class _Dom>
00074 void operator%=(const _Expr<_Dom,_Tp>&) const;
00075 template<class _Dom>
00076 void operator+=(const _Expr<_Dom,_Tp>&) const;
00077 template<class _Dom>
00078 void operator-=(const _Expr<_Dom,_Tp>&) const;
00079 template<class _Dom>
00080 void operator^=(const _Expr<_Dom,_Tp>&) const;
00081 template<class _Dom>
00082 void operator&=(const _Expr<_Dom,_Tp>&) const;
00083 template<class _Dom>
00084 void operator|=(const _Expr<_Dom,_Tp>&) const;
00085 template<class _Dom>
00086 void operator<<=(const _Expr<_Dom,_Tp>&) const;
00087 template<class _Dom>
00088 void operator>>=(const _Expr<_Dom,_Tp>&) const;
00089
00090 private:
00091 _Array<_Tp> _M_array;
00092 const valarray<size_t>& _M_index;
00093
00094 friend class valarray<_Tp>;
00095
00096 gslice_array(_Array<_Tp>, const valarray<size_t>&);
00097
00098
00099 gslice_array();
00100 };
00101
00102 template<typename _Tp>
00103 inline
00104 gslice_array<_Tp>::gslice_array(_Array<_Tp> __a,
00105 const valarray<size_t>& __i)
00106 : _M_array(__a), _M_index(__i) {}
00107
00108
00109 template<typename _Tp>
00110 inline
00111 gslice_array<_Tp>::gslice_array(const gslice_array<_Tp>& __a)
00112 : _M_array(__a._M_array), _M_index(__a._M_index) {}
00113
00114 template<typename _Tp>
00115 inline gslice_array<_Tp>&
00116 gslice_array<_Tp>::operator=(const gslice_array<_Tp>& __a)
00117 {
00118 __valarray_copy(_Array<_Tp>(__a._M_array),
00119 _Array<size_t>(__a._M_index), _M_index.size(),
00120 _M_array, _Array<size_t>(_M_index));
00121 return *this;
00122 }
00123
00124 template<typename _Tp>
00125 inline void
00126 gslice_array<_Tp>::operator=(const _Tp& __t) const
00127 {
00128 __valarray_fill(_M_array, _Array<size_t>(_M_index),
00129 _M_index.size(), __t);
00130 }
00131
00132 template<typename _Tp>
00133 inline void
00134 gslice_array<_Tp>::operator=(const valarray<_Tp>& __v) const
00135 {
00136 __valarray_copy(_Array<_Tp>(__v), __v.size(),
00137 _M_array, _Array<size_t>(_M_index));
00138 }
00139
00140 template<typename _Tp>
00141 template<class _Dom>
00142 inline void
00143 gslice_array<_Tp>::operator=(const _Expr<_Dom, _Tp>& __e) const
00144 {
00145 __valarray_copy (__e, _M_index.size(), _M_array,
00146 _Array<size_t>(_M_index));
00147 }
00148
00149 #undef _DEFINE_VALARRAY_OPERATOR
00150 #define _DEFINE_VALARRAY_OPERATOR(_Op, _Name) \
00151 template<typename _Tp> \
00152 inline void \
00153 gslice_array<_Tp>::operator _Op##=(const valarray<_Tp>& __v) const \
00154 { \
00155 _Array_augmented_##_Name(_M_array, _Array<size_t>(_M_index), \
00156 _Array<_Tp>(__v), __v.size()); \
00157 } \
00158 \
00159 template<typename _Tp> \
00160 template<class _Dom> \
00161 inline void \
00162 gslice_array<_Tp>::operator _Op##= (const _Expr<_Dom, _Tp>& __e) const\
00163 { \
00164 _Array_augmented_##_Name(_M_array, _Array<size_t>(_M_index), __e,\
00165 _M_index.size()); \
00166 }
00167
00168 _DEFINE_VALARRAY_OPERATOR(*, __multiplies)
00169 _DEFINE_VALARRAY_OPERATOR(/, __divides)
00170 _DEFINE_VALARRAY_OPERATOR(%, __modulus)
00171 _DEFINE_VALARRAY_OPERATOR(+, __plus)
00172 _DEFINE_VALARRAY_OPERATOR(-, __minus)
00173 _DEFINE_VALARRAY_OPERATOR(^, __bitwise_xor)
00174 _DEFINE_VALARRAY_OPERATOR(&, __bitwise_and)
00175 _DEFINE_VALARRAY_OPERATOR(|, __bitwise_or)
00176 _DEFINE_VALARRAY_OPERATOR(<<, __shift_left)
00177 _DEFINE_VALARRAY_OPERATOR(>>, __shift_right)
00178
00179 #undef _DEFINE_VALARRAY_OPERATOR
00180
00181 }
00182
00183 #endif
00184
00185
00186
00187