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