1bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant//===----------------------------------------------------------------------===// 2bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant// 3f5256e16dfc425c1d466f6308d4026d529ce9e0bHoward Hinnant// The LLVM Compiler Infrastructure 4bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant// 5b64f8b07c104c6cc986570ac8ee0ed16a9f23976Howard Hinnant// This file is dual licensed under the MIT and the University of Illinois Open 6b64f8b07c104c6cc986570ac8ee0ed16a9f23976Howard Hinnant// Source Licenses. See LICENSE.TXT for details. 7bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant// 8bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant//===----------------------------------------------------------------------===// 9bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 10bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant// <complex> 11bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 12bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant// template<class T> 13bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant// complex<T> 14bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant// cosh(const complex<T>& x); 15bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 16bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant#include <complex> 17bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant#include <cassert> 18bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 19bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant#include "../cases.h" 20bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 21bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnanttemplate <class T> 22bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnantvoid 23bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnanttest(const std::complex<T>& c, std::complex<T> x) 24bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant{ 25bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(cosh(c) == x); 26bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant} 27bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 28bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnanttemplate <class T> 29bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnantvoid 30bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnanttest() 31bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant{ 32bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant test(std::complex<T>(0, 0), std::complex<T>(1, 0)); 33bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant} 34bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 35bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnantvoid test_edges() 36bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant{ 37bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant typedef std::complex<double> C; 38bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant const double pi = std::atan2(+0., -0.); 39bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant const unsigned N = sizeof(x) / sizeof(x[0]); 40bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant for (unsigned i = 0; i < N; ++i) 41bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 42bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant std::complex<double> r = cosh(x[i]); 43bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant if (x[i].real() == 0 && x[i].imag() == 0) 44bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 45bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(r.real() == 1); 46bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(r.imag() == 0); 47bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::signbit(r.imag()) == std::signbit(x[i].imag())); 48bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 49bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (x[i].real() == 0 && std::isinf(x[i].imag())) 50bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 51bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.real())); 52bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(r.imag() == 0); 53bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 54bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (x[i].real() == 0 && std::isnan(x[i].imag())) 55bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 56bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.real())); 57bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(r.imag() == 0); 58bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 59bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (std::isfinite(x[i].real()) && std::isinf(x[i].imag())) 60bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 61bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.real())); 62bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.imag())); 63bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 64bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (std::isfinite(x[i].real()) && std::isnan(x[i].imag())) 65bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 66bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.real())); 67bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.imag())); 68bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 69bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (std::isinf(x[i].real()) && x[i].imag() == 0) 70bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 71bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isinf(r.real())); 72bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(!std::signbit(r.real())); 73bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(r.imag() == 0); 74bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::signbit(r.imag()) == std::signbit(x[i].imag())); 75bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 76bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (std::isinf(x[i].real()) && std::isfinite(x[i].imag())) 77bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 78bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isinf(r.real())); 79bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::signbit(r.real()) == std::signbit(cos(x[i].imag()))); 80bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isinf(r.imag())); 81bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::signbit(r.imag()) == std::signbit(x[i].real() * sin(x[i].imag()))); 82bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 83bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (std::isinf(x[i].real()) && std::isinf(x[i].imag())) 84bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 85bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isinf(r.real())); 86bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.imag())); 87bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 88bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (std::isinf(x[i].real()) && std::isnan(x[i].imag())) 89bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 90bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isinf(r.real())); 91bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(r.real() > 0); 92bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.imag())); 93bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 94bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (std::isnan(x[i].real()) && x[i].imag() == 0) 95bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 96bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.real())); 97bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(r.imag() == 0); 98bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 99bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (std::isnan(x[i].real()) && std::isfinite(x[i].imag())) 100bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 101bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.real())); 102bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.imag())); 103bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 104bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else if (std::isnan(x[i].real()) && std::isnan(x[i].imag())) 105bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 106bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.real())); 107bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(r.imag())); 108bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 109bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 110bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant} 111bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 112bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnantint main() 113bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant{ 114bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant test<float>(); 115bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant test<double>(); 116bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant test<long double>(); 117bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant test_edges(); 118bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant} 119