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