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 //   acos(const complex<T>& x);
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>& c, std::complex<T> x)
24 {
25     assert(acos(c) == x);
26 }
27 
28 template <class T>
29 void
30 test()
31 {
32     test(std::complex<T>(INFINITY, 1), std::complex<T>(0, -INFINITY));
33 }
34 
35 void test_edges()
36 {
37     typedef std::complex<double> C;
38     const double pi = std::atan2(+0., -0.);
39     const unsigned N = sizeof(x) / sizeof(x[0]);
40     for (unsigned i = 0; i < N; ++i)
41     {
42         std::complex<double> r = acos(x[i]);
43         if (x[i].real() == 0 && x[i].imag() == 0)
44         {
45             is_about(r.real(), pi/2);
46             assert(r.imag() == 0);
47             assert(std::signbit(x[i].imag()) != std::signbit(r.imag()));
48         }
49         else if (x[i].real() == 0 && std::isnan(x[i].imag()))
50         {
51             is_about(r.real(), pi/2);
52             assert(std::isnan(r.imag()));
53         }
54         else if (std::isfinite(x[i].real()) && std::isinf(x[i].imag()))
55         {
56             is_about(r.real(), pi/2);
57             assert(std::isinf(r.imag()));
58             assert(std::signbit(x[i].imag()) != std::signbit(r.imag()));
59         }
60         else if (std::isfinite(x[i].real()) && x[i].real() != 0 && std::isnan(x[i].imag()))
61         {
62             assert(std::isnan(r.real()));
63             assert(std::isnan(r.imag()));
64         }
65         else if (std::isinf(x[i].real()) && x[i].real() < 0 && std::isfinite(x[i].imag()))
66         {
67             is_about(r.real(), pi);
68             assert(std::isinf(r.imag()));
69             assert(std::signbit(x[i].imag()) != std::signbit(r.imag()));
70         }
71         else if (std::isinf(x[i].real()) && x[i].real() > 0 && std::isfinite(x[i].imag()))
72         {
73             assert(r.real() == 0);
74             assert(!std::signbit(r.real()));
75             assert(std::isinf(r.imag()));
76             assert(std::signbit(x[i].imag()) != std::signbit(r.imag()));
77         }
78         else if (std::isinf(x[i].real()) && x[i].real() < 0 && std::isinf(x[i].imag()))
79         {
80             is_about(r.real(), 0.75 * pi);
81             assert(std::isinf(r.imag()));
82             assert(std::signbit(x[i].imag()) != std::signbit(r.imag()));
83         }
84         else if (std::isinf(x[i].real()) && x[i].real() > 0 && std::isinf(x[i].imag()))
85         {
86             is_about(r.real(), 0.25 * pi);
87             assert(std::isinf(r.imag()));
88             assert(std::signbit(x[i].imag()) != std::signbit(r.imag()));
89         }
90         else if (std::isinf(x[i].real()) && std::isnan(x[i].imag()))
91         {
92             assert(std::isnan(r.real()));
93             assert(std::isinf(r.imag()));
94         }
95         else if (std::isnan(x[i].real()) && std::isfinite(x[i].imag()))
96         {
97             assert(std::isnan(r.real()));
98             assert(std::isnan(r.imag()));
99         }
100         else if (std::isnan(x[i].real()) && std::isinf(x[i].imag()))
101         {
102             assert(std::isnan(r.real()));
103             assert(std::isinf(r.imag()));
104             assert(std::signbit(x[i].imag()) != std::signbit(r.imag()));
105         }
106         else if (std::isnan(x[i].real()) && std::isnan(x[i].imag()))
107         {
108             assert(std::isnan(r.real()));
109             assert(std::isnan(r.imag()));
110         }
111         else if (!std::signbit(x[i].real()) && !std::signbit(x[i].imag()))
112         {
113             assert(!std::signbit(r.real()));
114             assert( std::signbit(r.imag()));
115         }
116         else if (std::signbit(x[i].real()) && !std::signbit(x[i].imag()))
117         {
118             assert(!std::signbit(r.real()));
119             assert( std::signbit(r.imag()));
120         }
121         else if (std::signbit(x[i].real()) && std::signbit(x[i].imag()))
122         {
123             assert(!std::signbit(r.real()));
124             assert(!std::signbit(r.imag()));
125         }
126         else if (!std::signbit(x[i].real()) && std::signbit(x[i].imag()))
127         {
128             assert(!std::signbit(r.real()));
129             assert(!std::signbit(r.imag()));
130         }
131     }
132 }
133 
134 int main()
135 {
136     test<float>();
137     test<double>();
138     test<long double>();
139     test_edges();
140 }
141