diff options
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" |