diff options
Diffstat (limited to 'test/std/numerics/complex.number/cases.h')
| -rw-r--r-- | test/std/numerics/complex.number/cases.h | 230 |
1 files changed, 230 insertions, 0 deletions
diff --git a/test/std/numerics/complex.number/cases.h b/test/std/numerics/complex.number/cases.h new file mode 100644 index 0000000000000..cd6bd165a9d16 --- /dev/null +++ b/test/std/numerics/complex.number/cases.h @@ -0,0 +1,230 @@ +//===----------------------------------------------------------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// + +// <complex> + +// test cases + +#ifndef CASES_H +#define CASES_H + +#include <complex> +#include <cassert> + +std::complex<double> x[] = +{ + std::complex<double>( 1.e-6, 1.e-6), + std::complex<double>(-1.e-6, 1.e-6), + std::complex<double>(-1.e-6, -1.e-6), + std::complex<double>( 1.e-6, -1.e-6), + + std::complex<double>( 1.e+6, 1.e-6), + std::complex<double>(-1.e+6, 1.e-6), + std::complex<double>(-1.e+6, -1.e-6), + std::complex<double>( 1.e+6, -1.e-6), + + std::complex<double>( 1.e-6, 1.e+6), + std::complex<double>(-1.e-6, 1.e+6), + std::complex<double>(-1.e-6, -1.e+6), + std::complex<double>( 1.e-6, -1.e+6), + + std::complex<double>( 1.e+6, 1.e+6), + std::complex<double>(-1.e+6, 1.e+6), + std::complex<double>(-1.e+6, -1.e+6), + std::complex<double>( 1.e+6, -1.e+6), + + std::complex<double>(NAN, NAN), + std::complex<double>(-INFINITY, NAN), + std::complex<double>(-2, NAN), + std::complex<double>(-1, NAN), + std::complex<double>(-0.5, NAN), + std::complex<double>(-0., NAN), + std::complex<double>(+0., NAN), + std::complex<double>(0.5, NAN), + std::complex<double>(1, NAN), + std::complex<double>(2, NAN), + std::complex<double>(INFINITY, NAN), + + std::complex<double>(NAN, -INFINITY), + std::complex<double>(-INFINITY, -INFINITY), + std::complex<double>(-2, -INFINITY), + std::complex<double>(-1, -INFINITY), + std::complex<double>(-0.5, -INFINITY), + std::complex<double>(-0., -INFINITY), + std::complex<double>(+0., -INFINITY), + std::complex<double>(0.5, -INFINITY), + std::complex<double>(1, -INFINITY), + std::complex<double>(2, -INFINITY), + std::complex<double>(INFINITY, -INFINITY), + + std::complex<double>(NAN, -2), + std::complex<double>(-INFINITY, -2), + std::complex<double>(-2, -2), + std::complex<double>(-1, -2), + std::complex<double>(-0.5, -2), + std::complex<double>(-0., -2), + std::complex<double>(+0., -2), + std::complex<double>(0.5, -2), + std::complex<double>(1, -2), + std::complex<double>(2, -2), + std::complex<double>(INFINITY, -2), + + std::complex<double>(NAN, -1), + std::complex<double>(-INFINITY, -1), + std::complex<double>(-2, -1), + std::complex<double>(-1, -1), + std::complex<double>(-0.5, -1), + std::complex<double>(-0., -1), + std::complex<double>(+0., -1), + std::complex<double>(0.5, -1), + std::complex<double>(1, -1), + std::complex<double>(2, -1), + std::complex<double>(INFINITY, -1), + + std::complex<double>(NAN, -0.5), + std::complex<double>(-INFINITY, -0.5), + std::complex<double>(-2, -0.5), + std::complex<double>(-1, -0.5), + std::complex<double>(-0.5, -0.5), + std::complex<double>(-0., -0.5), + std::complex<double>(+0., -0.5), + std::complex<double>(0.5, -0.5), + std::complex<double>(1, -0.5), + std::complex<double>(2, -0.5), + std::complex<double>(INFINITY, -0.5), + + std::complex<double>(NAN, -0.), + std::complex<double>(-INFINITY, -0.), + std::complex<double>(-2, -0.), + std::complex<double>(-1, -0.), + std::complex<double>(-0.5, -0.), + std::complex<double>(-0., -0.), + std::complex<double>(+0., -0.), + std::complex<double>(0.5, -0.), + std::complex<double>(1, -0.), + std::complex<double>(2, -0.), + std::complex<double>(INFINITY, -0.), + + std::complex<double>(NAN, +0.), + std::complex<double>(-INFINITY, +0.), + std::complex<double>(-2, +0.), + std::complex<double>(-1, +0.), + std::complex<double>(-0.5, +0.), + std::complex<double>(-0., +0.), + std::complex<double>(+0., +0.), + std::complex<double>(0.5, +0.), + std::complex<double>(1, +0.), + std::complex<double>(2, +0.), + std::complex<double>(INFINITY, +0.), + + std::complex<double>(NAN, 0.5), + std::complex<double>(-INFINITY, 0.5), + std::complex<double>(-2, 0.5), + std::complex<double>(-1, 0.5), + std::complex<double>(-0.5, 0.5), + std::complex<double>(-0., 0.5), + std::complex<double>(+0., 0.5), + std::complex<double>(0.5, 0.5), + std::complex<double>(1, 0.5), + std::complex<double>(2, 0.5), + std::complex<double>(INFINITY, 0.5), + + std::complex<double>(NAN, 1), + std::complex<double>(-INFINITY, 1), + std::complex<double>(-2, 1), + std::complex<double>(-1, 1), + std::complex<double>(-0.5, 1), + std::complex<double>(-0., 1), + std::complex<double>(+0., 1), + std::complex<double>(0.5, 1), + std::complex<double>(1, 1), + std::complex<double>(2, 1), + std::complex<double>(INFINITY, 1), + + std::complex<double>(NAN, 2), + std::complex<double>(-INFINITY, 2), + std::complex<double>(-2, 2), + std::complex<double>(-1, 2), + std::complex<double>(-0.5, 2), + std::complex<double>(-0., 2), + std::complex<double>(+0., 2), + std::complex<double>(0.5, 2), + std::complex<double>(1, 2), + std::complex<double>(2, 2), + std::complex<double>(INFINITY, 2), + + std::complex<double>(NAN, INFINITY), + std::complex<double>(-INFINITY, INFINITY), + std::complex<double>(-2, INFINITY), + std::complex<double>(-1, INFINITY), + std::complex<double>(-0.5, INFINITY), + std::complex<double>(-0., INFINITY), + std::complex<double>(+0., INFINITY), + std::complex<double>(0.5, INFINITY), + std::complex<double>(1, INFINITY), + std::complex<double>(2, INFINITY), + std::complex<double>(INFINITY, INFINITY) +}; + +enum {zero, non_zero, inf, NaN, non_zero_nan}; + +template <class T> +int +classify(const std::complex<T>& x) +{ + if (x == std::complex<T>()) + return zero; + if (std::isinf(x.real()) || std::isinf(x.imag())) + return inf; + if (std::isnan(x.real()) && std::isnan(x.imag())) + return NaN; + if (std::isnan(x.real())) + { + if (x.imag() == T(0)) + return NaN; + return non_zero_nan; + } + if (std::isnan(x.imag())) + { + if (x.real() == T(0)) + return NaN; + return non_zero_nan; + } + return non_zero; +} + +inline +int +classify(double x) +{ + if (x == 0) + return zero; + if (std::isinf(x)) + return inf; + if (std::isnan(x)) + return NaN; + return non_zero; +} + +void is_about(float x, float y) +{ + assert(std::abs((x-y)/(x+y)) < 1.e-6); +} + +void is_about(double x, double y) +{ + assert(std::abs((x-y)/(x+y)) < 1.e-14); +} + +void is_about(long double x, long double y) +{ + assert(std::abs((x-y)/(x+y)) < 1.e-14); +} + +#endif // CASES_H |
