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 // <complex> 11 12 // template<class T> 13 // complex<T> 14 // pow(const complex<T>& x, const T& y); 15 16 #include <complex> 17 #include <cassert> 18 19 #include "../cases.h" 20 21 template <class T> 22 void 23 test(const std::complex<T>& a, const T& b, std::complex<T> x) 24 { 25 std::complex<T> c = pow(a, b); 26 is_about(real(c), real(x)); 27 is_about(imag(c), imag(x)); 28 } 29 30 template <class T> 31 void 32 test() 33 { 34 test(std::complex<T>(2, 3), T(2), std::complex<T>(-5, 12)); 35 } 36 37 void test_edges() 38 { 39 const unsigned N = sizeof(testcases) / sizeof(testcases[0]); 40 for (unsigned i = 0; i < N; ++i) 41 { 42 for (unsigned j = 0; j < N; ++j) 43 { 44 std::complex<double> r = pow(testcases[i], real(testcases[j])); 45 std::complex<double> z = exp(std::complex<double>(real(testcases[j])) * log(testcases[i])); 46 if (std::isnan(real(r))) 47 assert(std::isnan(real(z))); 48 else 49 { 50 assert(real(r) == real(z)); 51 } 52 if (std::isnan(imag(r))) 53 assert(std::isnan(imag(z))); 54 else 55 { 56 assert(imag(r) == imag(z)); 57 } 58 } 59 } 60 } 61 62 int main() 63 { 64 test<float>(); 65 test<double>(); 66 test<long double>(); 67 test_edges(); 68 } 69
