1 /* Test of <math.h> substitute.
2 Copyright (C) 2007-2021 Free Software Foundation, Inc.
3
4 This program is free software: you can redistribute it and/or modify
5 it under the terms of the GNU General Public License as published by
6 the Free Software Foundation; either version 3 of the License, or
7 (at your option) any later version.
8
9 This program is distributed in the hope that it will be useful,
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 GNU General Public License for more details.
13
14 You should have received a copy of the GNU General Public License
15 along with this program. If not, see <https://www.gnu.org/licenses/>. */
16
17 /* Written by Bruno Haible <bruno@clisp.org>, 2007. */
18
19 #include <config.h>
20
21 #include <math.h>
22
23 #ifndef NAN
24 # error NAN should be defined
25 choke me
26 #endif
27
28 #ifndef HUGE_VALF
29 # error HUGE_VALF should be defined
30 choke me
31 #endif
32
33 #ifndef HUGE_VAL
34 # error HUGE_VAL should be defined
35 choke me
36 #endif
37
38 #ifndef HUGE_VALL
39 # error HUGE_VALL should be defined
40 choke me
41 #endif
42
43 #ifndef FP_ILOGB0
44 # error FP_ILOGB0 should be defined
45 choke me
46 #endif
47
48 #ifndef FP_ILOGBNAN
49 # error FP_ILOGBNAN should be defined
50 choke me
51 #endif
52
53 #include <limits.h>
54
55 #include "macros.h"
56
57 #if 0
58 /* Check that NAN expands into a constant expression. */
59 static float n = NAN;
60 #endif
61
62 /* Compare two numbers with ==.
63 This is a separate function because IRIX 6.5 "cc -O" miscompiles an
64 'x == x' test. */
65 static int
numeric_equalf(float x,float y)66 numeric_equalf (float x, float y)
67 {
68 return x == y;
69 }
70 static int
numeric_equald(double x,double y)71 numeric_equald (double x, double y)
72 {
73 return x == y;
74 }
75 static int
numeric_equall(long double x,long double y)76 numeric_equall (long double x, long double y)
77 {
78 return x == y;
79 }
80
81 int
main(void)82 main (void)
83 {
84 double d = NAN;
85 double zero = 0.0;
86 ASSERT (!numeric_equald (d, d));
87
88 d = HUGE_VAL;
89 ASSERT (numeric_equald (d, 1.0 / zero));
90
91 ASSERT (numeric_equalf (HUGE_VALF, HUGE_VALF + HUGE_VALF));
92
93 ASSERT (numeric_equald (HUGE_VAL, HUGE_VAL + HUGE_VAL));
94
95 ASSERT (numeric_equall (HUGE_VALL, HUGE_VALL + HUGE_VALL));
96
97 /* Check the value of FP_ILOGB0. */
98 ASSERT (FP_ILOGB0 == INT_MIN || FP_ILOGB0 == - INT_MAX);
99
100 /* Check the value of FP_ILOGBNAN. */
101 ASSERT (FP_ILOGBNAN == INT_MIN || FP_ILOGBNAN == INT_MAX);
102
103 return 0;
104 }
105