1*5a83710eSEric Fiselier //===----------------------------------------------------------------------===// 2*5a83710eSEric Fiselier // 3*5a83710eSEric Fiselier // The LLVM Compiler Infrastructure 4*5a83710eSEric Fiselier // 5*5a83710eSEric Fiselier // This file is dual licensed under the MIT and the University of Illinois Open 6*5a83710eSEric Fiselier // Source Licenses. See LICENSE.TXT for details. 7*5a83710eSEric Fiselier // 8*5a83710eSEric Fiselier //===----------------------------------------------------------------------===// 9*5a83710eSEric Fiselier 10*5a83710eSEric Fiselier // <complex> 11*5a83710eSEric Fiselier 12*5a83710eSEric Fiselier // template<class T> 13*5a83710eSEric Fiselier // T 14*5a83710eSEric Fiselier // norm(const complex<T>& x); 15*5a83710eSEric Fiselier 16*5a83710eSEric Fiselier #include <complex> 17*5a83710eSEric Fiselier #include <cassert> 18*5a83710eSEric Fiselier 19*5a83710eSEric Fiselier #include "../cases.h" 20*5a83710eSEric Fiselier 21*5a83710eSEric Fiselier template <class T> 22*5a83710eSEric Fiselier void 23*5a83710eSEric Fiselier test() 24*5a83710eSEric Fiselier { 25*5a83710eSEric Fiselier std::complex<T> z(3, 4); 26*5a83710eSEric Fiselier assert(norm(z) == 25); 27*5a83710eSEric Fiselier } 28*5a83710eSEric Fiselier 29*5a83710eSEric Fiselier void test_edges() 30*5a83710eSEric Fiselier { 31*5a83710eSEric Fiselier const unsigned N = sizeof(x) / sizeof(x[0]); 32*5a83710eSEric Fiselier for (unsigned i = 0; i < N; ++i) 33*5a83710eSEric Fiselier { 34*5a83710eSEric Fiselier double r = norm(x[i]); 35*5a83710eSEric Fiselier switch (classify(x[i])) 36*5a83710eSEric Fiselier { 37*5a83710eSEric Fiselier case zero: 38*5a83710eSEric Fiselier assert(r == 0); 39*5a83710eSEric Fiselier assert(!std::signbit(r)); 40*5a83710eSEric Fiselier break; 41*5a83710eSEric Fiselier case non_zero: 42*5a83710eSEric Fiselier assert(std::isfinite(r) && r > 0); 43*5a83710eSEric Fiselier break; 44*5a83710eSEric Fiselier case inf: 45*5a83710eSEric Fiselier assert(std::isinf(r) && r > 0); 46*5a83710eSEric Fiselier break; 47*5a83710eSEric Fiselier case NaN: 48*5a83710eSEric Fiselier assert(std::isnan(r)); 49*5a83710eSEric Fiselier break; 50*5a83710eSEric Fiselier case non_zero_nan: 51*5a83710eSEric Fiselier assert(std::isnan(r)); 52*5a83710eSEric Fiselier break; 53*5a83710eSEric Fiselier } 54*5a83710eSEric Fiselier } 55*5a83710eSEric Fiselier } 56*5a83710eSEric Fiselier 57*5a83710eSEric Fiselier int main() 58*5a83710eSEric Fiselier { 59*5a83710eSEric Fiselier test<float>(); 60*5a83710eSEric Fiselier test<double>(); 61*5a83710eSEric Fiselier test<long double>(); 62*5a83710eSEric Fiselier test_edges(); 63*5a83710eSEric Fiselier } 64