1// The template and inlines for the -*- C++ -*- valarray class. 2 3// Copyright (C) 1997-2021 Free Software Foundation, Inc. 4// 5// This file is part of the GNU ISO C++ Library. This library is free 6// software; you can redistribute it and/or modify it under the 7// terms of the GNU General Public License as published by the 8// Free Software Foundation; either version 3, or (at your option) 9// any later version. 10 11// This library is distributed in the hope that it will be useful, 12// but WITHOUT ANY WARRANTY; without even the implied warranty of 13// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 14// GNU General Public License for more details. 15 16// Under Section 7 of GPL version 3, you are granted additional 17// permissions described in the GCC Runtime Library Exception, version 18// 3.1, as published by the Free Software Foundation. 19 20// You should have received a copy of the GNU General Public License and 21// a copy of the GCC Runtime Library Exception along with this program; 22// see the files COPYING3 and COPYING.RUNTIME respectively. If not, see 23// <http://www.gnu.org/licenses/>. 24 25/** @file include/valarray 26 * This is a Standard C++ Library header. 27 */ 28 29// Written by Gabriel Dos Reis <Gabriel.Dos-Reis@DPTMaths.ENS-Cachan.Fr> 30 31#ifndef _GLIBCXX_VALARRAY 32#define _GLIBCXX_VALARRAY 1 33 34#pragma GCC system_header 35 36#include <bits/c++config.h> 37#include <cmath> 38#include <algorithm> 39#include <debug/debug.h> 40#if __cplusplus >= 201103L 41#include <initializer_list> 42#endif 43 44namespace std _GLIBCXX_VISIBILITY(default) 45{ 46_GLIBCXX_BEGIN_NAMESPACE_VERSION 47 48 template<class _Clos, typename _Tp> 49 class _Expr; 50 51 template<typename _Tp1, typename _Tp2> 52 class _ValArray; 53 54namespace __detail 55{ 56 template<class _Oper, template<class, class> class _Meta, class _Dom> 57 struct _UnClos; 58 59 template<class _Oper, 60 template<class, class> class _Meta1, 61 template<class, class> class _Meta2, 62 class _Dom1, class _Dom2> 63 class _BinClos; 64 65 template<template<class, class> class _Meta, class _Dom> 66 class _SClos; 67 68 template<template<class, class> class _Meta, class _Dom> 69 class _GClos; 70 71 template<template<class, class> class _Meta, class _Dom> 72 class _IClos; 73 74 template<template<class, class> class _Meta, class _Dom> 75 class _ValFunClos; 76 77 template<template<class, class> class _Meta, class _Dom> 78 class _RefFunClos; 79} // namespace __detail 80 81 using __detail::_UnClos; 82 using __detail::_BinClos; 83 using __detail::_SClos; 84 using __detail::_GClos; 85 using __detail::_IClos; 86 using __detail::_ValFunClos; 87 using __detail::_RefFunClos; 88 89 template<class _Tp> class valarray; // An array of type _Tp 90 class slice; // BLAS-like slice out of an array 91 template<class _Tp> class slice_array; 92 class gslice; // generalized slice out of an array 93 template<class _Tp> class gslice_array; 94 template<class _Tp> class mask_array; // masked array 95 template<class _Tp> class indirect_array; // indirected array 96 97_GLIBCXX_END_NAMESPACE_VERSION 98} // namespace 99 100#include <bits/valarray_array.h> 101#include <bits/valarray_before.h> 102 103namespace std _GLIBCXX_VISIBILITY(default) 104{ 105_GLIBCXX_BEGIN_NAMESPACE_VERSION 106 107 /** 108 * @defgroup numeric_arrays Numeric Arrays 109 * @ingroup numerics 110 * 111 * Classes and functions for representing and manipulating arrays of elements. 112 * @{ 113 */ 114 115 /** 116 * @brief Smart array designed to support numeric processing. 117 * 118 * A valarray is an array that provides constraints intended to allow for 119 * effective optimization of numeric array processing by reducing the 120 * aliasing that can result from pointer representations. It represents a 121 * one-dimensional array from which different multidimensional subsets can 122 * be accessed and modified. 123 * 124 * @tparam _Tp Type of object in the array. 125 */ 126 template<class _Tp> 127 class valarray 128 { 129 template<class _Op> 130 struct _UnaryOp 131 { 132 typedef typename __fun<_Op, _Tp>::result_type __rt; 133 typedef _Expr<_UnClos<_Op, _ValArray, _Tp>, __rt> _Rt; 134 }; 135 public: 136 typedef _Tp value_type; 137 138 // _lib.valarray.cons_ construct/destroy: 139 /// Construct an empty array. 140 valarray(); 141 142 /// Construct an array with @a n elements. 143 explicit valarray(size_t); 144 145 /// Construct an array with @a n elements initialized to @a t. 146 valarray(const _Tp&, size_t); 147 148 /// Construct an array initialized to the first @a n elements of @a t. 149 valarray(const _Tp* __restrict__, size_t); 150 151 /// Copy constructor. 152 valarray(const valarray&); 153 154#if __cplusplus >= 201103L 155 /// Move constructor. 156 valarray(valarray&&) noexcept; 157#endif 158 159 /// Construct an array with the same size and values in @a sa. 160 valarray(const slice_array<_Tp>&); 161 162 /// Construct an array with the same size and values in @a ga. 163 valarray(const gslice_array<_Tp>&); 164 165 /// Construct an array with the same size and values in @a ma. 166 valarray(const mask_array<_Tp>&); 167 168 /// Construct an array with the same size and values in @a ia. 169 valarray(const indirect_array<_Tp>&); 170 171#if __cplusplus >= 201103L 172 /// Construct an array with an initializer_list of values. 173 valarray(initializer_list<_Tp>); 174#endif 175 176 template<class _Dom> 177 valarray(const _Expr<_Dom, _Tp>& __e); 178 179 ~valarray() _GLIBCXX_NOEXCEPT; 180 181 // _lib.valarray.assign_ assignment: 182 /** 183 * @brief Assign elements to an array. 184 * 185 * Assign elements of array to values in @a v. 186 * 187 * @param __v Valarray to get values from. 188 */ 189 valarray<_Tp>& operator=(const valarray<_Tp>& __v); 190 191#if __cplusplus >= 201103L 192 /** 193 * @brief Move assign elements to an array. 194 * 195 * Move assign elements of array to values in @a v. 196 * 197 * @param __v Valarray to get values from. 198 */ 199 valarray<_Tp>& operator=(valarray<_Tp>&& __v) noexcept; 200#endif 201 202 /** 203 * @brief Assign elements to a value. 204 * 205 * Assign all elements of array to @a t. 206 * 207 * @param __t Value for elements. 208 */ 209 valarray<_Tp>& operator=(const _Tp& __t); 210 211 /** 212 * @brief Assign elements to an array subset. 213 * 214 * Assign elements of array to values in @a sa. Results are undefined 215 * if @a sa does not have the same size as this array. 216 * 217 * @param __sa Array slice to get values from. 218 */ 219 valarray<_Tp>& operator=(const slice_array<_Tp>& __sa); 220 221 /** 222 * @brief Assign elements to an array subset. 223 * 224 * Assign elements of array to values in @a ga. Results are undefined 225 * if @a ga does not have the same size as this array. 226 * 227 * @param __ga Array slice to get values from. 228 */ 229 valarray<_Tp>& operator=(const gslice_array<_Tp>& __ga); 230 231 /** 232 * @brief Assign elements to an array subset. 233 * 234 * Assign elements of array to values in @a ma. Results are undefined 235 * if @a ma does not have the same size as this array. 236 * 237 * @param __ma Array slice to get values from. 238 */ 239 valarray<_Tp>& operator=(const mask_array<_Tp>& __ma); 240 241 /** 242 * @brief Assign elements to an array subset. 243 * 244 * Assign elements of array to values in @a ia. Results are undefined 245 * if @a ia does not have the same size as this array. 246 * 247 * @param __ia Array slice to get values from. 248 */ 249 valarray<_Tp>& operator=(const indirect_array<_Tp>& __ia); 250 251#if __cplusplus >= 201103L 252 /** 253 * @brief Assign elements to an initializer_list. 254 * 255 * Assign elements of array to values in @a __l. Results are undefined 256 * if @a __l does not have the same size as this array. 257 * 258 * @param __l initializer_list to get values from. 259 */ 260 valarray& operator=(initializer_list<_Tp> __l); 261#endif 262 263 template<class _Dom> valarray<_Tp>& 264 operator= (const _Expr<_Dom, _Tp>&); 265 266 // _lib.valarray.access_ element access: 267 /** 268 * Return a reference to the i'th array element. 269 * 270 * @param __i Index of element to return. 271 * @return Reference to the i'th element. 272 */ 273 _Tp& operator[](size_t __i); 274 275 // _GLIBCXX_RESOLVE_LIB_DEFECTS 276 // 389. Const overload of valarray::operator[] returns by value. 277 const _Tp& operator[](size_t) const; 278 279 // _lib.valarray.sub_ subset operations: 280 /** 281 * @brief Return an array subset. 282 * 283 * Returns a new valarray containing the elements of the array 284 * indicated by the slice argument. The new valarray has the same size 285 * as the input slice. @see slice. 286 * 287 * @param __s The source slice. 288 * @return New valarray containing elements in @a __s. 289 */ 290 _Expr<_SClos<_ValArray, _Tp>, _Tp> operator[](slice __s) const; 291 292 /** 293 * @brief Return a reference to an array subset. 294 * 295 * Returns a new valarray containing the elements of the array 296 * indicated by the slice argument. The new valarray has the same size 297 * as the input slice. @see slice. 298 * 299 * @param __s The source slice. 300 * @return New valarray containing elements in @a __s. 301 */ 302 slice_array<_Tp> operator[](slice __s); 303 304 /** 305 * @brief Return an array subset. 306 * 307 * Returns a slice_array referencing the elements of the array 308 * indicated by the slice argument. @see gslice. 309 * 310 * @param __s The source slice. 311 * @return Slice_array referencing elements indicated by @a __s. 312 */ 313 _Expr<_GClos<_ValArray, _Tp>, _Tp> operator[](const gslice& __s) const; 314 315 /** 316 * @brief Return a reference to an array subset. 317 * 318 * Returns a new valarray containing the elements of the array 319 * indicated by the gslice argument. The new valarray has 320 * the same size as the input gslice. @see gslice. 321 * 322 * @param __s The source gslice. 323 * @return New valarray containing elements in @a __s. 324 */ 325 gslice_array<_Tp> operator[](const gslice& __s); 326 327 /** 328 * @brief Return an array subset. 329 * 330 * Returns a new valarray containing the elements of the array 331 * indicated by the argument. The input is a valarray of bool which 332 * represents a bitmask indicating which elements should be copied into 333 * the new valarray. Each element of the array is added to the return 334 * valarray if the corresponding element of the argument is true. 335 * 336 * @param __m The valarray bitmask. 337 * @return New valarray containing elements indicated by @a __m. 338 */ 339 valarray<_Tp> operator[](const valarray<bool>& __m) const; 340 341 /** 342 * @brief Return a reference to an array subset. 343 * 344 * Returns a new mask_array referencing the elements of the array 345 * indicated by the argument. The input is a valarray of bool which 346 * represents a bitmask indicating which elements are part of the 347 * subset. Elements of the array are part of the subset if the 348 * corresponding element of the argument is true. 349 * 350 * @param __m The valarray bitmask. 351 * @return New valarray containing elements indicated by @a __m. 352 */ 353 mask_array<_Tp> operator[](const valarray<bool>& __m); 354 355 /** 356 * @brief Return an array subset. 357 * 358 * Returns a new valarray containing the elements of the array 359 * indicated by the argument. The elements in the argument are 360 * interpreted as the indices of elements of this valarray to copy to 361 * the return valarray. 362 * 363 * @param __i The valarray element index list. 364 * @return New valarray containing elements in @a __s. 365 */ 366 _Expr<_IClos<_ValArray, _Tp>, _Tp> 367 operator[](const valarray<size_t>& __i) const; 368 369 /** 370 * @brief Return a reference to an array subset. 371 * 372 * Returns an indirect_array referencing the elements of the array 373 * indicated by the argument. The elements in the argument are 374 * interpreted as the indices of elements of this valarray to include 375 * in the subset. The returned indirect_array refers to these 376 * elements. 377 * 378 * @param __i The valarray element index list. 379 * @return Indirect_array referencing elements in @a __i. 380 */ 381 indirect_array<_Tp> operator[](const valarray<size_t>& __i); 382 383 // _lib.valarray.unary_ unary operators: 384 /// Return a new valarray by applying unary + to each element. 385 typename _UnaryOp<__unary_plus>::_Rt operator+() const; 386 387 /// Return a new valarray by applying unary - to each element. 388 typename _UnaryOp<__negate>::_Rt operator-() const; 389 390 /// Return a new valarray by applying unary ~ to each element. 391 typename _UnaryOp<__bitwise_not>::_Rt operator~() const; 392 393 /// Return a new valarray by applying unary ! to each element. 394 typename _UnaryOp<__logical_not>::_Rt operator!() const; 395 396 // _lib.valarray.cassign_ computed assignment: 397 /// Multiply each element of array by @a t. 398 valarray<_Tp>& operator*=(const _Tp&); 399 400 /// Divide each element of array by @a t. 401 valarray<_Tp>& operator/=(const _Tp&); 402 403 /// Set each element e of array to e % @a t. 404 valarray<_Tp>& operator%=(const _Tp&); 405 406 /// Add @a t to each element of array. 407 valarray<_Tp>& operator+=(const _Tp&); 408 409 /// Subtract @a t to each element of array. 410 valarray<_Tp>& operator-=(const _Tp&); 411 412 /// Set each element e of array to e ^ @a t. 413 valarray<_Tp>& operator^=(const _Tp&); 414 415 /// Set each element e of array to e & @a t. 416 valarray<_Tp>& operator&=(const _Tp&); 417 418 /// Set each element e of array to e | @a t. 419 valarray<_Tp>& operator|=(const _Tp&); 420 421 /// Left shift each element e of array by @a t bits. 422 valarray<_Tp>& operator<<=(const _Tp&); 423 424 /// Right shift each element e of array by @a t bits. 425 valarray<_Tp>& operator>>=(const _Tp&); 426 427 /// Multiply elements of array by corresponding elements of @a v. 428 valarray<_Tp>& operator*=(const valarray<_Tp>&); 429 430 /// Divide elements of array by corresponding elements of @a v. 431 valarray<_Tp>& operator/=(const valarray<_Tp>&); 432 433 /// Modulo elements of array by corresponding elements of @a v. 434 valarray<_Tp>& operator%=(const valarray<_Tp>&); 435 436 /// Add corresponding elements of @a v to elements of array. 437 valarray<_Tp>& operator+=(const valarray<_Tp>&); 438 439 /// Subtract corresponding elements of @a v from elements of array. 440 valarray<_Tp>& operator-=(const valarray<_Tp>&); 441 442 /// Logical xor corresponding elements of @a v with elements of array. 443 valarray<_Tp>& operator^=(const valarray<_Tp>&); 444 445 /// Logical or corresponding elements of @a v with elements of array. 446 valarray<_Tp>& operator|=(const valarray<_Tp>&); 447 448 /// Logical and corresponding elements of @a v with elements of array. 449 valarray<_Tp>& operator&=(const valarray<_Tp>&); 450 451 /// Left shift elements of array by corresponding elements of @a v. 452 valarray<_Tp>& operator<<=(const valarray<_Tp>&); 453 454 /// Right shift elements of array by corresponding elements of @a v. 455 valarray<_Tp>& operator>>=(const valarray<_Tp>&); 456 457 template<class _Dom> 458 valarray<_Tp>& operator*=(const _Expr<_Dom, _Tp>&); 459 template<class _Dom> 460 valarray<_Tp>& operator/=(const _Expr<_Dom, _Tp>&); 461 template<class _Dom> 462 valarray<_Tp>& operator%=(const _Expr<_Dom, _Tp>&); 463 template<class _Dom> 464 valarray<_Tp>& operator+=(const _Expr<_Dom, _Tp>&); 465 template<class _Dom> 466 valarray<_Tp>& operator-=(const _Expr<_Dom, _Tp>&); 467 template<class _Dom> 468 valarray<_Tp>& operator^=(const _Expr<_Dom, _Tp>&); 469 template<class _Dom> 470 valarray<_Tp>& operator|=(const _Expr<_Dom, _Tp>&); 471 template<class _Dom> 472 valarray<_Tp>& operator&=(const _Expr<_Dom, _Tp>&); 473 template<class _Dom> 474 valarray<_Tp>& operator<<=(const _Expr<_Dom, _Tp>&); 475 template<class _Dom> 476 valarray<_Tp>& operator>>=(const _Expr<_Dom, _Tp>&); 477 478 // _lib.valarray.members_ member functions: 479#if __cplusplus >= 201103L 480 /// Swap. 481 void swap(valarray<_Tp>& __v) noexcept; 482#endif 483 484 /// Return the number of elements in array. 485 size_t size() const; 486 487 /** 488 * @brief Return the sum of all elements in the array. 489 * 490 * Accumulates the sum of all elements into a Tp using +=. The order 491 * of adding the elements is unspecified. 492 */ 493 _Tp sum() const; 494 495 /// Return the minimum element using operator<(). 496 _Tp min() const; 497 498 /// Return the maximum element using operator<(). 499 _Tp max() const; 500 501 /** 502 * @brief Return a shifted array. 503 * 504 * A new valarray is constructed as a copy of this array with elements 505 * in shifted positions. For an element with index i, the new position 506 * is i - n. The new valarray has the same size as the current one. 507 * New elements without a value are set to 0. Elements whose new 508 * position is outside the bounds of the array are discarded. 509 * 510 * Positive arguments shift toward index 0, discarding elements [0, n). 511 * Negative arguments discard elements from the top of the array. 512 * 513 * @param __n Number of element positions to shift. 514 * @return New valarray with elements in shifted positions. 515 */ 516 valarray<_Tp> shift (int __n) const; 517 518 /** 519 * @brief Return a rotated array. 520 * 521 * A new valarray is constructed as a copy of this array with elements 522 * in shifted positions. For an element with index i, the new position 523 * is (i - n) % size(). The new valarray has the same size as the 524 * current one. Elements that are shifted beyond the array bounds are 525 * shifted into the other end of the array. No elements are lost. 526 * 527 * Positive arguments shift toward index 0, wrapping around the top. 528 * Negative arguments shift towards the top, wrapping around to 0. 529 * 530 * @param __n Number of element positions to rotate. 531 * @return New valarray with elements in shifted positions. 532 */ 533 valarray<_Tp> cshift(int __n) const; 534 535 /** 536 * @brief Apply a function to the array. 537 * 538 * Returns a new valarray with elements assigned to the result of 539 * applying func to the corresponding element of this array. The new 540 * array has the same size as this one. 541 * 542 * @param func Function of Tp returning Tp to apply. 543 * @return New valarray with transformed elements. 544 */ 545 _Expr<_ValFunClos<_ValArray, _Tp>, _Tp> apply(_Tp func(_Tp)) const; 546 547 /** 548 * @brief Apply a function to the array. 549 * 550 * Returns a new valarray with elements assigned to the result of 551 * applying func to the corresponding element of this array. The new 552 * array has the same size as this one. 553 * 554 * @param func Function of const Tp& returning Tp to apply. 555 * @return New valarray with transformed elements. 556 */ 557 _Expr<_RefFunClos<_ValArray, _Tp>, _Tp> apply(_Tp func(const _Tp&)) const; 558 559 /** 560 * @brief Resize array. 561 * 562 * Resize this array to @a size and set all elements to @a c. All 563 * references and iterators are invalidated. 564 * 565 * @param __size New array size. 566 * @param __c New value for all elements. 567 */ 568 void resize(size_t __size, _Tp __c = _Tp()); 569 570 private: 571 size_t _M_size; 572 _Tp* __restrict__ _M_data; 573 574 friend struct _Array<_Tp>; 575 }; 576 577#if __cpp_deduction_guides >= 201606 578 template<typename _Tp, size_t _Nm> 579 valarray(const _Tp(&)[_Nm], size_t) -> valarray<_Tp>; 580#endif 581 582 template<typename _Tp> 583 inline const _Tp& 584 valarray<_Tp>::operator[](size_t __i) const 585 { 586 __glibcxx_requires_subscript(__i); 587 return _M_data[__i]; 588 } 589 590 template<typename _Tp> 591 inline _Tp& 592 valarray<_Tp>::operator[](size_t __i) 593 { 594 __glibcxx_requires_subscript(__i); 595 return _M_data[__i]; 596 } 597 598 /// @} group numeric_arrays 599 600_GLIBCXX_END_NAMESPACE_VERSION 601} // namespace 602 603#include <bits/valarray_after.h> 604#include <bits/slice_array.h> 605#include <bits/gslice.h> 606#include <bits/gslice_array.h> 607#include <bits/mask_array.h> 608#include <bits/indirect_array.h> 609 610namespace std _GLIBCXX_VISIBILITY(default) 611{ 612_GLIBCXX_BEGIN_NAMESPACE_VERSION 613 614 /** 615 * @addtogroup numeric_arrays 616 * @{ 617 */ 618 619 template<typename _Tp> 620 inline 621 valarray<_Tp>::valarray() : _M_size(0), _M_data(0) {} 622 623 template<typename _Tp> 624 inline 625 valarray<_Tp>::valarray(size_t __n) 626 : _M_size(__n), _M_data(__valarray_get_storage<_Tp>(__n)) 627 { std::__valarray_default_construct(_M_data, _M_data + __n); } 628 629 template<typename _Tp> 630 inline 631 valarray<_Tp>::valarray(const _Tp& __t, size_t __n) 632 : _M_size(__n), _M_data(__valarray_get_storage<_Tp>(__n)) 633 { std::__valarray_fill_construct(_M_data, _M_data + __n, __t); } 634 635 template<typename _Tp> 636 inline 637 valarray<_Tp>::valarray(const _Tp* __restrict__ __p, size_t __n) 638 : _M_size(__n), _M_data(__valarray_get_storage<_Tp>(__n)) 639 { 640 __glibcxx_assert(__p != 0 || __n == 0); 641 std::__valarray_copy_construct(__p, __p + __n, _M_data); 642 } 643 644 template<typename _Tp> 645 inline 646 valarray<_Tp>::valarray(const valarray<_Tp>& __v) 647 : _M_size(__v._M_size), _M_data(__valarray_get_storage<_Tp>(__v._M_size)) 648 { std::__valarray_copy_construct(__v._M_data, __v._M_data + _M_size, 649 _M_data); } 650 651#if __cplusplus >= 201103L 652 template<typename _Tp> 653 inline 654 valarray<_Tp>::valarray(valarray<_Tp>&& __v) noexcept 655 : _M_size(__v._M_size), _M_data(__v._M_data) 656 { 657 __v._M_size = 0; 658 __v._M_data = 0; 659 } 660#endif 661 662 template<typename _Tp> 663 inline 664 valarray<_Tp>::valarray(const slice_array<_Tp>& __sa) 665 : _M_size(__sa._M_sz), _M_data(__valarray_get_storage<_Tp>(__sa._M_sz)) 666 { 667 std::__valarray_copy_construct 668 (__sa._M_array, __sa._M_sz, __sa._M_stride, _Array<_Tp>(_M_data)); 669 } 670 671 template<typename _Tp> 672 inline 673 valarray<_Tp>::valarray(const gslice_array<_Tp>& __ga) 674 : _M_size(__ga._M_index.size()), 675 _M_data(__valarray_get_storage<_Tp>(_M_size)) 676 { 677 std::__valarray_copy_construct 678 (__ga._M_array, _Array<size_t>(__ga._M_index), 679 _Array<_Tp>(_M_data), _M_size); 680 } 681 682 template<typename _Tp> 683 inline 684 valarray<_Tp>::valarray(const mask_array<_Tp>& __ma) 685 : _M_size(__ma._M_sz), _M_data(__valarray_get_storage<_Tp>(__ma._M_sz)) 686 { 687 std::__valarray_copy_construct 688 (__ma._M_array, __ma._M_mask, _Array<_Tp>(_M_data), _M_size); 689 } 690 691 template<typename _Tp> 692 inline 693 valarray<_Tp>::valarray(const indirect_array<_Tp>& __ia) 694 : _M_size(__ia._M_sz), _M_data(__valarray_get_storage<_Tp>(__ia._M_sz)) 695 { 696 std::__valarray_copy_construct 697 (__ia._M_array, __ia._M_index, _Array<_Tp>(_M_data), _M_size); 698 } 699 700#if __cplusplus >= 201103L 701 template<typename _Tp> 702 inline 703 valarray<_Tp>::valarray(initializer_list<_Tp> __l) 704 : _M_size(__l.size()), _M_data(__valarray_get_storage<_Tp>(__l.size())) 705 { std::__valarray_copy_construct(__l.begin(), __l.end(), _M_data); } 706#endif 707 708 template<typename _Tp> template<class _Dom> 709 inline 710 valarray<_Tp>::valarray(const _Expr<_Dom, _Tp>& __e) 711 : _M_size(__e.size()), _M_data(__valarray_get_storage<_Tp>(_M_size)) 712 { std::__valarray_copy_construct(__e, _M_size, _Array<_Tp>(_M_data)); } 713 714 template<typename _Tp> 715 inline 716 valarray<_Tp>::~valarray() _GLIBCXX_NOEXCEPT 717 { 718 std::__valarray_destroy_elements(_M_data, _M_data + _M_size); 719 std::__valarray_release_memory(_M_data); 720 } 721 722 template<typename _Tp> 723 inline valarray<_Tp>& 724 valarray<_Tp>::operator=(const valarray<_Tp>& __v) 725 { 726 // _GLIBCXX_RESOLVE_LIB_DEFECTS 727 // 630. arrays of valarray. 728 if (_M_size == __v._M_size) 729 std::__valarray_copy(__v._M_data, _M_size, _M_data); 730 else 731 { 732 if (_M_data) 733 { 734 std::__valarray_destroy_elements(_M_data, _M_data + _M_size); 735 std::__valarray_release_memory(_M_data); 736 } 737 _M_size = __v._M_size; 738 _M_data = __valarray_get_storage<_Tp>(_M_size); 739 std::__valarray_copy_construct(__v._M_data, __v._M_data + _M_size, 740 _M_data); 741 } 742 return *this; 743 } 744 745#if __cplusplus >= 201103L 746 template<typename _Tp> 747 inline valarray<_Tp>& 748 valarray<_Tp>::operator=(valarray<_Tp>&& __v) noexcept 749 { 750 if (_M_data) 751 { 752 std::__valarray_destroy_elements(_M_data, _M_data + _M_size); 753 std::__valarray_release_memory(_M_data); 754 } 755 _M_size = __v._M_size; 756 _M_data = __v._M_data; 757 __v._M_size = 0; 758 __v._M_data = 0; 759 return *this; 760 } 761 762 template<typename _Tp> 763 inline valarray<_Tp>& 764 valarray<_Tp>::operator=(initializer_list<_Tp> __l) 765 { 766 // _GLIBCXX_RESOLVE_LIB_DEFECTS 767 // 630. arrays of valarray. 768 if (_M_size == __l.size()) 769 std::__valarray_copy(__l.begin(), __l.size(), _M_data); 770 else 771 { 772 if (_M_data) 773 { 774 std::__valarray_destroy_elements(_M_data, _M_data + _M_size); 775 std::__valarray_release_memory(_M_data); 776 } 777 _M_size = __l.size(); 778 _M_data = __valarray_get_storage<_Tp>(_M_size); 779 std::__valarray_copy_construct(__l.begin(), __l.begin() + _M_size, 780 _M_data); 781 } 782 return *this; 783 } 784#endif 785 786 template<typename _Tp> 787 inline valarray<_Tp>& 788 valarray<_Tp>::operator=(const _Tp& __t) 789 { 790 std::__valarray_fill(_M_data, _M_size, __t); 791 return *this; 792 } 793 794 template<typename _Tp> 795 inline valarray<_Tp>& 796 valarray<_Tp>::operator=(const slice_array<_Tp>& __sa) 797 { 798 __glibcxx_assert(_M_size == __sa._M_sz); 799 std::__valarray_copy(__sa._M_array, __sa._M_sz, 800 __sa._M_stride, _Array<_Tp>(_M_data)); 801 return *this; 802 } 803 804 template<typename _Tp> 805 inline valarray<_Tp>& 806 valarray<_Tp>::operator=(const gslice_array<_Tp>& __ga) 807 { 808 __glibcxx_assert(_M_size == __ga._M_index.size()); 809 std::__valarray_copy(__ga._M_array, _Array<size_t>(__ga._M_index), 810 _Array<_Tp>(_M_data), _M_size); 811 return *this; 812 } 813 814 template<typename _Tp> 815 inline valarray<_Tp>& 816 valarray<_Tp>::operator=(const mask_array<_Tp>& __ma) 817 { 818 __glibcxx_assert(_M_size == __ma._M_sz); 819 std::__valarray_copy(__ma._M_array, __ma._M_mask, 820 _Array<_Tp>(_M_data), _M_size); 821 return *this; 822 } 823 824 template<typename _Tp> 825 inline valarray<_Tp>& 826 valarray<_Tp>::operator=(const indirect_array<_Tp>& __ia) 827 { 828 __glibcxx_assert(_M_size == __ia._M_sz); 829 std::__valarray_copy(__ia._M_array, __ia._M_index, 830 _Array<_Tp>(_M_data), _M_size); 831 return *this; 832 } 833 834 template<typename _Tp> template<class _Dom> 835 inline valarray<_Tp>& 836 valarray<_Tp>::operator=(const _Expr<_Dom, _Tp>& __e) 837 { 838 // _GLIBCXX_RESOLVE_LIB_DEFECTS 839 // 630. arrays of valarray. 840 if (_M_size == __e.size()) 841 std::__valarray_copy(__e, _M_size, _Array<_Tp>(_M_data)); 842 else 843 { 844 if (_M_data) 845 { 846 std::__valarray_destroy_elements(_M_data, _M_data + _M_size); 847 std::__valarray_release_memory(_M_data); 848 } 849 _M_size = __e.size(); 850 _M_data = __valarray_get_storage<_Tp>(_M_size); 851 std::__valarray_copy_construct(__e, _M_size, _Array<_Tp>(_M_data)); 852 } 853 return *this; 854 } 855 856 template<typename _Tp> 857 inline _Expr<_SClos<_ValArray,_Tp>, _Tp> 858 valarray<_Tp>::operator[](slice __s) const 859 { 860 typedef _SClos<_ValArray,_Tp> _Closure; 861 return _Expr<_Closure, _Tp>(_Closure (_Array<_Tp>(_M_data), __s)); 862 } 863 864 template<typename _Tp> 865 inline slice_array<_Tp> 866 valarray<_Tp>::operator[](slice __s) 867 { return slice_array<_Tp>(_Array<_Tp>(_M_data), __s); } 868 869 template<typename _Tp> 870 inline _Expr<_GClos<_ValArray,_Tp>, _Tp> 871 valarray<_Tp>::operator[](const gslice& __gs) const 872 { 873 typedef _GClos<_ValArray,_Tp> _Closure; 874 return _Expr<_Closure, _Tp> 875 (_Closure(_Array<_Tp>(_M_data), __gs._M_index->_M_index)); 876 } 877 878 template<typename _Tp> 879 inline gslice_array<_Tp> 880 valarray<_Tp>::operator[](const gslice& __gs) 881 { 882 return gslice_array<_Tp> 883 (_Array<_Tp>(_M_data), __gs._M_index->_M_index); 884 } 885 886 template<typename _Tp> 887 inline valarray<_Tp> 888 valarray<_Tp>::operator[](const valarray<bool>& __m) const 889 { 890 size_t __s = 0; 891 size_t __e = __m.size(); 892 for (size_t __i=0; __i<__e; ++__i) 893 if (__m[__i]) ++__s; 894 return valarray<_Tp>(mask_array<_Tp>(_Array<_Tp>(_M_data), __s, 895 _Array<bool> (__m))); 896 } 897 898 template<typename _Tp> 899 inline mask_array<_Tp> 900 valarray<_Tp>::operator[](const valarray<bool>& __m) 901 { 902 size_t __s = 0; 903 size_t __e = __m.size(); 904 for (size_t __i=0; __i<__e; ++__i) 905 if (__m[__i]) ++__s; 906 return mask_array<_Tp>(_Array<_Tp>(_M_data), __s, _Array<bool>(__m)); 907 } 908 909 template<typename _Tp> 910 inline _Expr<_IClos<_ValArray,_Tp>, _Tp> 911 valarray<_Tp>::operator[](const valarray<size_t>& __i) const 912 { 913 typedef _IClos<_ValArray,_Tp> _Closure; 914 return _Expr<_Closure, _Tp>(_Closure(*this, __i)); 915 } 916 917 template<typename _Tp> 918 inline indirect_array<_Tp> 919 valarray<_Tp>::operator[](const valarray<size_t>& __i) 920 { 921 return indirect_array<_Tp>(_Array<_Tp>(_M_data), __i.size(), 922 _Array<size_t>(__i)); 923 } 924 925#if __cplusplus >= 201103L 926 template<class _Tp> 927 inline void 928 valarray<_Tp>::swap(valarray<_Tp>& __v) noexcept 929 { 930 std::swap(_M_size, __v._M_size); 931 std::swap(_M_data, __v._M_data); 932 } 933#endif 934 935 template<class _Tp> 936 inline size_t 937 valarray<_Tp>::size() const 938 { return _M_size; } 939 940 template<class _Tp> 941 inline _Tp 942 valarray<_Tp>::sum() const 943 { 944 __glibcxx_assert(_M_size > 0); 945 return std::__valarray_sum(_M_data, _M_data + _M_size); 946 } 947 948 template<class _Tp> 949 inline valarray<_Tp> 950 valarray<_Tp>::shift(int __n) const 951 { 952 valarray<_Tp> __ret; 953 954 if (_M_size == 0) 955 return __ret; 956 957 _Tp* __restrict__ __tmp_M_data = 958 std::__valarray_get_storage<_Tp>(_M_size); 959 960 if (__n == 0) 961 std::__valarray_copy_construct(_M_data, 962 _M_data + _M_size, __tmp_M_data); 963 else if (__n > 0) // shift left 964 { 965 if (size_t(__n) > _M_size) 966 __n = int(_M_size); 967 968 std::__valarray_copy_construct(_M_data + __n, 969 _M_data + _M_size, __tmp_M_data); 970 std::__valarray_default_construct(__tmp_M_data + _M_size - __n, 971 __tmp_M_data + _M_size); 972 } 973 else // shift right 974 { 975 if (-size_t(__n) > _M_size) 976 __n = -int(_M_size); 977 978 std::__valarray_copy_construct(_M_data, _M_data + _M_size + __n, 979 __tmp_M_data - __n); 980 std::__valarray_default_construct(__tmp_M_data, 981 __tmp_M_data - __n); 982 } 983 984 __ret._M_size = _M_size; 985 __ret._M_data = __tmp_M_data; 986 return __ret; 987 } 988 989 template<class _Tp> 990 inline valarray<_Tp> 991 valarray<_Tp>::cshift(int __n) const 992 { 993 valarray<_Tp> __ret; 994 995 if (_M_size == 0) 996 return __ret; 997 998 _Tp* __restrict__ __tmp_M_data = 999 std::__valarray_get_storage<_Tp>(_M_size); 1000 1001 if (__n == 0) 1002 std::__valarray_copy_construct(_M_data, 1003 _M_data + _M_size, __tmp_M_data); 1004 else if (__n > 0) // cshift left 1005 { 1006 if (size_t(__n) > _M_size) 1007 __n = int(__n % _M_size); 1008 1009 std::__valarray_copy_construct(_M_data, _M_data + __n, 1010 __tmp_M_data + _M_size - __n); 1011 std::__valarray_copy_construct(_M_data + __n, _M_data + _M_size, 1012 __tmp_M_data); 1013 } 1014 else // cshift right 1015 { 1016 if (-size_t(__n) > _M_size) 1017 __n = -int(-size_t(__n) % _M_size); 1018 1019 std::__valarray_copy_construct(_M_data + _M_size + __n, 1020 _M_data + _M_size, __tmp_M_data); 1021 std::__valarray_copy_construct(_M_data, _M_data + _M_size + __n, 1022 __tmp_M_data - __n); 1023 } 1024 1025 __ret._M_size = _M_size; 1026 __ret._M_data = __tmp_M_data; 1027 return __ret; 1028 } 1029 1030 template<class _Tp> 1031 inline void 1032 valarray<_Tp>::resize(size_t __n, _Tp __c) 1033 { 1034 // This complication is so to make valarray<valarray<T> > work 1035 // even though it is not required by the standard. Nobody should 1036 // be saying valarray<valarray<T> > anyway. See the specs. 1037 std::__valarray_destroy_elements(_M_data, _M_data + _M_size); 1038 if (_M_size != __n) 1039 { 1040 std::__valarray_release_memory(_M_data); 1041 _M_size = __n; 1042 _M_data = __valarray_get_storage<_Tp>(__n); 1043 } 1044 std::__valarray_fill_construct(_M_data, _M_data + __n, __c); 1045 } 1046 1047 template<typename _Tp> 1048 inline _Tp 1049 valarray<_Tp>::min() const 1050 { 1051 __glibcxx_assert(_M_size > 0); 1052 return *std::min_element(_M_data, _M_data + _M_size); 1053 } 1054 1055 template<typename _Tp> 1056 inline _Tp 1057 valarray<_Tp>::max() const 1058 { 1059 __glibcxx_assert(_M_size > 0); 1060 return *std::max_element(_M_data, _M_data + _M_size); 1061 } 1062 1063 template<class _Tp> 1064 inline _Expr<_ValFunClos<_ValArray, _Tp>, _Tp> 1065 valarray<_Tp>::apply(_Tp func(_Tp)) const 1066 { 1067 typedef _ValFunClos<_ValArray, _Tp> _Closure; 1068 return _Expr<_Closure, _Tp>(_Closure(*this, func)); 1069 } 1070 1071 template<class _Tp> 1072 inline _Expr<_RefFunClos<_ValArray, _Tp>, _Tp> 1073 valarray<_Tp>::apply(_Tp func(const _Tp &)) const 1074 { 1075 typedef _RefFunClos<_ValArray, _Tp> _Closure; 1076 return _Expr<_Closure, _Tp>(_Closure(*this, func)); 1077 } 1078 1079#define _DEFINE_VALARRAY_UNARY_OPERATOR(_Op, _Name) \ 1080 template<typename _Tp> \ 1081 inline typename valarray<_Tp>::template _UnaryOp<_Name>::_Rt \ 1082 valarray<_Tp>::operator _Op() const \ 1083 { \ 1084 typedef _UnClos<_Name, _ValArray, _Tp> _Closure; \ 1085 typedef typename __fun<_Name, _Tp>::result_type _Rt; \ 1086 return _Expr<_Closure, _Rt>(_Closure(*this)); \ 1087 } 1088 1089 _DEFINE_VALARRAY_UNARY_OPERATOR(+, __unary_plus) 1090 _DEFINE_VALARRAY_UNARY_OPERATOR(-, __negate) 1091 _DEFINE_VALARRAY_UNARY_OPERATOR(~, __bitwise_not) 1092 _DEFINE_VALARRAY_UNARY_OPERATOR (!, __logical_not) 1093 1094#undef _DEFINE_VALARRAY_UNARY_OPERATOR 1095 1096#define _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(_Op, _Name) \ 1097 template<class _Tp> \ 1098 inline valarray<_Tp>& \ 1099 valarray<_Tp>::operator _Op##=(const _Tp &__t) \ 1100 { \ 1101 _Array_augmented_##_Name(_Array<_Tp>(_M_data), _M_size, __t); \ 1102 return *this; \ 1103 } \ 1104 \ 1105 template<class _Tp> \ 1106 inline valarray<_Tp>& \ 1107 valarray<_Tp>::operator _Op##=(const valarray<_Tp> &__v) \ 1108 { \ 1109 __glibcxx_assert(_M_size == __v._M_size); \ 1110 _Array_augmented_##_Name(_Array<_Tp>(_M_data), _M_size, \ 1111 _Array<_Tp>(__v._M_data)); \ 1112 return *this; \ 1113 } 1114 1115_DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(+, __plus) 1116_DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(-, __minus) 1117_DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(*, __multiplies) 1118_DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(/, __divides) 1119_DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(%, __modulus) 1120_DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(^, __bitwise_xor) 1121_DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(&, __bitwise_and) 1122_DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(|, __bitwise_or) 1123_DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(<<, __shift_left) 1124_DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT(>>, __shift_right) 1125 1126#undef _DEFINE_VALARRAY_AUGMENTED_ASSIGNMENT 1127 1128#define _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(_Op, _Name) \ 1129 template<class _Tp> template<class _Dom> \ 1130 inline valarray<_Tp>& \ 1131 valarray<_Tp>::operator _Op##=(const _Expr<_Dom, _Tp>& __e) \ 1132 { \ 1133 _Array_augmented_##_Name(_Array<_Tp>(_M_data), __e, _M_size); \ 1134 return *this; \ 1135 } 1136 1137_DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(+, __plus) 1138_DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(-, __minus) 1139_DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(*, __multiplies) 1140_DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(/, __divides) 1141_DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(%, __modulus) 1142_DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(^, __bitwise_xor) 1143_DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(&, __bitwise_and) 1144_DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(|, __bitwise_or) 1145_DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(<<, __shift_left) 1146_DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT(>>, __shift_right) 1147 1148#undef _DEFINE_VALARRAY_EXPR_AUGMENTED_ASSIGNMENT 1149 1150 1151#define _DEFINE_BINARY_OPERATOR(_Op, _Name) \ 1152 template<typename _Tp> \ 1153 inline _Expr<_BinClos<_Name, _ValArray, _ValArray, _Tp, _Tp>, \ 1154 typename __fun<_Name, _Tp>::result_type> \ 1155 operator _Op(const valarray<_Tp>& __v, const valarray<_Tp>& __w) \ 1156 { \ 1157 __glibcxx_assert(__v.size() == __w.size()); \ 1158 typedef _BinClos<_Name, _ValArray, _ValArray, _Tp, _Tp> _Closure; \ 1159 typedef typename __fun<_Name, _Tp>::result_type _Rt; \ 1160 return _Expr<_Closure, _Rt>(_Closure(__v, __w)); \ 1161 } \ 1162 \ 1163 template<typename _Tp> \ 1164 inline _Expr<_BinClos<_Name, _ValArray,_Constant, _Tp, _Tp>, \ 1165 typename __fun<_Name, _Tp>::result_type> \ 1166 operator _Op(const valarray<_Tp>& __v, \ 1167 const typename valarray<_Tp>::value_type& __t) \ 1168 { \ 1169 typedef _BinClos<_Name, _ValArray, _Constant, _Tp, _Tp> _Closure; \ 1170 typedef typename __fun<_Name, _Tp>::result_type _Rt; \ 1171 return _Expr<_Closure, _Rt>(_Closure(__v, __t)); \ 1172 } \ 1173 \ 1174 template<typename _Tp> \ 1175 inline _Expr<_BinClos<_Name, _Constant, _ValArray, _Tp, _Tp>, \ 1176 typename __fun<_Name, _Tp>::result_type> \ 1177 operator _Op(const typename valarray<_Tp>::value_type& __t, \ 1178 const valarray<_Tp>& __v) \ 1179 { \ 1180 typedef _BinClos<_Name, _Constant, _ValArray, _Tp, _Tp> _Closure; \ 1181 typedef typename __fun<_Name, _Tp>::result_type _Rt; \ 1182 return _Expr<_Closure, _Rt>(_Closure(__t, __v)); \ 1183 } 1184 1185_DEFINE_BINARY_OPERATOR(+, __plus) 1186_DEFINE_BINARY_OPERATOR(-, __minus) 1187_DEFINE_BINARY_OPERATOR(*, __multiplies) 1188_DEFINE_BINARY_OPERATOR(/, __divides) 1189_DEFINE_BINARY_OPERATOR(%, __modulus) 1190_DEFINE_BINARY_OPERATOR(^, __bitwise_xor) 1191_DEFINE_BINARY_OPERATOR(&, __bitwise_and) 1192_DEFINE_BINARY_OPERATOR(|, __bitwise_or) 1193_DEFINE_BINARY_OPERATOR(<<, __shift_left) 1194_DEFINE_BINARY_OPERATOR(>>, __shift_right) 1195_DEFINE_BINARY_OPERATOR(&&, __logical_and) 1196_DEFINE_BINARY_OPERATOR(||, __logical_or) 1197_DEFINE_BINARY_OPERATOR(==, __equal_to) 1198_DEFINE_BINARY_OPERATOR(!=, __not_equal_to) 1199_DEFINE_BINARY_OPERATOR(<, __less) 1200_DEFINE_BINARY_OPERATOR(>, __greater) 1201_DEFINE_BINARY_OPERATOR(<=, __less_equal) 1202_DEFINE_BINARY_OPERATOR(>=, __greater_equal) 1203 1204#undef _DEFINE_BINARY_OPERATOR 1205 1206#if __cplusplus >= 201103L 1207 /** 1208 * @brief Return an iterator pointing to the first element of 1209 * the valarray. 1210 * @param __va valarray. 1211 */ 1212 template<class _Tp> 1213 inline _Tp* 1214 begin(valarray<_Tp>& __va) 1215 { return std::__addressof(__va[0]); } 1216 1217 /** 1218 * @brief Return an iterator pointing to the first element of 1219 * the const valarray. 1220 * @param __va valarray. 1221 */ 1222 template<class _Tp> 1223 inline const _Tp* 1224 begin(const valarray<_Tp>& __va) 1225 { return std::__addressof(__va[0]); } 1226 1227 /** 1228 * @brief Return an iterator pointing to one past the last element of 1229 * the valarray. 1230 * @param __va valarray. 1231 */ 1232 template<class _Tp> 1233 inline _Tp* 1234 end(valarray<_Tp>& __va) 1235 { return std::__addressof(__va[0]) + __va.size(); } 1236 1237 /** 1238 * @brief Return an iterator pointing to one past the last element of 1239 * the const valarray. 1240 * @param __va valarray. 1241 */ 1242 template<class _Tp> 1243 inline const _Tp* 1244 end(const valarray<_Tp>& __va) 1245 { return std::__addressof(__va[0]) + __va.size(); } 1246#endif // C++11 1247 1248 /// @} group numeric_arrays 1249 1250_GLIBCXX_END_NAMESPACE_VERSION 1251} // namespace 1252 1253#endif /* _GLIBCXX_VALARRAY */ 1254