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 valarray;
13
14 // template<class T>
15 // valarray<T>
16 // pow(const valarray<T>& x, const valarray<T>& y);
17
18 #include <valarray>
19 #include <cassert>
20 #include <sstream>
21 #include <cstddef>
22
is_about(double x,double y,int p)23 bool is_about(double x, double y, int p)
24 {
25 std::ostringstream o;
26 o.precision(p);
27 scientific(o);
28 o << x;
29 std::string a = o.str();
30 o.str("");
31 o << y;
32 return a == o.str();
33 }
34
main()35 int main()
36 {
37 {
38 typedef double T;
39 T a1[] = {.9, .5, 0., .5, .75};
40 T a2[] = {-.8, .25, 0.375, -.5, .75};
41 T a3[] = {1.0879426248455297e+00,
42 8.4089641525371450e-01,
43 0.0000000000000000e+00,
44 1.4142135623730949e+00,
45 8.0592744886765644e-01};
46 const unsigned N = sizeof(a1)/sizeof(a1[0]);
47 std::valarray<T> v1(a1, N);
48 std::valarray<T> v2(a2, N);
49 std::valarray<T> v3 = pow(v1, v2);
50 assert(v3.size() == v1.size());
51 for (std::size_t i = 0; i < v3.size(); ++i)
52 assert(is_about(v3[i], a3[i], 10));
53 }
54 }
55