mirror of
https://sourceware.org/git/glibc.git
synced 2025-01-11 11:50:06 +00:00
ec94343f59
TS 18661-1 defines a type femode_t to represent the set of dynamic floating-point control modes (such as the rounding mode and trap enablement modes), and functions fegetmode and fesetmode to manipulate those modes (without affecting other state such as the raised exception flags) and a corresponding macro FE_DFL_MODE. This patch series implements those interfaces for glibc. This first patch adds the architecture-independent pieces, the x86 and x86_64 implementations, and the <bits/fenv.h> and ABI baseline updates for all architectures so glibc keeps building and passing the ABI tests on all architectures. Subsequent patches add the fegetmode and fesetmode implementations for other architectures. femode_t is generally an integer type - the same type as fenv_t, or as the single element of fenv_t where fenv_t is a structure containing a single integer (or the single relevant element, where it has elements for both status and control registers) - except where architecture properties or consistency with the fenv_t implementation indicate otherwise. FE_DFL_MODE follows FE_DFL_ENV in whether it's a magic pointer value (-1 cast to const femode_t *), a value that can be distinguished from valid pointers by its high bits but otherwise contains a representation of the desired register contents, or a pointer to a constant variable (the powerpc case; __fe_dfl_mode is added as an exported constant object, an alias to __fe_dfl_env). Note that where architectures (that share a register between control and status bits) gain definitions of new floating-point control or status bits in future, the implementations of fesetmode for those architectures may need updating (depending on whether the new bits are control or status bits and what the implementation does with previously unknown bits), just like existing implementations of <fenv.h> functions that take care not to touch reserved bits may need updating when the set of reserved bits changes. (As any new bits are outside the scope of ISO C, that's just a quality-of-implementation issue for supporting them, not a conformance issue.) As with fenv_t, femode_t should properly include any software DFP rounding mode (and for both fenv_t and femode_t I'd consider that fragment of DFP support appropriate for inclusion in glibc even in the absence of the rest of libdfp; hardware DFP rounding modes should already be included if the definitions of which bits are status / control bits are correct). Tested for x86_64, x86, mips64 (hard float, and soft float to test the fallback version), arm (hard float) and powerpc (hard float, soft float and e500). Other architecture versions are untested. * math/fegetmode.c: New file. * math/fesetmode.c: Likewise. * sysdeps/i386/fpu/fegetmode.c: Likewise. * sysdeps/i386/fpu/fesetmode.c: Likewise. * sysdeps/x86_64/fpu/fegetmode.c: Likewise. * sysdeps/x86_64/fpu/fesetmode.c: Likewise. * math/fenv.h: Update comment on inclusion of <bits/fenv.h>. [__GLIBC_USE (IEC_60559_BFP_EXT)] (fegetmode): New function declaration. [__GLIBC_USE (IEC_60559_BFP_EXT)] (fesetmode): Likewise. * bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/aarch64/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/alpha/fpu/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/arm/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/hppa/fpu/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/ia64/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/m68k/fpu/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/microblaze/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/mips/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/nios2/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/powerpc/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (__fe_dfl_mode): New variable declaration. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/s390/fpu/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/sh/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/sparc/fpu/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/tile/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * sysdeps/x86/fpu/bits/fenv.h [__GLIBC_USE (IEC_60559_BFP_EXT)] (femode_t): New typedef. [__GLIBC_USE (IEC_60559_BFP_EXT)] (FE_DFL_MODE): New macro. * manual/arith.texi (FE_DFL_MODE): Document macro. (fegetmode): Document function. (fesetmode): Likewise. * math/Versions (fegetmode): New libm symbol at version GLIBC_2.25. (fesetmode): Likewise. * math/Makefile (libm-support): Add fegetmode and fesetmode. (tests): Add test-femode and test-femode-traps. * math/test-femode-traps.c: New file. * math/test-femode.c: Likewise. * sysdeps/powerpc/fpu/fenv_const.c (__fe_dfl_mode): Declare as alias for __fe_dfl_env. * sysdeps/powerpc/nofpu/fenv_const.c (__fe_dfl_mode): Likewise. * sysdeps/powerpc/powerpc32/e500/nofpu/fenv_const.c (__fe_dfl_mode): Likewise. * sysdeps/powerpc/Versions (__fe_dfl_mode): New libm symbol at version GLIBC_2.25. * sysdeps/nacl/libm.abilist: Update. * sysdeps/unix/sysv/linux/aarch64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/alpha/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/arm/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/hppa/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/i386/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/ia64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/m68k/coldfire/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/m68k/m680x0/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/microblaze/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/mips/mips32/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/mips/mips64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/nios2/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/powerpc/powerpc32/fpu/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/powerpc/powerpc32/nofpu/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/powerpc/powerpc64/libm-le.abilist: Likewise. * sysdeps/unix/sysv/linux/powerpc/powerpc64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/s390/s390-32/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/s390/s390-64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/sh/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/sparc/sparc32/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/sparc/sparc64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/tile/tilegx/tilegx32/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/tile/tilegx/tilegx64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/tile/tilepro/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/x86_64/64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/x86_64/x32/libm.abilist: Likewise.
181 lines
6.0 KiB
C
181 lines
6.0 KiB
C
/* Copyright (C) 1997-2016 Free Software Foundation, Inc.
|
|
This file is part of the GNU C Library.
|
|
|
|
The GNU C 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, see
|
|
<http://www.gnu.org/licenses/>. */
|
|
|
|
#ifndef _FENV_H
|
|
# error "Never use <bits/fenv.h> directly; include <fenv.h> instead."
|
|
#endif
|
|
|
|
|
|
/* Define bits representing the exception. We use the bit positions of
|
|
the appropriate bits in the FPSCR... */
|
|
enum
|
|
{
|
|
FE_INEXACT =
|
|
#define FE_INEXACT (1 << (31 - 6))
|
|
FE_INEXACT,
|
|
FE_DIVBYZERO =
|
|
#define FE_DIVBYZERO (1 << (31 - 5))
|
|
FE_DIVBYZERO,
|
|
FE_UNDERFLOW =
|
|
#define FE_UNDERFLOW (1 << (31 - 4))
|
|
FE_UNDERFLOW,
|
|
FE_OVERFLOW =
|
|
#define FE_OVERFLOW (1 << (31 - 3))
|
|
FE_OVERFLOW,
|
|
|
|
/* ... except for FE_INVALID, for which we use bit 31. FE_INVALID
|
|
actually corresponds to bits 7 through 12 and 21 through 23
|
|
in the FPSCR, but we can't use that because the current draft
|
|
says that it must be a power of 2. Instead we use bit 2 which
|
|
is the summary bit for all the FE_INVALID exceptions, which
|
|
kind of makes sense. */
|
|
FE_INVALID =
|
|
#define FE_INVALID (1 << (31 - 2))
|
|
FE_INVALID,
|
|
|
|
#ifdef __USE_GNU
|
|
/* Breakdown of the FE_INVALID bits. Setting FE_INVALID on an
|
|
input to a routine is equivalent to setting all of these bits;
|
|
FE_INVALID will be set on output from a routine iff one of
|
|
these bits is set. Note, though, that you can't disable or
|
|
enable these exceptions individually. */
|
|
|
|
/* Operation with a sNaN. */
|
|
FE_INVALID_SNAN =
|
|
# define FE_INVALID_SNAN (1 << (31 - 7))
|
|
FE_INVALID_SNAN,
|
|
|
|
/* Inf - Inf */
|
|
FE_INVALID_ISI =
|
|
# define FE_INVALID_ISI (1 << (31 - 8))
|
|
FE_INVALID_ISI,
|
|
|
|
/* Inf / Inf */
|
|
FE_INVALID_IDI =
|
|
# define FE_INVALID_IDI (1 << (31 - 9))
|
|
FE_INVALID_IDI,
|
|
|
|
/* 0 / 0 */
|
|
FE_INVALID_ZDZ =
|
|
# define FE_INVALID_ZDZ (1 << (31 - 10))
|
|
FE_INVALID_ZDZ,
|
|
|
|
/* Inf * 0 */
|
|
FE_INVALID_IMZ =
|
|
# define FE_INVALID_IMZ (1 << (31 - 11))
|
|
FE_INVALID_IMZ,
|
|
|
|
/* Comparison with a NaN. */
|
|
FE_INVALID_COMPARE =
|
|
# define FE_INVALID_COMPARE (1 << (31 - 12))
|
|
FE_INVALID_COMPARE,
|
|
|
|
/* Invalid operation flag for software (not set by hardware). */
|
|
/* Note that some chips don't have this implemented, presumably
|
|
because no-one expected anyone to write software for them %-). */
|
|
FE_INVALID_SOFTWARE =
|
|
# define FE_INVALID_SOFTWARE (1 << (31 - 21))
|
|
FE_INVALID_SOFTWARE,
|
|
|
|
/* Square root of negative number (including -Inf). */
|
|
/* Note that some chips don't have this implemented. */
|
|
FE_INVALID_SQRT =
|
|
# define FE_INVALID_SQRT (1 << (31 - 22))
|
|
FE_INVALID_SQRT,
|
|
|
|
/* Conversion-to-integer of a NaN or a number too large or too small. */
|
|
FE_INVALID_INTEGER_CONVERSION =
|
|
# define FE_INVALID_INTEGER_CONVERSION (1 << (31 - 23))
|
|
FE_INVALID_INTEGER_CONVERSION
|
|
|
|
# define FE_ALL_INVALID \
|
|
(FE_INVALID_SNAN | FE_INVALID_ISI | FE_INVALID_IDI | FE_INVALID_ZDZ \
|
|
| FE_INVALID_IMZ | FE_INVALID_COMPARE | FE_INVALID_SOFTWARE \
|
|
| FE_INVALID_SQRT | FE_INVALID_INTEGER_CONVERSION)
|
|
#endif
|
|
};
|
|
|
|
#define FE_ALL_EXCEPT \
|
|
(FE_INEXACT | FE_DIVBYZERO | FE_UNDERFLOW | FE_OVERFLOW | FE_INVALID)
|
|
|
|
/* PowerPC chips support all of the four defined rounding modes. We
|
|
use the bit pattern in the FPSCR as the values for the
|
|
appropriate macros. */
|
|
enum
|
|
{
|
|
FE_TONEAREST =
|
|
#define FE_TONEAREST 0
|
|
FE_TONEAREST,
|
|
FE_TOWARDZERO =
|
|
#define FE_TOWARDZERO 1
|
|
FE_TOWARDZERO,
|
|
FE_UPWARD =
|
|
#define FE_UPWARD 2
|
|
FE_UPWARD,
|
|
FE_DOWNWARD =
|
|
#define FE_DOWNWARD 3
|
|
FE_DOWNWARD
|
|
};
|
|
|
|
/* Type representing exception flags. */
|
|
typedef unsigned int fexcept_t;
|
|
|
|
/* Type representing floating-point environment. We leave it as 'double'
|
|
for efficiency reasons (rather than writing it to a 32-bit integer). */
|
|
typedef double fenv_t;
|
|
|
|
/* If the default argument is used we use this value. */
|
|
extern const fenv_t __fe_dfl_env;
|
|
#define FE_DFL_ENV (&__fe_dfl_env)
|
|
|
|
#ifdef __USE_GNU
|
|
/* Floating-point environment where all exceptions are enabled. Note that
|
|
this is not sufficient to give you SIGFPE. */
|
|
extern const fenv_t __fe_enabled_env;
|
|
# define FE_ENABLED_ENV (&__fe_enabled_env)
|
|
|
|
/* Floating-point environment with (processor-dependent) non-IEEE floating
|
|
point. */
|
|
extern const fenv_t __fe_nonieee_env;
|
|
# define FE_NONIEEE_ENV (&__fe_nonieee_env)
|
|
|
|
/* Floating-point environment with all exceptions enabled. Note that
|
|
just evaluating this value does not change the processor exception mode.
|
|
Passing this mask to fesetenv will result in a prctl syscall to change
|
|
the MSR FE0/FE1 bits to "Precise Mode". On some processors this will
|
|
result in slower floating point execution. This will last until an
|
|
fenv or exception mask is installed that disables all FP exceptions. */
|
|
# define FE_NOMASK_ENV FE_ENABLED_ENV
|
|
|
|
/* Floating-point environment with all exceptions disabled. Note that
|
|
just evaluating this value does not change the processor exception mode.
|
|
Passing this mask to fesetenv will result in a prctl syscall to change
|
|
the MSR FE0/FE1 bits to "Ignore Exceptions Mode". On most processors
|
|
this allows the fastest possible floating point execution.*/
|
|
# define FE_MASK_ENV FE_DFL_ENV
|
|
|
|
#endif
|
|
|
|
#if __GLIBC_USE (IEC_60559_BFP_EXT)
|
|
/* Type representing floating-point control modes. */
|
|
typedef double femode_t;
|
|
|
|
/* Default floating-point control modes. */
|
|
extern const femode_t __fe_dfl_mode;
|
|
# define FE_DFL_MODE (&__fe_dfl_mode)
|
|
#endif
|