1//===----------------------------------------------------------------------===//
2//
3//                     The LLVM Compiler Infrastructure
4//
5// This file is dual licensed under the MIT and the University of Illinois Open
6// Source Licenses. See LICENSE.TXT for details.
7//
8//===----------------------------------------------------------------------===//
9
10// <valarray>
11
12// template <class T> class indirect_array
13
14// void operator*=(const valarray<value_type>& v) const;
15
16#include <valarray>
17#include <cassert>
18
19int main()
20{
21    int a1[] = { 0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11,
22                12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23,
23                24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,
24                36, 37, 38, 39, 40};
25    const std::size_t N1 = sizeof(a1)/sizeof(a1[0]);
26    int a2[] = {  1,   2,   3,   4,   5,   6,   7,   8,   9,  10,  11, 12,
27                 13,  14,  15,  16,  17,  18,  19,  20,  21,  22,  23, 24};
28    const std::size_t N2 = sizeof(a2)/sizeof(a2[0]);
29    std::size_t s[N2] = { 3,  4,  5,  7,  8,  9, 11, 12, 13, 15, 16, 17,
30                         22, 23, 24, 26, 27, 28, 30, 31, 32, 34, 35, 36};
31    std::valarray<int> v1(a1, N1);
32    std::valarray<int> v2(a2, N2);
33    std::valarray<std::size_t> ia(s, N2);
34    v1[ia] *= v2;
35    assert(v1.size() == 41);
36    assert(v1[ 0] ==  0);
37    assert(v1[ 1] ==  1);
38    assert(v1[ 2] ==  2);
39    assert(v1[ 3] ==  3);
40    assert(v1[ 4] ==  8);
41    assert(v1[ 5] == 15);
42    assert(v1[ 6] ==  6);
43    assert(v1[ 7] == 28);
44    assert(v1[ 8] == 40);
45    assert(v1[ 9] == 54);
46    assert(v1[10] == 10);
47    assert(v1[11] == 77);
48    assert(v1[12] == 96);
49    assert(v1[13] == 117);
50    assert(v1[14] == 14);
51    assert(v1[15] == 150);
52    assert(v1[16] == 176);
53    assert(v1[17] == 204);
54    assert(v1[18] == 18);
55    assert(v1[19] == 19);
56    assert(v1[20] == 20);
57    assert(v1[21] == 21);
58    assert(v1[22] == 286);
59    assert(v1[23] == 322);
60    assert(v1[24] == 360);
61    assert(v1[25] == 25);
62    assert(v1[26] == 416);
63    assert(v1[27] == 459);
64    assert(v1[28] == 504);
65    assert(v1[29] == 29);
66    assert(v1[30] == 570);
67    assert(v1[31] == 620);
68    assert(v1[32] == 672);
69    assert(v1[33] == 33);
70    assert(v1[34] == 748);
71    assert(v1[35] == 805);
72    assert(v1[36] == 864);
73    assert(v1[37] == 37);
74    assert(v1[38] == 38);
75    assert(v1[39] == 39);
76    assert(v1[40] == 40);
77}
78