Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 1 | //===----------------------------------------------------------------------===// |
| 2 | // |
Howard Hinnant | f5256e1 | 2010-05-11 21:36:01 +0000 | [diff] [blame] | 3 | // The LLVM Compiler Infrastructure |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 4 | // |
Howard Hinnant | b64f8b0 | 2010-11-16 22:09:02 +0000 | [diff] [blame] | 5 | // This file is dual licensed under the MIT and the University of Illinois Open |
| 6 | // Source Licenses. See LICENSE.TXT for details. |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 7 | // |
| 8 | //===----------------------------------------------------------------------===// |
| 9 | |
| 10 | // <complex> |
| 11 | |
| 12 | // template<class T> |
| 13 | // complex<T> |
| 14 | // asinh(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(asinh(c) == x); |
| 26 | } |
| 27 | |
| 28 | template <class T> |
| 29 | void |
| 30 | test() |
| 31 | { |
| 32 | test(std::complex<T>(0, 0), std::complex<T>(0, 0)); |
| 33 | } |
| 34 | |
| 35 | void test_edges() |
| 36 | { |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 37 | const double pi = std::atan2(+0., -0.); |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 38 | const unsigned N = sizeof(testcases) / sizeof(testcases[0]); |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 39 | for (unsigned i = 0; i < N; ++i) |
| 40 | { |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 41 | std::complex<double> r = asinh(testcases[i]); |
| 42 | if (testcases[i].real() == 0 && testcases[i].imag() == 0) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 43 | { |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 44 | assert(std::signbit(r.real()) == std::signbit(testcases[i].real())); |
| 45 | assert(std::signbit(r.imag()) == std::signbit(testcases[i].imag())); |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 46 | } |
Mikhail Maltsev | c4658ab | 2018-02-19 15:41:36 +0000 | [diff] [blame^] | 47 | else if (testcases[i].real() == 0 && std::abs(testcases[i].imag()) == 1) |
| 48 | { |
| 49 | assert(r.real() == 0); |
| 50 | assert(std::signbit(testcases[i].imag()) == std::signbit(r.imag())); |
| 51 | if (std::signbit(testcases[i].imag())) |
| 52 | is_about(r.imag(), -pi/2); |
| 53 | else |
| 54 | is_about(r.imag(), pi/2); |
| 55 | } |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 56 | else if (std::isfinite(testcases[i].real()) && std::isinf(testcases[i].imag())) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 57 | { |
| 58 | assert(std::isinf(r.real())); |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 59 | assert(std::signbit(testcases[i].real()) == std::signbit(r.real())); |
| 60 | if (std::signbit(testcases[i].imag())) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 61 | is_about(r.imag(), -pi/2); |
| 62 | else |
| 63 | is_about(r.imag(), pi/2); |
| 64 | } |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 65 | else if (std::isfinite(testcases[i].real()) && std::isnan(testcases[i].imag())) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 66 | { |
| 67 | assert(std::isnan(r.real())); |
| 68 | assert(std::isnan(r.imag())); |
| 69 | } |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 70 | else if (std::isinf(testcases[i].real()) && std::isfinite(testcases[i].imag())) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 71 | { |
| 72 | assert(std::isinf(r.real())); |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 73 | assert(std::signbit(testcases[i].real()) == std::signbit(r.real())); |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 74 | assert(r.imag() == 0); |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 75 | assert(std::signbit(testcases[i].imag()) == std::signbit(r.imag())); |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 76 | } |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 77 | else if (std::isinf(testcases[i].real()) && std::isinf(testcases[i].imag())) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 78 | { |
| 79 | assert(std::isinf(r.real())); |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 80 | assert(std::signbit(testcases[i].real()) == std::signbit(r.real())); |
| 81 | if (std::signbit(testcases[i].imag())) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 82 | is_about(r.imag(), -pi/4); |
| 83 | else |
| 84 | is_about(r.imag(), pi/4); |
| 85 | } |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 86 | else if (std::isinf(testcases[i].real()) && std::isnan(testcases[i].imag())) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 87 | { |
| 88 | assert(std::isinf(r.real())); |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 89 | assert(std::signbit(testcases[i].real()) == std::signbit(r.real())); |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 90 | assert(std::isnan(r.imag())); |
| 91 | } |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 92 | else if (std::isnan(testcases[i].real()) && testcases[i].imag() == 0) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 93 | { |
| 94 | assert(std::isnan(r.real())); |
| 95 | assert(r.imag() == 0); |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 96 | assert(std::signbit(testcases[i].imag()) == std::signbit(r.imag())); |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 97 | } |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 98 | else if (std::isnan(testcases[i].real()) && std::isfinite(testcases[i].imag())) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 99 | { |
| 100 | assert(std::isnan(r.real())); |
| 101 | assert(std::isnan(r.imag())); |
| 102 | } |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 103 | else if (std::isnan(testcases[i].real()) && std::isinf(testcases[i].imag())) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 104 | { |
| 105 | assert(std::isinf(r.real())); |
| 106 | assert(std::isnan(r.imag())); |
| 107 | } |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 108 | else if (std::isnan(testcases[i].real()) && std::isnan(testcases[i].imag())) |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 109 | { |
| 110 | assert(std::isnan(r.real())); |
| 111 | assert(std::isnan(r.imag())); |
| 112 | } |
| 113 | else |
| 114 | { |
Marshall Clow | fb59f5c | 2016-06-06 18:24:08 +0000 | [diff] [blame] | 115 | assert(std::signbit(r.real()) == std::signbit(testcases[i].real())); |
| 116 | assert(std::signbit(r.imag()) == std::signbit(testcases[i].imag())); |
Howard Hinnant | bc8d3f9 | 2010-05-11 19:42:16 +0000 | [diff] [blame] | 117 | } |
| 118 | } |
| 119 | } |
| 120 | |
| 121 | int main() |
| 122 | { |
| 123 | test<float>(); |
| 124 | test<double>(); |
| 125 | test<long double>(); |
| 126 | test_edges(); |
| 127 | } |