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 // sinh(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(sinh(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 { 37 typedef std::complex<double> C; 38 const unsigned N = sizeof(testcases) / sizeof(testcases[0]); 39 for (unsigned i = 0; i < N; ++i) 40 { 41 std::complex<double> r = sinh(testcases[i]); 42 if (testcases[i].real() == 0 && testcases[i].imag() == 0) 43 { 44 assert(r.real() == 0); 45 assert(std::signbit(r.real()) == std::signbit(testcases[i].real())); 46 assert(r.imag() == 0); 47 assert(std::signbit(r.imag()) == std::signbit(testcases[i].imag())); 48 } 49 else if (testcases[i].real() == 0 && std::isinf(testcases[i].imag())) 50 { 51 assert(r.real() == 0); 52 assert(std::isnan(r.imag())); 53 } 54 else if (std::isfinite(testcases[i].real()) && std::isinf(testcases[i].imag())) 55 { 56 assert(std::isnan(r.real())); 57 assert(std::isnan(r.imag())); 58 } 59 else if (testcases[i].real() == 0 && std::isnan(testcases[i].imag())) 60 { 61 assert(r.real() == 0); 62 assert(std::isnan(r.imag())); 63 } 64 else if (std::isfinite(testcases[i].real()) && std::isnan(testcases[i].imag())) 65 { 66 assert(std::isnan(r.real())); 67 assert(std::isnan(r.imag())); 68 } 69 else if (std::isinf(testcases[i].real()) && testcases[i].imag() == 0) 70 { 71 assert(std::isinf(r.real())); 72 assert(std::signbit(r.real()) == std::signbit(testcases[i].real())); 73 assert(r.imag() == 0); 74 assert(std::signbit(r.imag()) == std::signbit(testcases[i].imag())); 75 } 76 else if (std::isinf(testcases[i].real()) && std::isfinite(testcases[i].imag())) 77 { 78 assert(std::isinf(r.real())); 79 assert(std::signbit(r.real()) == std::signbit(testcases[i].real() * cos(testcases[i].imag()))); 80 assert(std::isinf(r.imag())); 81 assert(std::signbit(r.imag()) == std::signbit(sin(testcases[i].imag()))); 82 } 83 else if (std::isinf(testcases[i].real()) && std::isinf(testcases[i].imag())) 84 { 85 assert(std::isinf(r.real())); 86 assert(std::isnan(r.imag())); 87 } 88 else if (std::isinf(testcases[i].real()) && std::isnan(testcases[i].imag())) 89 { 90 assert(std::isinf(r.real())); 91 assert(std::isnan(r.imag())); 92 } 93 else if (std::isnan(testcases[i].real()) && testcases[i].imag() == 0) 94 { 95 assert(std::isnan(r.real())); 96 assert(r.imag() == 0); 97 assert(std::signbit(r.imag()) == std::signbit(testcases[i].imag())); 98 } 99 else if (std::isnan(testcases[i].real()) && std::isfinite(testcases[i].imag())) 100 { 101 assert(std::isnan(r.real())); 102 assert(std::isnan(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 } 110 } 111 112 int main() 113 { 114 test<float>(); 115 test<double>(); 116 test<long double>(); 117 test_edges(); 118 } 119