| 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 | |
|---|
| 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(v: v1, w: 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 | |
|---|
