1 /* Test for powl 2 Copyright (C) 2011-2022 Free Software Foundation, Inc. 3 This file is part of the GNU C Library. 4 5 The GNU C Library is free software; you can redistribute it and/or 6 modify it under the terms of the GNU Lesser General Public 7 License as published by the Free Software Foundation; either 8 version 2.1 of the License, or (at your option) any later version. 9 10 The GNU C Library is distributed in the hope that it will be useful, 11 but WITHOUT ANY WARRANTY; without even the implied warranty of 12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 13 Lesser General Public License for more details. 14 15 You should have received a copy of the GNU Lesser General Public 16 License along with the GNU C Library; if not, see 17 <https://www.gnu.org/licenses/>. */ 18 19 #include <stdio.h> 20 #include <math.h> 21 #include <float.h> 22 #include <ieee754.h> 23 24 static int do_test(void)25do_test (void) 26 { 27 int result = 0; 28 29 #if LDBL_MANT_DIG == 64 30 { 31 long double x = 1e-20; 32 union ieee854_long_double u; 33 u.ieee.mantissa0 = 1; 34 u.ieee.mantissa1 = 1; 35 u.ieee.exponent = 0; 36 u.ieee.negative = 0; 37 (void) powl (0.2, u.d); 38 x = powl (x, 1.5); 39 if (fabsl (x - 1e-30) > 1e-10) 40 { 41 printf ("powl (1e-20, 1.5): wrong result: %Lg\n", x); 42 result = 1; 43 } 44 } 45 #endif 46 47 return result; 48 } 49 50 #define TEST_FUNCTION do_test () 51 #include "../test-skeleton.c" 52