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
24static int
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 (x: 0.2, y: u.d);
38 x = powl (x: x, y: 1.5);
39 if (fabsl (x: x - 1e-30) > 1e-10)
40 {
41 printf (format: "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

source code of glibc/math/test-powl.c