//===----------------------------------------------------------------------===// // // The LLVM Compiler Infrastructure // // This file is dual licensed under the MIT and the University of Illinois Open // Source Licenses. See LICENSE.TXT for details. // //===----------------------------------------------------------------------===// // // template class gslice_array // void operator*= (const valarray& v) const; #include #include int main() { int a1[] = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40}; int a2[] = { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24}; std::valarray v1(a1, sizeof(a1)/sizeof(a1[0])); std::valarray v2(a2, sizeof(a2)/sizeof(a2[0])); std::size_t sz[] = {2, 4, 3}; std::size_t st[] = {19, 4, 1}; typedef std::valarray sizes; typedef std::valarray strides; v1[std::gslice(3, sizes(sz, sizeof(sz)/sizeof(sz[0])), strides(st, sizeof(st)/sizeof(st[0])))] *= v2; assert(v1.size() == 41); assert(v1[ 0] == 0); assert(v1[ 1] == 1); assert(v1[ 2] == 2); assert(v1[ 3] == 3); assert(v1[ 4] == 8); assert(v1[ 5] == 15); assert(v1[ 6] == 6); assert(v1[ 7] == 28); assert(v1[ 8] == 40); assert(v1[ 9] == 54); assert(v1[10] == 10); assert(v1[11] == 77); assert(v1[12] == 96); assert(v1[13] == 117); assert(v1[14] == 14); assert(v1[15] == 150); assert(v1[16] == 176); assert(v1[17] == 204); assert(v1[18] == 18); assert(v1[19] == 19); assert(v1[20] == 20); assert(v1[21] == 21); assert(v1[22] == 286); assert(v1[23] == 322); assert(v1[24] == 360); assert(v1[25] == 25); assert(v1[26] == 416); assert(v1[27] == 459); assert(v1[28] == 504); assert(v1[29] == 29); assert(v1[30] == 570); assert(v1[31] == 620); assert(v1[32] == 672); assert(v1[33] == 33); assert(v1[34] == 748); assert(v1[35] == 805); assert(v1[36] == 864); assert(v1[37] == 37); assert(v1[38] == 38); assert(v1[39] == 39); assert(v1[40] == 40); }