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