diff options
author | Waldemar Brodkorb <wbx@uclibc-ng.org> | 2016-11-03 13:03:55 +0100 |
---|---|---|
committer | Waldemar Brodkorb <wbx@uclibc-ng.org> | 2016-11-03 13:03:55 +0100 |
commit | 3384c45e66ddf18f235654b67ae34ac7dcb07534 (patch) | |
tree | edbd0f03b72acb5961f7919e13da34a34a84fc23 /test/math/test-powl.c | |
parent | 29d1023614fcfef8520edec8b64e9f8953d9a121 (diff) |
math: sync with GNU libc
The format of the ULPS files have changed, non-glibc architecture
files needs to be updated later.
Add all math tests from latest GNU libc and allow to compile
and run them on uClibc-ng and GNU libc systems.
Diffstat (limited to 'test/math/test-powl.c')
-rw-r--r-- | test/math/test-powl.c | 53 |
1 files changed, 53 insertions, 0 deletions
diff --git a/test/math/test-powl.c b/test/math/test-powl.c new file mode 100644 index 0000000..e67f052 --- /dev/null +++ b/test/math/test-powl.c @@ -0,0 +1,53 @@ +/* Test for powl + Copyright (C) 2011-2016 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + <http://www.gnu.org/licenses/>. */ + +#include <stdio.h> +#include <math.h> +#include <float.h> +#include <ieee754.h> + +static int +do_test (void) +{ + int result = 0; + +#ifndef NO_LONG_DOUBLE +# if LDBL_MANT_DIG == 64 + { + long double x = 1e-20; + union ieee854_long_double u; + u.ieee.mantissa0 = 1; + u.ieee.mantissa1 = 1; + u.ieee.exponent = 0; + u.ieee.negative = 0; + (void) powl (0.2, u.d); + x = powl (x, 1.5); + if (fabsl (x - 1e-30) > 1e-10) + { + printf ("powl (1e-20, 1.5): wrong result: %Lg\n", x); + result = 1; + } + } +# endif +#endif + + return result; +} + +#define TEST_FUNCTION do_test () +#include "../test-skeleton.c" |