libstdc++
slice_array.h
Go to the documentation of this file.
1 // The template and inlines for the -*- C++ -*- slice_array class.
2 
3 // Copyright (C) 1997, 1998, 1999, 2000, 2001, 2002, 2004, 2005, 2006, 2009,
4 // 2010, 2011 Free Software Foundation, Inc.
5 //
6 // This file is part of the GNU ISO C++ Library. This library is free
7 // software; you can redistribute it and/or modify it under the
8 // terms of the GNU General Public License as published by the
9 // Free Software Foundation; either version 3, or (at your option)
10 // any later version.
11 
12 // This library is distributed in the hope that it will be useful,
13 // but WITHOUT ANY WARRANTY; without even the implied warranty of
14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 // GNU General Public License for more details.
16 
17 // Under Section 7 of GPL version 3, you are granted additional
18 // permissions described in the GCC Runtime Library Exception, version
19 // 3.1, as published by the Free Software Foundation.
20 
21 // You should have received a copy of the GNU General Public License and
22 // a copy of the GCC Runtime Library Exception along with this program;
23 // see the files COPYING3 and COPYING.RUNTIME respectively. If not, see
24 // <http://www.gnu.org/licenses/>.
25 
26 /** @file bits/slice_array.h
27  * This is an internal header file, included by other library headers.
28  * Do not attempt to use it directly. @headername{valarray}
29  */
30 
31 // Written by Gabriel Dos Reis <Gabriel.Dos-Reis@DPTMaths.ENS-Cachan.Fr>
32 
33 #ifndef _SLICE_ARRAY_H
34 #define _SLICE_ARRAY_H 1
35 
36 #pragma GCC system_header
37 
38 namespace std _GLIBCXX_VISIBILITY(default)
39 {
40 _GLIBCXX_BEGIN_NAMESPACE_VERSION
41 
42  /**
43  * @addtogroup numeric_arrays
44  * @{
45  */
46 
47  /**
48  * @brief Class defining one-dimensional subset of an array.
49  *
50  * The slice class represents a one-dimensional subset of an array,
51  * specified by three parameters: start offset, size, and stride. The
52  * start offset is the index of the first element of the array that is part
53  * of the subset. The size is the total number of elements in the subset.
54  * Stride is the distance between each successive array element to include
55  * in the subset.
56  *
57  * For example, with an array of size 10, and a slice with offset 1, size 3
58  * and stride 2, the subset consists of array elements 1, 3, and 5.
59  */
60  class slice
61  {
62  public:
63  /// Construct an empty slice.
64  slice();
65 
66  /**
67  * @brief Construct a slice.
68  *
69  * @param __o Offset in array of first element.
70  * @param __d Number of elements in slice.
71  * @param __s Stride between array elements.
72  */
73  slice(size_t __o, size_t __d, size_t __s);
74 
75  /// Return array offset of first slice element.
76  size_t start() const;
77  /// Return size of slice.
78  size_t size() const;
79  /// Return array stride of slice.
80  size_t stride() const;
81 
82  private:
83  size_t _M_off; // offset
84  size_t _M_sz; // size
85  size_t _M_st; // stride unit
86  };
87 
88  // _GLIBCXX_RESOLVE_LIB_DEFECTS
89  // 543. valarray slice default constructor
90  inline
92  : _M_off(0), _M_sz(0), _M_st(0) {}
93 
94  inline
95  slice::slice(size_t __o, size_t __d, size_t __s)
96  : _M_off(__o), _M_sz(__d), _M_st(__s) {}
97 
98  inline size_t
99  slice::start() const
100  { return _M_off; }
101 
102  inline size_t
103  slice::size() const
104  { return _M_sz; }
105 
106  inline size_t
108  { return _M_st; }
109 
110  /**
111  * @brief Reference to one-dimensional subset of an array.
112  *
113  * A slice_array is a reference to the actual elements of an array
114  * specified by a slice. The way to get a slice_array is to call
115  * operator[](slice) on a valarray. The returned slice_array then permits
116  * carrying operations out on the referenced subset of elements in the
117  * original valarray. For example, operator+=(valarray) will add values
118  * to the subset of elements in the underlying valarray this slice_array
119  * refers to.
120  *
121  * @param Tp Element type.
122  */
123  template<typename _Tp>
125  {
126  public:
127  typedef _Tp value_type;
128 
129  // _GLIBCXX_RESOLVE_LIB_DEFECTS
130  // 253. valarray helper functions are almost entirely useless
131 
132  /// Copy constructor. Both slices refer to the same underlying array.
133  slice_array(const slice_array&);
134 
135  /// Assignment operator. Assigns slice elements to corresponding
136  /// elements of @a a.
138 
139  /// Assign slice elements to corresponding elements of @a v.
140  void operator=(const valarray<_Tp>&) const;
141  /// Multiply slice elements by corresponding elements of @a v.
142  void operator*=(const valarray<_Tp>&) const;
143  /// Divide slice elements by corresponding elements of @a v.
144  void operator/=(const valarray<_Tp>&) const;
145  /// Modulo slice elements by corresponding elements of @a v.
146  void operator%=(const valarray<_Tp>&) const;
147  /// Add corresponding elements of @a v to slice elements.
148  void operator+=(const valarray<_Tp>&) const;
149  /// Subtract corresponding elements of @a v from slice elements.
150  void operator-=(const valarray<_Tp>&) const;
151  /// Logical xor slice elements with corresponding elements of @a v.
152  void operator^=(const valarray<_Tp>&) const;
153  /// Logical and slice elements with corresponding elements of @a v.
154  void operator&=(const valarray<_Tp>&) const;
155  /// Logical or slice elements with corresponding elements of @a v.
156  void operator|=(const valarray<_Tp>&) const;
157  /// Left shift slice elements by corresponding elements of @a v.
158  void operator<<=(const valarray<_Tp>&) const;
159  /// Right shift slice elements by corresponding elements of @a v.
160  void operator>>=(const valarray<_Tp>&) const;
161  /// Assign all slice elements to @a t.
162  void operator=(const _Tp &) const;
163  // ~slice_array ();
164 
165  template<class _Dom>
166  void operator=(const _Expr<_Dom, _Tp>&) const;
167  template<class _Dom>
168  void operator*=(const _Expr<_Dom, _Tp>&) const;
169  template<class _Dom>
170  void operator/=(const _Expr<_Dom, _Tp>&) const;
171  template<class _Dom>
172  void operator%=(const _Expr<_Dom, _Tp>&) const;
173  template<class _Dom>
174  void operator+=(const _Expr<_Dom, _Tp>&) const;
175  template<class _Dom>
176  void operator-=(const _Expr<_Dom, _Tp>&) const;
177  template<class _Dom>
178  void operator^=(const _Expr<_Dom, _Tp>&) const;
179  template<class _Dom>
180  void operator&=(const _Expr<_Dom, _Tp>&) const;
181  template<class _Dom>
182  void operator|=(const _Expr<_Dom, _Tp>&) const;
183  template<class _Dom>
184  void operator<<=(const _Expr<_Dom, _Tp>&) const;
185  template<class _Dom>
186  void operator>>=(const _Expr<_Dom, _Tp>&) const;
187 
188  private:
189  friend class valarray<_Tp>;
190  slice_array(_Array<_Tp>, const slice&);
191 
192  const size_t _M_sz;
193  const size_t _M_stride;
194  const _Array<_Tp> _M_array;
195 
196  // not implemented
197  slice_array();
198  };
199 
200  template<typename _Tp>
201  inline
202  slice_array<_Tp>::slice_array(_Array<_Tp> __a, const slice& __s)
203  : _M_sz(__s.size()), _M_stride(__s.stride()),
204  _M_array(__a.begin() + __s.start()) {}
205 
206  template<typename _Tp>
207  inline
209  : _M_sz(a._M_sz), _M_stride(a._M_stride), _M_array(a._M_array) {}
210 
211  // template<typename _Tp>
212  // inline slice_array<_Tp>::~slice_array () {}
213 
214  template<typename _Tp>
215  inline slice_array<_Tp>&
217  {
218  std::__valarray_copy(__a._M_array, __a._M_sz, __a._M_stride,
219  _M_array, _M_stride);
220  return *this;
221  }
222 
223  template<typename _Tp>
224  inline void
225  slice_array<_Tp>::operator=(const _Tp& __t) const
226  { std::__valarray_fill(_M_array, _M_sz, _M_stride, __t); }
227 
228  template<typename _Tp>
229  inline void
231  { std::__valarray_copy(_Array<_Tp>(__v), _M_array, _M_sz, _M_stride); }
232 
233  template<typename _Tp>
234  template<class _Dom>
235  inline void
236  slice_array<_Tp>::operator=(const _Expr<_Dom,_Tp>& __e) const
237  { std::__valarray_copy(__e, _M_sz, _M_array, _M_stride); }
238 
239 #undef _DEFINE_VALARRAY_OPERATOR
240 #define _DEFINE_VALARRAY_OPERATOR(_Op,_Name) \
241  template<typename _Tp> \
242  inline void \
243  slice_array<_Tp>::operator _Op##=(const valarray<_Tp>& __v) const \
244  { \
245  _Array_augmented_##_Name(_M_array, _M_sz, _M_stride, _Array<_Tp>(__v));\
246  } \
247  \
248  template<typename _Tp> \
249  template<class _Dom> \
250  inline void \
251  slice_array<_Tp>::operator _Op##=(const _Expr<_Dom,_Tp>& __e) const\
252  { \
253  _Array_augmented_##_Name(_M_array, _M_stride, __e, _M_sz); \
254  }
255 
256 
257 _DEFINE_VALARRAY_OPERATOR(*, __multiplies)
258 _DEFINE_VALARRAY_OPERATOR(/, __divides)
259 _DEFINE_VALARRAY_OPERATOR(%, __modulus)
260 _DEFINE_VALARRAY_OPERATOR(+, __plus)
261 _DEFINE_VALARRAY_OPERATOR(-, __minus)
262 _DEFINE_VALARRAY_OPERATOR(^, __bitwise_xor)
263 _DEFINE_VALARRAY_OPERATOR(&, __bitwise_and)
264 _DEFINE_VALARRAY_OPERATOR(|, __bitwise_or)
265 _DEFINE_VALARRAY_OPERATOR(<<, __shift_left)
266 _DEFINE_VALARRAY_OPERATOR(>>, __shift_right)
267 
268 #undef _DEFINE_VALARRAY_OPERATOR
269 
270  // @} group numeric_arrays
271 
272 _GLIBCXX_END_NAMESPACE_VERSION
273 } // namespace
274 
275 #endif /* _SLICE_ARRAY_H */