blob: eaedcc701214c472fb21819c03f5c81f4675d53c [file] [log] [blame]
lh9ed821d2023-04-07 01:36:19 -07001/* Test for powl
2 Copyright (C) 2011-2015 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 <http://www.gnu.org/licenses/>. */
18
19#include <stdio.h>
20#include <math.h>
21#include <float.h>
22#include <ieee754.h>
23
24static int
25do_test (void)
26{
27 int result = 0;
28
29#ifndef NO_LONG_DOUBLE
30# if LDBL_MANT_DIG == 64
31 {
32 long double x = 1e-20;
33 union ieee854_long_double u;
34 u.ieee.mantissa0 = 1;
35 u.ieee.mantissa1 = 1;
36 u.ieee.exponent = 0;
37 u.ieee.negative = 0;
38 (void) powl (0.2, u.d);
39 x = powl (x, 1.5);
40 if (fabsl (x - 1e-30) > 1e-10)
41 {
42 printf ("powl (1e-20, 1.5): wrong result: %Lg\n", x);
43 result = 1;
44 }
45 }
46# endif
47#endif
48
49 return result;
50}
51
52#define TEST_FUNCTION do_test ()
53#include "../test-skeleton.c"