133 lines
4.1 KiB
C
133 lines
4.1 KiB
C
/* 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 Library General Public
|
|
* License as published by the Free Software Foundation; either
|
|
* version 2 of the License, or (at your option) any later version.
|
|
*
|
|
* This 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
|
|
* Library General Public License for more details.
|
|
*
|
|
* You should have received a copy of the GNU Library General Public
|
|
* License along with this library; if not, write to the Free
|
|
* Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, 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 */
|
|
|
|
extern __fpmax_t __strtofpmax(const char *str, char **endptr, int exp_adjust);
|
|
|
|
#ifdef __UCLIBC_HAS_XLOCALE__
|
|
extern __fpmax_t __strtofpmax_l(const char *str, char **endptr, int exp_adjust,
|
|
__locale_t locale_arg);
|
|
#endif
|
|
|
|
#ifdef __UCLIBC_HAS_WCHAR__
|
|
extern __fpmax_t __wcstofpmax(const wchar_t *wcs, wchar_t **endptr,
|
|
int exp_adjust);
|
|
|
|
#ifdef __UCLIBC_HAS_XLOCALE__
|
|
extern __fpmax_t __wcstofpmax_l(const wchar_t *wcs, wchar_t **endptr,
|
|
int exp_adjust, __locale_t locale_arg);
|
|
#endif
|
|
#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 */
|