1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
|
/* Copyright (C) 2003 Manuel Novoa III
*
* This 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, write to the Free
* Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
* 02111-1307 USA.
*/
/* Define a maximal floating point type, and the associated constants
* that are defined for the floating point types in float.h.
*
* This is to support archs that are missing long double, or even double.
*/
#ifndef _UCLIBC_FPMAX_H
#define _UCLIBC_FPMAX_H
#ifndef _ISOC99_SOURCE
#define _ISOC99_SOURCE 1
#endif
#include <features.h>
#include <float.h>
#ifdef __UCLIBC_HAS_FLOATS__
#if defined(LDBL_MANT_DIG)
typedef long double __fpmax_t;
#define FPMAX_TYPE 3
#define FPMAX_MANT_DIG LDBL_MANT_DIG
#define FPMAX_DIG LDBL_DIG
#define FPMAX_EPSILON LDBL_EPSILON
#define FPMAX_MIN_EXP LDBL_MIN_EXP
#define FPMAX_MIN LDBL_MIN
#define FPMAX_MIN_10_EXP LDBL_MIN_10_EXP
#define FPMAX_MAX_EXP LDBL_MAX_EXP
#define FPMAX_MAX LDBL_MAX
#define FPMAX_MAX_10_EXP LDBL_MAX_10_EXP
#elif defined(DBL_MANT_DIG)
typedef double __fpmax_t;
#define FPMAX_TYPE 2
#define FPMAX_MANT_DIG DBL_MANT_DIG
#define FPMAX_DIG DBL_DIG
#define FPMAX_EPSILON DBL_EPSILON
#define FPMAX_MIN_EXP DBL_MIN_EXP
#define FPMAX_MIN DBL_MIN
#define FPMAX_MIN_10_EXP DBL_MIN_10_EXP
#define FPMAX_MAX_EXP DBL_MAX_EXP
#define FPMAX_MAX DBL_MAX
#define FPMAX_MAX_10_EXP DBL_MAX_10_EXP
#elif defined(FLT_MANT_DIG)
typedef float __fpmax_t;
#define FPMAX_TYPE 1
#define FPMAX_MANT_DIG FLT_MANT_DIG
#define FPMAX_DIG FLT_DIG
#define FPMAX_EPSILON FLT_EPSILON
#define FPMAX_MIN_EXP FLT_MIN_EXP
#define FPMAX_MIN FLT_MIN
#define FPMAX_MIN_10_EXP FLT_MIN_10_EXP
#define FPMAX_MAX_EXP FLT_MAX_EXP
#define FPMAX_MAX FLT_MAX
#define FPMAX_MAX_10_EXP FLT_MAX_10_EXP
#else
#error unable to determine appropriate type for __fpmax_t!
#endif
#ifndef DECIMAL_DIG
#ifdef L___strtofpmax
/* Emit warning only once. */
#warning DECIMAL_DIG is not defined! If you are using gcc, it may not be defining __STDC_VERSION__ as it should.
#endif
#if !defined(FLT_RADIX) || (FLT_RADIX != 2)
#error unable to compensate for missing DECIMAL_DIG!
#endif
/* ceil (1 + #mantissa * log10 (FLT_RADIX)) */
#define DECIMAL_DIG (1 + (((FPMAX_MANT_DIG * 100) + 331) / 332))
#endif /* DECIMAL_DIG */
#if defined _LIBC && defined IS_IN_libc
extern __fpmax_t __strtofpmax(const char *str, char **endptr, int exp_adjust) attribute_hidden;
#ifdef __UCLIBC_HAS_XLOCALE__
extern __fpmax_t __strtofpmax_l(const char *str, char **endptr, int exp_adjust,
__locale_t locale_arg) attribute_hidden;
#endif
#ifdef __UCLIBC_HAS_WCHAR__
extern __fpmax_t __wcstofpmax(const wchar_t *wcs, wchar_t **endptr,
int exp_adjust) attribute_hidden;
#ifdef __UCLIBC_HAS_XLOCALE__
extern __fpmax_t __wcstofpmax_l(const wchar_t *wcs, wchar_t **endptr,
int exp_adjust, __locale_t locale_arg) attribute_hidden;
#endif
#endif /* _LIBC */
#endif /* __UCLIBC_HAS_WCHAR__ */
/* The following checks in an __fpmax_t is either 0 or +/- infinity.
*
* WARNING!!! WARNING!!! WARNING!!! WARNING!!! WARNING!!! WARNING!!!
*
* This only works if __fpmax_t is the actual maximal floating point type used
* in intermediate calculations. Otherwise, excess precision in the
* intermediate values can cause the test to fail.
*
* WARNING!!! WARNING!!! WARNING!!! WARNING!!! WARNING!!! WARNING!!!
*/
#define __FPMAX_ZERO_OR_INF_CHECK(x) ((x) == ((x)/4) )
#endif /* __UCLIBC_HAS_FLOATS__ */
#endif /* _UCLIBC_FPMAX_H */
|
