46 lines · cpp
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-exception6//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 24int 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