15a83710eSEric Fiselier //===----------------------------------------------------------------------===//
25a83710eSEric Fiselier //
35a83710eSEric Fiselier //                     The LLVM Compiler Infrastructure
45a83710eSEric Fiselier //
55a83710eSEric Fiselier // This file is dual licensed under the MIT and the University of Illinois Open
65a83710eSEric Fiselier // Source Licenses. See LICENSE.TXT for details.
75a83710eSEric Fiselier //
85a83710eSEric Fiselier //===----------------------------------------------------------------------===//
95a83710eSEric Fiselier 
105a83710eSEric Fiselier // <complex>
115a83710eSEric Fiselier 
125a83710eSEric Fiselier // template<class T>
135a83710eSEric Fiselier //   T
145a83710eSEric Fiselier //   norm(const complex<T>& x);
155a83710eSEric Fiselier 
165a83710eSEric Fiselier #include <complex>
175a83710eSEric Fiselier #include <cassert>
185a83710eSEric Fiselier 
195a83710eSEric Fiselier #include "../cases.h"
205a83710eSEric Fiselier 
215a83710eSEric Fiselier template <class T>
225a83710eSEric Fiselier void
235a83710eSEric Fiselier test()
245a83710eSEric Fiselier {
255a83710eSEric Fiselier     std::complex<T> z(3, 4);
265a83710eSEric Fiselier     assert(norm(z) == 25);
275a83710eSEric Fiselier }
285a83710eSEric Fiselier 
295a83710eSEric Fiselier void test_edges()
305a83710eSEric Fiselier {
31*217cf69bSMarshall Clow     const unsigned N = sizeof(testcases) / sizeof(testcases[0]);
325a83710eSEric Fiselier     for (unsigned i = 0; i < N; ++i)
335a83710eSEric Fiselier     {
34*217cf69bSMarshall Clow         double r = norm(testcases[i]);
35*217cf69bSMarshall Clow         switch (classify(testcases[i]))
365a83710eSEric Fiselier         {
375a83710eSEric Fiselier         case zero:
385a83710eSEric Fiselier             assert(r == 0);
395a83710eSEric Fiselier             assert(!std::signbit(r));
405a83710eSEric Fiselier             break;
415a83710eSEric Fiselier         case non_zero:
425a83710eSEric Fiselier             assert(std::isfinite(r) && r > 0);
435a83710eSEric Fiselier             break;
445a83710eSEric Fiselier         case inf:
455a83710eSEric Fiselier             assert(std::isinf(r) && r > 0);
465a83710eSEric Fiselier             break;
475a83710eSEric Fiselier         case NaN:
485a83710eSEric Fiselier             assert(std::isnan(r));
495a83710eSEric Fiselier             break;
505a83710eSEric Fiselier         case non_zero_nan:
515a83710eSEric Fiselier             assert(std::isnan(r));
525a83710eSEric Fiselier             break;
535a83710eSEric Fiselier         }
545a83710eSEric Fiselier     }
555a83710eSEric Fiselier }
565a83710eSEric Fiselier 
575a83710eSEric Fiselier int main()
585a83710eSEric Fiselier {
595a83710eSEric Fiselier     test<float>();
605a83710eSEric Fiselier     test<double>();
615a83710eSEric Fiselier     test<long double>();
625a83710eSEric Fiselier     test_edges();
635a83710eSEric Fiselier }
64