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