1 // The template and inlines for the -*- C++ -*- valarray class.
3 // Copyright (C) 1997, 1998, 1999, 2000, 2001, 2002, 2004
4 // Free Software Foundation, Inc.
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 2, or (at your option)
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.
17 // You should have received a copy of the GNU General Public License along
18 // with this library; see the file COPYING. If not, write to the Free
19 // Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307,
22 // As a special exception, you may use this file as part of a free software
23 // library without restriction. Specifically, if other files instantiate
24 // templates or use macros or inline functions from this file, or you compile
25 // this file and link it with other files to produce an executable, this
26 // file does not by itself cause the resulting executable to be covered by
27 // the GNU General Public License. This exception does not however
28 // invalidate any other reasons why the executable file might be covered by
29 // the GNU General Public License.
31 // Written by Gabriel Dos Reis <Gabriel.Dos-Reis@DPTMaths.ENS-Cachan.Fr>
34 * This is a Standard C++ Library header. You should @c #include this header
35 * in your programs, rather than any of the "st[dl]_*.h" implementation files.
38 #ifndef _GLIBCXX_VALARRAY
39 #define _GLIBCXX_VALARRAY 1
41 #pragma GCC system_header
43 #include <bits/c++config.h>
49 #include <debug/debug.h>
53 template<class _Clos
, typename _Tp
>
56 template<typename _Tp1
, typename _Tp2
>
59 template<class _Oper
, template<class, class> class _Meta
, class _Dom
>
63 template<class, class> class _Meta1
,
64 template<class, class> class _Meta2
,
65 class _Dom1
, class _Dom2
>
68 template<template<class, class> class _Meta
, class _Dom
>
71 template<template<class, class> class _Meta
, class _Dom
>
74 template<template<class, class> class _Meta
, class _Dom
>
77 template<template<class, class> class _Meta
, class _Dom
>
80 template<template<class, class> class _Meta
, class _Dom
>
83 template<class _Tp
> class valarray
; // An array of type _Tp
84 class slice
; // BLAS-like slice out of an array
85 template<class _Tp
> class slice_array
;
86 class gslice
; // generalized slice out of an array
87 template<class _Tp
> class gslice_array
;
88 template<class _Tp
> class mask_array
; // masked array
89 template<class _Tp
> class indirect_array
; // indirected array
93 #include <bits/valarray_array.h>
94 #include <bits/valarray_before.h>
99 * @brief Smart array designed to support numeric processing.
101 * A valarray is an array that provides constraints intended to allow for
102 * effective optimization of numeric array processing by reducing the
103 * aliasing that can result from pointer representations. It represents a
104 * one-dimensional array from which different multidimensional subsets can
105 * be accessed and modified.
107 * @param Tp Type of object in the array.
115 typedef typename __fun
<_Op
, _Tp
>::result_type __rt
;
116 typedef _Expr
<_UnClos
<_Op
, _ValArray
, _Tp
>, __rt
> _Rt
;
119 typedef _Tp value_type
;
121 // _lib.valarray.cons_ construct/destroy:
122 /// Construct an empty array.
125 /// Construct an array with @a n elements.
126 explicit valarray(size_t);
128 /// Construct an array with @a n elements initialized to @a t.
129 valarray(const _Tp
&, size_t);
131 /// Construct an array initialized to the first @a n elements of @a t.
132 valarray(const _Tp
* __restrict__
, size_t);
134 /// Copy constructor.
135 valarray(const valarray
&);
137 /// Construct an array with the same size and values in @a sa.
138 valarray(const slice_array
<_Tp
>&);
140 /// Construct an array with the same size and values in @a ga.
141 valarray(const gslice_array
<_Tp
>&);
143 /// Construct an array with the same size and values in @a ma.
144 valarray(const mask_array
<_Tp
>&);
146 /// Construct an array with the same size and values in @a ia.
147 valarray(const indirect_array
<_Tp
>&);
150 valarray(const _Expr
<_Dom
,_Tp
>& __e
);
153 // _lib.valarray.assign_ assignment:
155 * @brief Assign elements to an array.
157 * Assign elements of array to values in @a v. Results are undefined
158 * if @a v is not the same size as this array.
160 * @param v Valarray to get values from.
162 valarray
<_Tp
>& operator=(const valarray
<_Tp
>&);
165 * @brief Assign elements to a value.
167 * Assign all elements of array to @a t.
169 * @param t Value for elements.
171 valarray
<_Tp
>& operator=(const _Tp
&);
174 * @brief Assign elements to an array subset.
176 * Assign elements of array to values in @a sa. Results are undefined
177 * if @a sa is not the same size as this array.
179 * @param sa Array slice to get values from.
181 valarray
<_Tp
>& operator=(const slice_array
<_Tp
>&);
184 * @brief Assign elements to an array subset.
186 * Assign elements of array to values in @a ga. Results are undefined
187 * if @a ga is not the same size as this array.
189 * @param ga Array slice to get values from.
191 valarray
<_Tp
>& operator=(const gslice_array
<_Tp
>&);
194 * @brief Assign elements to an array subset.
196 * Assign elements of array to values in @a ma. Results are undefined
197 * if @a ma is not the same size as this array.
199 * @param ma Array slice to get values from.
201 valarray
<_Tp
>& operator=(const mask_array
<_Tp
>&);
204 * @brief Assign elements to an array subset.
206 * Assign elements of array to values in @a ia. Results are undefined
207 * if @a ia is not the same size as this array.
209 * @param ia Array slice to get values from.
211 valarray
<_Tp
>& operator=(const indirect_array
<_Tp
>&);
213 template<class _Dom
> valarray
<_Tp
>&
214 operator= (const _Expr
<_Dom
,_Tp
>&);
216 // _lib.valarray.access_ element access:
218 * Return a reference to the i'th array element.
220 * @param i Index of element to return.
221 * @return Reference to the i'th element.
223 _Tp
& operator[](size_t);
225 // _GLIBCXX_RESOLVE_LIB_DEFECTS
226 // 389. Const overload of valarray::operator[] returns by value.
227 const _Tp
& operator[](size_t) const;
229 // _lib.valarray.sub_ subset operations:
231 * @brief Return an array subset.
233 * Returns a new valarray containing the elements of the array
234 * indicated by the slice argument. The new valarray is the size of
235 * the input slice. @see slice.
237 * @param s The source slice.
238 * @return New valarray containing elements in @a s.
240 _Expr
<_SClos
<_ValArray
,_Tp
>, _Tp
> operator[](slice
) const;
243 * @brief Return a reference to an array subset.
245 * Returns a new valarray containing the elements of the array
246 * indicated by the slice argument. The new valarray is the size of
247 * the input slice. @see slice.
249 * @param s The source slice.
250 * @return New valarray containing elements in @a s.
252 slice_array
<_Tp
> operator[](slice
);
255 * @brief Return an array subset.
257 * Returns a slice_array referencing the elements of the array
258 * indicated by the slice argument. @see gslice.
260 * @param s The source slice.
261 * @return Slice_array referencing elements indicated by @a s.
263 _Expr
<_GClos
<_ValArray
,_Tp
>, _Tp
> operator[](const gslice
&) const;
266 * @brief Return a reference to an array subset.
268 * Returns a new valarray containing the elements of the array
269 * indicated by the gslice argument. The new valarray is
270 * the size of the input gslice. @see gslice.
272 * @param s The source gslice.
273 * @return New valarray containing elements in @a s.
275 gslice_array
<_Tp
> operator[](const gslice
&);
278 * @brief Return an array subset.
280 * Returns a new valarray containing the elements of the array
281 * indicated by the argument. The input is a valarray of bool which
282 * represents a bitmask indicating which elements should be copied into
283 * the new valarray. Each element of the array is added to the return
284 * valarray if the corresponding element of the argument is true.
286 * @param m The valarray bitmask.
287 * @return New valarray containing elements indicated by @a m.
289 valarray
<_Tp
> operator[](const valarray
<bool>&) const;
292 * @brief Return a reference to an array subset.
294 * Returns a new mask_array referencing the elements of the array
295 * indicated by the argument. The input is a valarray of bool which
296 * represents a bitmask indicating which elements are part of the
297 * subset. Elements of the array are part of the subset if the
298 * corresponding element of the argument is true.
300 * @param m The valarray bitmask.
301 * @return New valarray containing elements indicated by @a m.
303 mask_array
<_Tp
> operator[](const valarray
<bool>&);
306 * @brief Return an array subset.
308 * Returns a new valarray containing the elements of the array
309 * indicated by the argument. The elements in the argument are
310 * interpreted as the indices of elements of this valarray to copy to
311 * the return valarray.
313 * @param i The valarray element index list.
314 * @return New valarray containing elements in @a s.
316 _Expr
<_IClos
<_ValArray
, _Tp
>, _Tp
>
317 operator[](const valarray
<size_t>&) const;
320 * @brief Return a reference to an array subset.
322 * Returns an indirect_array referencing the elements of the array
323 * indicated by the argument. The elements in the argument are
324 * interpreted as the indices of elements of this valarray to include
325 * in the subset. The returned indirect_array refers to these
328 * @param i The valarray element index list.
329 * @return Indirect_array referencing elements in @a i.
331 indirect_array
<_Tp
> operator[](const valarray
<size_t>&);
333 // _lib.valarray.unary_ unary operators:
334 /// Return a new valarray by applying unary + to each element.
335 typename _UnaryOp
<__unary_plus
>::_Rt
operator+() const;
337 /// Return a new valarray by applying unary - to each element.
338 typename _UnaryOp
<__negate
>::_Rt
operator-() const;
340 /// Return a new valarray by applying unary ~ to each element.
341 typename _UnaryOp
<__bitwise_not
>::_Rt
operator~() const;
343 /// Return a new valarray by applying unary ! to each element.
344 typename _UnaryOp
<__logical_not
>::_Rt
operator!() const;
346 // _lib.valarray.cassign_ computed assignment:
347 /// Multiply each element of array by @a t.
348 valarray
<_Tp
>& operator*=(const _Tp
&);
350 /// Divide each element of array by @a t.
351 valarray
<_Tp
>& operator/=(const _Tp
&);
353 /// Set each element e of array to e % @a t.
354 valarray
<_Tp
>& operator%=(const _Tp
&);
356 /// Add @a t to each element of array.
357 valarray
<_Tp
>& operator+=(const _Tp
&);
359 /// Subtract @a t to each element of array.
360 valarray
<_Tp
>& operator-=(const _Tp
&);
362 /// Set each element e of array to e ^ @a t.
363 valarray
<_Tp
>& operator^=(const _Tp
&);
365 /// Set each element e of array to e & @a t.
366 valarray
<_Tp
>& operator&=(const _Tp
&);
368 /// Set each element e of array to e | @a t.
369 valarray
<_Tp
>& operator|=(const _Tp
&);
371 /// Left shift each element e of array by @a t bits.
372 valarray
<_Tp
>& operator<<=(const _Tp
&);
374 /// Right shift each element e of array by @a t bits.
375 valarray
<_Tp
>& operator>>=(const _Tp
&);
377 /// Multiply elements of array by corresponding elements of @a v.
378 valarray
<_Tp
>& operator*=(const valarray
<_Tp
>&);
380 /// Divide elements of array by corresponding elements of @a v.
381 valarray
<_Tp
>& operator/=(const valarray
<_Tp
>&);
383 /// Modulo elements of array by corresponding elements of @a v.
384 valarray
<_Tp
>& operator%=(const valarray
<_Tp
>&);
386 /// Add corresponding elements of @a v to elements of array.
387 valarray
<_Tp
>& operator+=(const valarray
<_Tp
>&);
389 /// Subtract corresponding elements of @a v from elements of array.
390 valarray
<_Tp
>& operator-=(const valarray
<_Tp
>&);
392 /// Logical xor corresponding elements of @a v with elements of array.
393 valarray
<_Tp
>& operator^=(const valarray
<_Tp
>&);
395 /// Logical or corresponding elements of @a v with elements of array.
396 valarray
<_Tp
>& operator|=(const valarray
<_Tp
>&);
398 /// Logical and corresponding elements of @a v with elements of array.
399 valarray
<_Tp
>& operator&=(const valarray
<_Tp
>&);
401 /// Left shift elements of array by corresponding elements of @a v.
402 valarray
<_Tp
>& operator<<=(const valarray
<_Tp
>&);
404 /// Right shift elements of array by corresponding elements of @a v.
405 valarray
<_Tp
>& operator>>=(const valarray
<_Tp
>&);
408 valarray
<_Tp
>& operator*=(const _Expr
<_Dom
,_Tp
>&);
410 valarray
<_Tp
>& operator/=(const _Expr
<_Dom
,_Tp
>&);
412 valarray
<_Tp
>& operator%=(const _Expr
<_Dom
,_Tp
>&);
414 valarray
<_Tp
>& operator+=(const _Expr
<_Dom
,_Tp
>&);
416 valarray
<_Tp
>& operator-=(const _Expr
<_Dom
,_Tp
>&);
418 valarray
<_Tp
>& operator^=(const _Expr
<_Dom
,_Tp
>&);
420 valarray
<_Tp
>& operator|=(const _Expr
<_Dom
,_Tp
>&);
422 valarray
<_Tp
>& operator&=(const _Expr
<_Dom
,_Tp
>&);
424 valarray
<_Tp
>& operator<<=(const _Expr
<_Dom
,_Tp
>&);
426 valarray
<_Tp
>& operator>>=(const _Expr
<_Dom
,_Tp
>&);
429 // _lib.valarray.members_ member functions:
430 /// Return the number of elements in array.
434 * @brief Return the sum of all elements in the array.
436 * Accumulates the sum of all elements into a Tp using +=. The order
437 * of adding the elements is unspecified.
441 /// Return the minimum element using operator<().
444 /// Return the maximum element using operator<().
447 // // FIXME: Extension
448 // _Tp product () const;
451 * @brief Return a shifted array.
453 * A new valarray is constructed as a copy of this array with elements
454 * in shifted positions. For an element with index i, the new position
455 * is i - n. The new valarray is the same size as the current one.
456 * New elements without a value are set to 0. Elements whos new
457 * position is outside the bounds of the array are discarded.
459 * Positive arguments shift toward index 0, discarding elements [0, n).
460 * Negative arguments discard elements from the top of the array.
462 * @param n Number of element positions to shift.
463 * @return New valarray with elements in shifted positions.
465 valarray
<_Tp
> shift (int) const;
468 * @brief Return a rotated array.
470 * A new valarray is constructed as a copy of this array with elements
471 * in shifted positions. For an element with index i, the new position
472 * is (i - n) % size(). The new valarray is the same size as the
473 * current one. Elements that are shifted beyond the array bounds are
474 * shifted into the other end of the array. No elements are lost.
476 * Positive arguments shift toward index 0, wrapping around the top.
477 * Negative arguments shift towards the top, wrapping around to 0.
479 * @param n Number of element positions to rotate.
480 * @return New valarray with elements in shifted positions.
482 valarray
<_Tp
> cshift(int) const;
485 * @brief Apply a function to the array.
487 * Returns a new valarray with elements assigned to the result of
488 * applying func to the corresponding element of this array. The new
489 * array is the same size as this one.
491 * @param func Function of Tp returning Tp to apply.
492 * @return New valarray with transformed elements.
494 _Expr
<_ValFunClos
<_ValArray
,_Tp
>,_Tp
> apply(_Tp
func(_Tp
)) const;
497 * @brief Apply a function to the array.
499 * Returns a new valarray with elements assigned to the result of
500 * applying func to the corresponding element of this array. The new
501 * array is the same size as this one.
503 * @param func Function of const Tp& returning Tp to apply.
504 * @return New valarray with transformed elements.
506 _Expr
<_RefFunClos
<_ValArray
,_Tp
>,_Tp
> apply(_Tp
func(const _Tp
&)) const;
509 * @brief Resize array.
511 * Resize this array to be @a size and set all elements to @a c. All
512 * references and iterators are invalidated.
514 * @param size New array size.
515 * @param c New value for all elements.
517 void resize(size_t __size
, _Tp __c
= _Tp());
521 _Tp
* __restrict__ _M_data
;
523 friend class _Array
<_Tp
>;
526 template<typename _Tp
>
528 valarray
<_Tp
>::operator[](size_t __i
) const
530 __glibcxx_requires_subscript(__i
);
534 template<typename _Tp
>
536 valarray
<_Tp
>::operator[](size_t __i
)
538 __glibcxx_requires_subscript(__i
);
544 #include <bits/valarray_after.h>
546 #include <bits/slice_array.h>
547 #include <bits/gslice.h>
548 #include <bits/gslice_array.h>
549 #include <bits/mask_array.h>
550 #include <bits/indirect_array.h>
554 template<typename _Tp
>
556 valarray
<_Tp
>::valarray() : _M_size(0), _M_data(0) {}
558 template<typename _Tp
>
560 valarray
<_Tp
>::valarray(size_t __n
)
561 : _M_size(__n
), _M_data(__valarray_get_storage
<_Tp
>(__n
))
562 { std::__valarray_default_construct(_M_data
, _M_data
+ __n
); }
564 template<typename _Tp
>
566 valarray
<_Tp
>::valarray(const _Tp
& __t
, size_t __n
)
567 : _M_size(__n
), _M_data(__valarray_get_storage
<_Tp
>(__n
))
568 { std::__valarray_fill_construct(_M_data
, _M_data
+ __n
, __t
); }
570 template<typename _Tp
>
572 valarray
<_Tp
>::valarray(const _Tp
* __restrict__ __p
, size_t __n
)
573 : _M_size(__n
), _M_data(__valarray_get_storage
<_Tp
>(__n
))
575 _GLIBCXX_DEBUG_ASSERT(__p
!= 0 || __n
== 0);
576 std::__valarray_copy_construct(__p
, __p
+ __n
, _M_data
);
579 template<typename _Tp
>
581 valarray
<_Tp
>::valarray(const valarray
<_Tp
>& __v
)
582 : _M_size(__v
._M_size
), _M_data(__valarray_get_storage
<_Tp
>(__v
._M_size
))
583 { std::__valarray_copy_construct(__v
._M_data
, __v
._M_data
+ _M_size
, _M_data
); }
585 template<typename _Tp
>
587 valarray
<_Tp
>::valarray(const slice_array
<_Tp
>& __sa
)
588 : _M_size(__sa
._M_sz
), _M_data(__valarray_get_storage
<_Tp
>(__sa
._M_sz
))
591 (__sa
._M_array
, __sa
._M_sz
, __sa
._M_stride
, _Array
<_Tp
>(_M_data
));
594 template<typename _Tp
>
596 valarray
<_Tp
>::valarray(const gslice_array
<_Tp
>& __ga
)
597 : _M_size(__ga
._M_index
.size()),
598 _M_data(__valarray_get_storage
<_Tp
>(_M_size
))
601 (__ga
._M_array
, _Array
<size_t>(__ga
._M_index
),
602 _Array
<_Tp
>(_M_data
), _M_size
);
605 template<typename _Tp
>
607 valarray
<_Tp
>::valarray(const mask_array
<_Tp
>& __ma
)
608 : _M_size(__ma
._M_sz
), _M_data(__valarray_get_storage
<_Tp
>(__ma
._M_sz
))
611 (__ma
._M_array
, __ma
._M_mask
, _Array
<_Tp
>(_M_data
), _M_size
);
614 template<typename _Tp
>
616 valarray
<_Tp
>::valarray(const indirect_array
<_Tp
>& __ia
)
617 : _M_size(__ia
._M_sz
), _M_data(__valarray_get_storage
<_Tp
>(__ia
._M_sz
))
620 (__ia
._M_array
, __ia
._M_index
, _Array
<_Tp
>(_M_data
), _M_size
);
623 template<typename _Tp
> template<class _Dom
>
625 valarray
<_Tp
>::valarray(const _Expr
<_Dom
, _Tp
>& __e
)
626 : _M_size(__e
.size()), _M_data(__valarray_get_storage
<_Tp
>(_M_size
))
627 { std::__valarray_copy(__e
, _M_size
, _Array
<_Tp
>(_M_data
)); }
629 template<typename _Tp
>
631 valarray
<_Tp
>::~valarray()
633 std::__valarray_destroy_elements(_M_data
, _M_data
+ _M_size
);
634 std::__valarray_release_memory(_M_data
);
637 template<typename _Tp
>
638 inline valarray
<_Tp
>&
639 valarray
<_Tp
>::operator=(const valarray
<_Tp
>& __v
)
641 _GLIBCXX_DEBUG_ASSERT(_M_size
== __v
._M_size
);
642 std::__valarray_copy(__v
._M_data
, _M_size
, _M_data
);
646 template<typename _Tp
>
647 inline valarray
<_Tp
>&
648 valarray
<_Tp
>::operator=(const _Tp
& __t
)
650 std::__valarray_fill(_M_data
, _M_size
, __t
);
654 template<typename _Tp
>
655 inline valarray
<_Tp
>&
656 valarray
<_Tp
>::operator=(const slice_array
<_Tp
>& __sa
)
658 _GLIBCXX_DEBUG_ASSERT(_M_size
== __sa
._M_sz
);
659 std::__valarray_copy(__sa
._M_array
, __sa
._M_sz
,
660 __sa
._M_stride
, _Array
<_Tp
>(_M_data
));
664 template<typename _Tp
>
665 inline valarray
<_Tp
>&
666 valarray
<_Tp
>::operator=(const gslice_array
<_Tp
>& __ga
)
668 _GLIBCXX_DEBUG_ASSERT(_M_size
== __ga
._M_index
.size());
669 std::__valarray_copy(__ga
._M_array
, _Array
<size_t>(__ga
._M_index
),
670 _Array
<_Tp
>(_M_data
), _M_size
);
674 template<typename _Tp
>
675 inline valarray
<_Tp
>&
676 valarray
<_Tp
>::operator=(const mask_array
<_Tp
>& __ma
)
678 _GLIBCXX_DEBUG_ASSERT(_M_size
== __ma
._M_sz
);
679 std::__valarray_copy(__ma
._M_array
, __ma
._M_mask
,
680 _Array
<_Tp
>(_M_data
), _M_size
);
684 template<typename _Tp
>
685 inline valarray
<_Tp
>&
686 valarray
<_Tp
>::operator=(const indirect_array
<_Tp
>& __ia
)
688 _GLIBCXX_DEBUG_ASSERT(_M_size
== __ia
._M_sz
);
689 std::__valarray_copy(__ia
._M_array
, __ia
._M_index
,
690 _Array
<_Tp
>(_M_data
), _M_size
);
694 template<typename _Tp
> template<class _Dom
>
695 inline valarray
<_Tp
>&
696 valarray
<_Tp
>::operator=(const _Expr
<_Dom
, _Tp
>& __e
)
698 _GLIBCXX_DEBUG_ASSERT(_M_size
== __e
.size());
699 std::__valarray_copy(__e
, _M_size
, _Array
<_Tp
>(_M_data
));
703 template<typename _Tp
>
704 inline _Expr
<_SClos
<_ValArray
,_Tp
>, _Tp
>
705 valarray
<_Tp
>::operator[](slice __s
) const
707 typedef _SClos
<_ValArray
,_Tp
> _Closure
;
708 return _Expr
<_Closure
, _Tp
>(_Closure (_Array
<_Tp
>(_M_data
), __s
));
711 template<typename _Tp
>
712 inline slice_array
<_Tp
>
713 valarray
<_Tp
>::operator[](slice __s
)
715 return slice_array
<_Tp
>(_Array
<_Tp
>(_M_data
), __s
);
718 template<typename _Tp
>
719 inline _Expr
<_GClos
<_ValArray
,_Tp
>, _Tp
>
720 valarray
<_Tp
>::operator[](const gslice
& __gs
) const
722 typedef _GClos
<_ValArray
,_Tp
> _Closure
;
723 return _Expr
<_Closure
, _Tp
>
724 (_Closure(_Array
<_Tp
>(_M_data
), __gs
._M_index
->_M_index
));
727 template<typename _Tp
>
728 inline gslice_array
<_Tp
>
729 valarray
<_Tp
>::operator[](const gslice
& __gs
)
731 return gslice_array
<_Tp
>
732 (_Array
<_Tp
>(_M_data
), __gs
._M_index
->_M_index
);
735 template<typename _Tp
>
737 valarray
<_Tp
>::operator[](const valarray
<bool>& __m
) const
740 size_t __e
= __m
.size();
741 for (size_t __i
=0; __i
<__e
; ++__i
)
743 return valarray
<_Tp
>(mask_array
<_Tp
>(_Array
<_Tp
>(_M_data
), __s
,
744 _Array
<bool> (__m
)));
747 template<typename _Tp
>
748 inline mask_array
<_Tp
>
749 valarray
<_Tp
>::operator[](const valarray
<bool>& __m
)
752 size_t __e
= __m
.size();
753 for (size_t __i
=0; __i
<__e
; ++__i
)
755 return mask_array
<_Tp
>(_Array
<_Tp
>(_M_data
), __s
, _Array
<bool>(__m
));
758 template<typename _Tp
>
759 inline _Expr
<_IClos
<_ValArray
,_Tp
>, _Tp
>
760 valarray
<_Tp
>::operator[](const valarray
<size_t>& __i
) const
762 typedef _IClos
<_ValArray
,_Tp
> _Closure
;
763 return _Expr
<_Closure
, _Tp
>(_Closure(*this, __i
));
766 template<typename _Tp
>
767 inline indirect_array
<_Tp
>
768 valarray
<_Tp
>::operator[](const valarray
<size_t>& __i
)
770 return indirect_array
<_Tp
>(_Array
<_Tp
>(_M_data
), __i
.size(),
771 _Array
<size_t>(__i
));
776 valarray
<_Tp
>::size() const
781 valarray
<_Tp
>::sum() const
783 _GLIBCXX_DEBUG_ASSERT(_M_size
> 0);
784 return std::__valarray_sum(_M_data
, _M_data
+ _M_size
);
787 // template<typename _Tp>
789 // valarray<_Tp>::product () const
791 // return __valarray_product(_M_data, _M_data + _M_size);
796 valarray
<_Tp
>::shift(int __n
) const
798 _Tp
* const __a
= static_cast<_Tp
*>
799 (__builtin_alloca(sizeof(_Tp
) * _M_size
));
800 if (__n
== 0) // no shift
801 std::__valarray_copy_construct(_M_data
, _M_data
+ _M_size
, __a
);
802 else if (__n
> 0) // __n > 0: shift left
804 if (size_t(__n
) > _M_size
)
805 std::__valarray_default_construct(__a
, __a
+ __n
);
808 std::__valarray_copy_construct(_M_data
+__n
, _M_data
+ _M_size
, __a
);
809 std::__valarray_default_construct(__a
+_M_size
-__n
, __a
+ _M_size
);
812 else // __n < 0: shift right
814 std::__valarray_copy_construct (_M_data
, _M_data
+_M_size
+__n
, __a
-__n
);
815 std::__valarray_default_construct(__a
, __a
- __n
);
817 return valarray
<_Tp
> (__a
, _M_size
);
822 valarray
<_Tp
>::cshift (int __n
) const
824 _Tp
* const __a
= static_cast<_Tp
*>
825 (__builtin_alloca (sizeof(_Tp
) * _M_size
));
826 if (__n
== 0) // no cshift
827 std::__valarray_copy_construct(_M_data
, _M_data
+ _M_size
, __a
);
828 else if (__n
> 0) // cshift left
830 std::__valarray_copy_construct(_M_data
, _M_data
+__n
, __a
+_M_size
-__n
);
831 std::__valarray_copy_construct(_M_data
+__n
, _M_data
+ _M_size
, __a
);
835 std::__valarray_copy_construct
836 (_M_data
+ _M_size
+__n
, _M_data
+ _M_size
, __a
);
837 std::__valarray_copy_construct
838 (_M_data
, _M_data
+ _M_size
+__n
, __a
- __n
);
840 return valarray
<_Tp
>(__a
, _M_size
);
845 valarray
<_Tp
>::resize (size_t __n
, _Tp __c
)
847 // This complication is so to make valarray<valarray<T> > work
848 // even though it is not required by the standard. Nobody should
849 // be saying valarray<valarray<T> > anyway. See the specs.
850 std::__valarray_destroy_elements(_M_data
, _M_data
+ _M_size
);
853 std::__valarray_release_memory(_M_data
);
855 _M_data
= __valarray_get_storage
<_Tp
>(__n
);
857 std::__valarray_fill_construct(_M_data
, _M_data
+ __n
, __c
);
860 template<typename _Tp
>
862 valarray
<_Tp
>::min() const
864 _GLIBCXX_DEBUG_ASSERT(_M_size
> 0);
865 return *std::min_element (_M_data
, _M_data
+_M_size
);
868 template<typename _Tp
>
870 valarray
<_Tp
>::max() const
872 _GLIBCXX_DEBUG_ASSERT(_M_size
> 0);
873 return *std::max_element (_M_data
, _M_data
+_M_size
);
877 inline _Expr
<_ValFunClos
<_ValArray
,_Tp
>,_Tp
>
878 valarray
<_Tp
>::apply(_Tp
func(_Tp
)) const
880 typedef _ValFunClos
<_ValArray
,_Tp
> _Closure
;
881 return _Expr
<_Closure
,_Tp
>(_Closure(*this, func
));
885 inline _Expr
<_RefFunClos
<_ValArray
,_Tp
>,_Tp
>
886 valarray
<_Tp
>::apply(_Tp
func(const _Tp
&)) const
888 typedef _RefFunClos
<_ValArray
,_Tp
> _Closure
;
889 return _Expr
<_Closure
,_Tp
>(_Closure(*this, func
));
892 #define _DEFINE_VALARRAY_UNARY_OPERATOR(_Op, _Name) \
893 template<typename _Tp> \
894 inline typename valarray<_Tp>::template _UnaryOp<_Name>::_Rt \
895 valarray<_Tp>::operator _Op() const \
897 typedef _UnClos<_Name,_ValArray,_Tp> _Closure; \
898 typedef typename __fun<_Name, _Tp>::result_type _Rt; \
899 return _Expr<_Closure, _Rt>(_Closure(*this)); \
902 _DEFINE_VALARRAY_UNARY_OPERATOR(+, __unary_plus
)
903 _DEFINE_VALARRAY_UNARY_OPERATOR(-, __negate
)
904 _DEFINE_VALARRAY_UNARY_OPERATOR(~, __bitwise_not
)
905 _DEFINE_VALARRAY_UNARY_OPERATOR (!, __logical_not
)
907 #undef _DEFINE_VALARRAY_UNARY_OPERATOR
909 #define _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(_Op, _Name) \
910 template<class _Tp> \
911 inline valarray<_Tp>& \
912 valarray<_Tp>::operator _Op##=(const _Tp &__t) \
914 _Array_augmented_##_Name(_Array<_Tp>(_M_data), _M_size, __t); \
918 template<class _Tp> \
919 inline valarray<_Tp>& \
920 valarray<_Tp>::operator _Op##=(const valarray<_Tp> &__v) \
922 _GLIBCXX_DEBUG_ASSERT(_M_size == __v._M_size); \
923 _Array_augmented_##_Name(_Array<_Tp>(_M_data), _M_size, \
924 _Array<_Tp>(__v._M_data)); \
928 _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(+, __plus
)
929 _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(-, __minus
)
930 _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(*, __multiplies
)
931 _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(/, __divides
)
932 _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(%, __modulus
)
933 _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(^, __bitwise_xor
)
934 _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(&, __bitwise_and
)
935 _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(|, __bitwise_or
)
936 _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(<<, __shift_left
)
937 _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(>>, __shift_right
)
939 #undef _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT
941 #define _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(_Op, _Name) \
942 template<class _Tp> template<class _Dom> \
943 inline valarray<_Tp>& \
944 valarray<_Tp>::operator _Op##=(const _Expr<_Dom,_Tp>& __e) \
946 _Array_augmented_##_Name(_Array<_Tp>(_M_data), __e, _M_size); \
950 _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(+, __plus
)
951 _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(-, __minus
)
952 _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(*, __multiplies
)
953 _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(/, __divides
)
954 _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(%, __modulus
)
955 _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(^, __bitwise_xor
)
956 _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(&, __bitwise_and
)
957 _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(|, __bitwise_or
)
958 _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(<<, __shift_left
)
959 _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(>>, __shift_right
)
961 #undef _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT
964 #define _DEFINE_BINARY_OPERATOR(_Op, _Name) \
965 template<typename _Tp> \
966 inline _Expr<_BinClos<_Name,_ValArray,_ValArray,_Tp,_Tp>, \
967 typename __fun<_Name, _Tp>::result_type> \
968 operator _Op(const valarray<_Tp>& __v, const valarray<_Tp>& __w) \
970 _GLIBCXX_DEBUG_ASSERT(__v.size() == __w.size()); \
971 typedef _BinClos<_Name,_ValArray,_ValArray,_Tp,_Tp> _Closure; \
972 typedef typename __fun<_Name, _Tp>::result_type _Rt; \
973 return _Expr<_Closure, _Rt>(_Closure(__v, __w)); \
976 template<typename _Tp> \
977 inline _Expr<_BinClos<_Name,_ValArray,_Constant,_Tp,_Tp>, \
978 typename __fun<_Name, _Tp>::result_type> \
979 operator _Op(const valarray<_Tp>& __v, const _Tp& __t) \
981 typedef _BinClos<_Name,_ValArray,_Constant,_Tp,_Tp> _Closure; \
982 typedef typename __fun<_Name, _Tp>::result_type _Rt; \
983 return _Expr<_Closure, _Rt>(_Closure(__v, __t)); \
986 template<typename _Tp> \
987 inline _Expr<_BinClos<_Name,_Constant,_ValArray,_Tp,_Tp>, \
988 typename __fun<_Name, _Tp>::result_type> \
989 operator _Op(const _Tp& __t, const valarray<_Tp>& __v) \
991 typedef _BinClos<_Name,_Constant,_ValArray,_Tp,_Tp> _Closure; \
992 typedef typename __fun<_Name, _Tp>::result_type _Rt; \
993 return _Expr<_Closure, _Tp>(_Closure(__t, __v)); \
996 _DEFINE_BINARY_OPERATOR(+, __plus
)
997 _DEFINE_BINARY_OPERATOR(-, __minus
)
998 _DEFINE_BINARY_OPERATOR(*, __multiplies
)
999 _DEFINE_BINARY_OPERATOR(/, __divides
)
1000 _DEFINE_BINARY_OPERATOR(%, __modulus
)
1001 _DEFINE_BINARY_OPERATOR(^, __bitwise_xor
)
1002 _DEFINE_BINARY_OPERATOR(&, __bitwise_and
)
1003 _DEFINE_BINARY_OPERATOR(|, __bitwise_or
)
1004 _DEFINE_BINARY_OPERATOR(<<, __shift_left
)
1005 _DEFINE_BINARY_OPERATOR(>>, __shift_right
)
1006 _DEFINE_BINARY_OPERATOR(&&, __logical_and
)
1007 _DEFINE_BINARY_OPERATOR(||, __logical_or
)
1008 _DEFINE_BINARY_OPERATOR(==, __equal_to
)
1009 _DEFINE_BINARY_OPERATOR(!=, __not_equal_to
)
1010 _DEFINE_BINARY_OPERATOR(<, __less
)
1011 _DEFINE_BINARY_OPERATOR(>, __greater
)
1012 _DEFINE_BINARY_OPERATOR(<=, __less_equal
)
1013 _DEFINE_BINARY_OPERATOR(>=, __greater_equal
)
1017 #endif /* _GLIBCXX_VALARRAY */