15a83710eSEric Fiselier //===----------------------------------------------------------------------===//
25a83710eSEric Fiselier //
357b08b09SChandler Carruth // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
457b08b09SChandler Carruth // See https://llvm.org/LICENSE.txt for license information.
557b08b09SChandler Carruth // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
65a83710eSEric Fiselier //
75a83710eSEric Fiselier //===----------------------------------------------------------------------===//
85a83710eSEric Fiselier 
95a83710eSEric Fiselier // <complex>
105a83710eSEric Fiselier 
115a83710eSEric Fiselier // template<class T>
125a83710eSEric Fiselier //   T
135a83710eSEric Fiselier //   norm(const complex<T>& x);
145a83710eSEric Fiselier 
155a83710eSEric Fiselier #include <complex>
165a83710eSEric Fiselier #include <cassert>
175a83710eSEric Fiselier 
185a83710eSEric Fiselier #include "../cases.h"
195a83710eSEric Fiselier 
205a83710eSEric Fiselier template <class T>
215a83710eSEric Fiselier void
225a83710eSEric Fiselier test()
235a83710eSEric Fiselier {
245a83710eSEric Fiselier     std::complex<T> z(3, 4);
255a83710eSEric Fiselier     assert(norm(z) == 25);
265a83710eSEric Fiselier }
275a83710eSEric Fiselier 
285a83710eSEric Fiselier void test_edges()
295a83710eSEric Fiselier {
30217cf69bSMarshall Clow     const unsigned N = sizeof(testcases) / sizeof(testcases[0]);
315a83710eSEric Fiselier     for (unsigned i = 0; i < N; ++i)
325a83710eSEric Fiselier     {
33217cf69bSMarshall Clow         double r = norm(testcases[i]);
34217cf69bSMarshall Clow         switch (classify(testcases[i]))
355a83710eSEric Fiselier         {
365a83710eSEric Fiselier         case zero:
375a83710eSEric Fiselier             assert(r == 0);
385a83710eSEric Fiselier             assert(!std::signbit(r));
395a83710eSEric Fiselier             break;
405a83710eSEric Fiselier         case non_zero:
415a83710eSEric Fiselier             assert(std::isfinite(r) && r > 0);
425a83710eSEric Fiselier             break;
435a83710eSEric Fiselier         case inf:
445a83710eSEric Fiselier             assert(std::isinf(r) && r > 0);
455a83710eSEric Fiselier             break;
465a83710eSEric Fiselier         case NaN:
475a83710eSEric Fiselier             assert(std::isnan(r));
485a83710eSEric Fiselier             break;
495a83710eSEric Fiselier         case non_zero_nan:
505a83710eSEric Fiselier             assert(std::isnan(r));
515a83710eSEric Fiselier             break;
525a83710eSEric Fiselier         }
535a83710eSEric Fiselier     }
545a83710eSEric Fiselier }
555a83710eSEric Fiselier 
56*2df59c50SJF Bastien int main(int, char**)
575a83710eSEric Fiselier {
585a83710eSEric Fiselier     test<float>();
595a83710eSEric Fiselier     test<double>();
605a83710eSEric Fiselier     test<long double>();
615a83710eSEric Fiselier     test_edges();
62*2df59c50SJF Bastien 
63*2df59c50SJF Bastien   return 0;
645a83710eSEric Fiselier }
65