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 // XFAIL: LIBCXX-AIX-FIXME 10 11 // <complex> 12 13 // template<class T> 14 // complex<T> 15 // asin(const complex<T>& x); 16 17 #include <complex> 18 #include <cassert> 19 20 #include "test_macros.h" 21 #include "../cases.h" 22 23 template <class T> 24 void 25 test(const std::complex<T>& c, std::complex<T> x) 26 { 27 assert(asin(c) == x); 28 } 29 30 template <class T> 31 void 32 test() 33 { 34 test(std::complex<T>(0, 0), std::complex<T>(0, 0)); 35 } 36 37 void test_edges() 38 { 39 const double pi = std::atan2(+0., -0.); 40 const unsigned N = sizeof(testcases) / sizeof(testcases[0]); 41 for (unsigned i = 0; i < N; ++i) 42 { 43 std::complex<double> r = asin(testcases[i]); 44 if (testcases[i].real() == 0 && testcases[i].imag() == 0) 45 { 46 assert(std::signbit(r.real()) == std::signbit(testcases[i].real())); 47 assert(std::signbit(r.imag()) == std::signbit(testcases[i].imag())); 48 } 49 else if (std::isfinite(testcases[i].real()) && std::isinf(testcases[i].imag())) 50 { 51 assert(r.real() == 0); 52 assert(std::signbit(testcases[i].real()) == std::signbit(r.real())); 53 assert(std::isinf(r.imag())); 54 assert(std::signbit(testcases[i].imag()) == std::signbit(r.imag())); 55 } 56 else if ( testcases[i].real() == 0 && std::isnan(testcases[i].imag())) 57 { 58 assert(r.real() == 0); 59 assert(std::signbit(testcases[i].real()) == std::signbit(r.real())); 60 assert(std::isnan(r.imag())); 61 } 62 else if (std::isfinite(testcases[i].real()) && std::isnan(testcases[i].imag())) 63 { 64 assert(std::isnan(r.real())); 65 assert(std::isnan(r.imag())); 66 } 67 else if (std::isinf(testcases[i].real()) && std::isfinite(testcases[i].imag())) 68 { 69 if (testcases[i].real() > 0) 70 is_about(r.real(), pi/2); 71 else 72 is_about(r.real(), - pi/2); 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()) && std::isinf(testcases[i].imag())) 77 { 78 if (std::signbit(testcases[i].real())) 79 is_about(r.real(), -pi/4); 80 else 81 is_about(r.real(), pi/4); 82 assert(std::isinf(r.imag())); 83 assert(std::signbit(testcases[i].imag()) == std::signbit(r.imag())); 84 } 85 else if (std::isinf(testcases[i].real()) && std::isnan(testcases[i].imag())) 86 { 87 assert(std::isnan(r.real())); 88 assert(std::isinf(r.imag())); 89 assert(std::signbit(testcases[i].real()) != std::signbit(r.imag())); 90 } 91 else if (std::isnan(testcases[i].real()) && std::isfinite(testcases[i].imag())) 92 { 93 assert(std::isnan(r.real())); 94 assert(std::isnan(r.imag())); 95 } 96 else if (std::isnan(testcases[i].real()) && std::isinf(testcases[i].imag())) 97 { 98 assert(std::isnan(r.real())); 99 assert(std::isinf(r.imag())); 100 } 101 else if (std::isnan(testcases[i].real()) && std::isnan(testcases[i].imag())) 102 { 103 assert(std::isnan(r.real())); 104 assert(std::isnan(r.imag())); 105 } 106 else 107 { 108 assert(std::signbit(r.real()) == std::signbit(testcases[i].real())); 109 assert(std::signbit(r.imag()) == std::signbit(testcases[i].imag())); 110 } 111 } 112 } 113 114 int main(int, char**) 115 { 116 test<float>(); 117 test<double>(); 118 test<long double>(); 119 test_edges(); 120 121 return 0; 122 } 123