1 //===----------------------------------------------------------------------===//
2 //
3 // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
4 // See https://llvm.org/LICENSE.txt for license information.
5 // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
6 //
7 //===----------------------------------------------------------------------===//
8
9 // <valarray>
10
11 // template<class T> class valarray;
12
13 // template<class T>
14 // valarray<T>
15 // pow(const valarray<T>& x, const valarray<T>& y);
16
17 #include <valarray>
18 #include <cassert>
19 #include <cstddef>
20
21 #include "test_macros.h"
22 #include "valarray_helper.h"
23
main(int,char **)24 int main(int, char**)
25 {
26 {
27 typedef double T;
28 T a1[] = {.9, .5, 0., .5, .75};
29 T a2[] = {-.8, .25, 0.375, -.5, .75};
30 T a3[] = {1.0879426248455297e+00,
31 8.4089641525371450e-01,
32 0.0000000000000000e+00,
33 1.4142135623730949e+00,
34 8.0592744886765644e-01};
35 const unsigned N = sizeof(a1)/sizeof(a1[0]);
36 std::valarray<T> v1(a1, N);
37 std::valarray<T> v2(a2, N);
38 std::valarray<T> v3 = pow(v1, v2);
39 assert(v3.size() == v1.size());
40 for (std::size_t i = 0; i < v3.size(); ++i)
41 assert(is_about(v3[i], a3[i], 10));
42 }
43
44 return 0;
45 }
46