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 120e20cae1a5be18fba591cd884aa2a389b66a3f49Howard Hinnant// template<class T> 13bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant// complex<T> 14bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant// pow(const complex<T>& x, const complex<T>& y); 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>& a, const std::complex<T>& b, std::complex<T> x) 24bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant{ 25bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant std::complex<T> c = pow(a, b); 26bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant is_about(real(c), real(x)); 27bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant is_about(imag(c), imag(x)); 28bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant} 29bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 30bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnanttemplate <class T> 31bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnantvoid 32bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnanttest() 33bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant{ 34bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant test(std::complex<T>(2, 3), std::complex<T>(2, 0), std::complex<T>(-5, 12)); 35bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant} 36bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 37bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnantvoid test_edges() 38bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant{ 39bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant const double pi = std::atan2(+0., -0.); 40bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant const unsigned N = sizeof(x) / sizeof(x[0]); 41bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant for (unsigned i = 0; i < N; ++i) 42bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 43bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant for (unsigned j = 0; j < N; ++j) 44bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 45bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant std::complex<double> r = pow(x[i], x[j]); 46bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant std::complex<double> z = exp(x[j] * log(x[i])); 47bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant if (std::isnan(real(r))) 48bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(real(z))); 49bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else 50bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 51bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(real(r) == real(z)); 52bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::signbit(real(r)) == std::signbit(real(z))); 53bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 54bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant if (std::isnan(imag(r))) 55bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::isnan(imag(z))); 56bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant else 57bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant { 58bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(imag(r) == imag(z)); 59bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant assert(std::signbit(imag(r)) == std::signbit(imag(z))); 60bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 61bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 62bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant } 63bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant} 64bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant 65bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnantint main() 66bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant{ 67bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant test<float>(); 68bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant test<double>(); 69bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant test<long double>(); 70bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant test_edges(); 71bc8d3f97eb5c958007f2713238472e0c1c8fe02Howard Hinnant} 72