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