2013-04-02 11:51:02 +00:00
|
|
|
/* Internal math stuff. MIPS version.
|
2015-01-02 16:28:19 +00:00
|
|
|
Copyright (C) 2013-2015 Free Software Foundation, Inc.
|
2013-04-02 11:51:02 +00:00
|
|
|
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 _MATH_PRIVATE_H
|
|
|
|
|
MIPS: IEEE 754-2008 NaN encoding support
It has been a long practice for software using IEEE 754 floating-point
arithmetic run on MIPS processors to use an encoding of Not-a-Number
(NaN) data different to one used by software run on other processors.
And as of IEEE 754-2008 revision [1] this encoding does not follow one
recommended in the standard, as specified in section 6.2.1, where it
is stated that quiet NaNs should have the first bit (d1) of their
significand set to 1 while signalling NaNs should have that bit set to
0, but MIPS software interprets the two bits in the opposite manner.
As from revision 3.50 [2][3] the MIPS Architecture provides for
processors that support the IEEE 754-2008 preferred NaN encoding format.
As the two formats (further referred to as "legacy NaN" and "2008 NaN")
are incompatible to each other, tools have to provide support for the
two formats to help people avoid using incompatible binary modules.
The change is comprised of two functional groups of features, both of
which are required for correct support.
1. Dynamic linker support.
To enforce the NaN encoding requirement in dynamic linking a new ELF
file header flag has been defined. This flag is set for 2008-NaN
shared modules and executables and clear for legacy-NaN ones. The
dynamic linker silently ignores any incompatible modules it
encounters in dependency processing.
To avoid unnecessary processing of incompatible modules in the
presence of a shared module cache, a set of new cache flags has been
defined to mark 2008-NaN modules for the three ABIs supported.
Changes to sysdeps/unix/sysv/linux/mips/readelflib.c have been made
following an earlier code quality suggestion made here:
http://sourceware.org/ml/libc-ports/2009-03/msg00036.html
and are therefore a little bit more extensive than the minimum
required.
Finally a new name has been defined for the dynamic linker so that
2008-NaN and legacy-NaN binaries can coexist on a single system that
supports dual-mode operation and that a legacy dynamic linker that
does not support verifying the 2008-NaN ELF file header flag is not
chosen to interpret a 2008-NaN binary by accident.
2. Floating environment support.
IEEE 754-2008 features are controlled in the Floating-Point Control
and Status (FCSR) register and updates are needed to floating
environment support so that the 2008-NaN flag is set correctly and
the kernel default, inferred from the 2008-NaN ELF file header flag
at the time an executable is loaded, respected.
As the NaN encoding format is a property of GCC code generation that is
both a user-selected GCC configuration default and can be overridden
with GCC options, code that needs to know what NaN encoding standard it
has been configured for checks for the __mips_nan2008 macro that is
defined internally by GCC whenever the 2008-NaN mode has been selected.
This mode is determined at the glibc configuration time and therefore a
few consistency checks have been added to catch cases where compilation
flags have been overridden by the user.
The 2008 NaN set of features relies on kernel support as the in-kernel
floating-point emulator needs to be aware of the NaN encoding used even
on hard-float processors and configure the FPU context according to the
value of the 2008 NaN ELF file header flag of the executable being
started. As at this time work on kernel support is still in progress
and the relevant changes have not made their way yet to linux.org master
repository.
Therefore the minimum version supported has been artificially set to
10.0.0 so that 2008-NaN code is not accidentally run on a Linux kernel
that does not suppport it. It is anticipated that the version is
adjusted later on to the actual initial linux.org kernel version to
support this feature. Legacy NaN encoding support is unaffected, older
kernel versions remain supported.
[1] "IEEE Standard for Floating-Point Arithmetic", IEEE Computer
Society, IEEE Std 754-2008, 29 August 2008
[2] "MIPS Architecture For Programmers, Volume I-A: Introduction to the
MIPS32 Architecture", MIPS Technologies, Inc., Document Number:
MD00082, Revision 3.50, September 20, 2012
[3] "MIPS Architecture For Programmers, Volume I-A: Introduction to the
MIPS64 Architecture", MIPS Technologies, Inc., Document Number:
MD00083, Revision 3.50, September 20, 2012
2013-09-18 20:04:27 +00:00
|
|
|
#ifdef __mips_nan2008
|
|
|
|
/* MIPS aligned to IEEE 754-2008. */
|
|
|
|
#else
|
2013-04-02 11:51:02 +00:00
|
|
|
/* One of the few architectures where the meaning of the quiet/signaling bit is
|
|
|
|
inverse to IEEE 754-2008 (as well as common practice for IEEE 754-1985). */
|
MIPS: IEEE 754-2008 NaN encoding support
It has been a long practice for software using IEEE 754 floating-point
arithmetic run on MIPS processors to use an encoding of Not-a-Number
(NaN) data different to one used by software run on other processors.
And as of IEEE 754-2008 revision [1] this encoding does not follow one
recommended in the standard, as specified in section 6.2.1, where it
is stated that quiet NaNs should have the first bit (d1) of their
significand set to 1 while signalling NaNs should have that bit set to
0, but MIPS software interprets the two bits in the opposite manner.
As from revision 3.50 [2][3] the MIPS Architecture provides for
processors that support the IEEE 754-2008 preferred NaN encoding format.
As the two formats (further referred to as "legacy NaN" and "2008 NaN")
are incompatible to each other, tools have to provide support for the
two formats to help people avoid using incompatible binary modules.
The change is comprised of two functional groups of features, both of
which are required for correct support.
1. Dynamic linker support.
To enforce the NaN encoding requirement in dynamic linking a new ELF
file header flag has been defined. This flag is set for 2008-NaN
shared modules and executables and clear for legacy-NaN ones. The
dynamic linker silently ignores any incompatible modules it
encounters in dependency processing.
To avoid unnecessary processing of incompatible modules in the
presence of a shared module cache, a set of new cache flags has been
defined to mark 2008-NaN modules for the three ABIs supported.
Changes to sysdeps/unix/sysv/linux/mips/readelflib.c have been made
following an earlier code quality suggestion made here:
http://sourceware.org/ml/libc-ports/2009-03/msg00036.html
and are therefore a little bit more extensive than the minimum
required.
Finally a new name has been defined for the dynamic linker so that
2008-NaN and legacy-NaN binaries can coexist on a single system that
supports dual-mode operation and that a legacy dynamic linker that
does not support verifying the 2008-NaN ELF file header flag is not
chosen to interpret a 2008-NaN binary by accident.
2. Floating environment support.
IEEE 754-2008 features are controlled in the Floating-Point Control
and Status (FCSR) register and updates are needed to floating
environment support so that the 2008-NaN flag is set correctly and
the kernel default, inferred from the 2008-NaN ELF file header flag
at the time an executable is loaded, respected.
As the NaN encoding format is a property of GCC code generation that is
both a user-selected GCC configuration default and can be overridden
with GCC options, code that needs to know what NaN encoding standard it
has been configured for checks for the __mips_nan2008 macro that is
defined internally by GCC whenever the 2008-NaN mode has been selected.
This mode is determined at the glibc configuration time and therefore a
few consistency checks have been added to catch cases where compilation
flags have been overridden by the user.
The 2008 NaN set of features relies on kernel support as the in-kernel
floating-point emulator needs to be aware of the NaN encoding used even
on hard-float processors and configure the FPU context according to the
value of the 2008 NaN ELF file header flag of the executable being
started. As at this time work on kernel support is still in progress
and the relevant changes have not made their way yet to linux.org master
repository.
Therefore the minimum version supported has been artificially set to
10.0.0 so that 2008-NaN code is not accidentally run on a Linux kernel
that does not suppport it. It is anticipated that the version is
adjusted later on to the actual initial linux.org kernel version to
support this feature. Legacy NaN encoding support is unaffected, older
kernel versions remain supported.
[1] "IEEE Standard for Floating-Point Arithmetic", IEEE Computer
Society, IEEE Std 754-2008, 29 August 2008
[2] "MIPS Architecture For Programmers, Volume I-A: Introduction to the
MIPS32 Architecture", MIPS Technologies, Inc., Document Number:
MD00082, Revision 3.50, September 20, 2012
[3] "MIPS Architecture For Programmers, Volume I-A: Introduction to the
MIPS64 Architecture", MIPS Technologies, Inc., Document Number:
MD00083, Revision 3.50, September 20, 2012
2013-09-18 20:04:27 +00:00
|
|
|
# define HIGH_ORDER_BIT_IS_SET_FOR_SNAN
|
|
|
|
#endif
|
2013-04-02 11:51:02 +00:00
|
|
|
|
2013-09-23 16:58:30 +00:00
|
|
|
/* Inline functions to speed up the math library implementation. The
|
|
|
|
default versions of these routines are in generic/math_private.h
|
|
|
|
and call fesetround, feholdexcept, etc. These routines use inlined
|
|
|
|
code instead. */
|
|
|
|
|
|
|
|
#ifdef __mips_hard_float
|
|
|
|
|
|
|
|
# include <fenv.h>
|
|
|
|
# include <fenv_libc.h>
|
|
|
|
# include <fpu_control.h>
|
|
|
|
|
2014-02-26 22:19:20 +00:00
|
|
|
# define _FPU_MASK_ALL (_FPU_MASK_V | _FPU_MASK_Z | _FPU_MASK_O \
|
|
|
|
|_FPU_MASK_U | _FPU_MASK_I | FE_ALL_EXCEPT)
|
|
|
|
|
2013-09-23 16:58:30 +00:00
|
|
|
static __always_inline void
|
|
|
|
libc_feholdexcept_mips (fenv_t *envp)
|
|
|
|
{
|
|
|
|
fpu_control_t cw;
|
|
|
|
|
|
|
|
/* Save the current state. */
|
|
|
|
_FPU_GETCW (cw);
|
|
|
|
envp->__fp_control_register = cw;
|
|
|
|
|
|
|
|
/* Clear all exception enable bits and flags. */
|
2014-02-26 22:19:20 +00:00
|
|
|
cw &= ~(_FPU_MASK_ALL);
|
2013-09-23 16:58:30 +00:00
|
|
|
_FPU_SETCW (cw);
|
|
|
|
}
|
|
|
|
# define libc_feholdexcept libc_feholdexcept_mips
|
|
|
|
# define libc_feholdexceptf libc_feholdexcept_mips
|
|
|
|
# define libc_feholdexceptl libc_feholdexcept_mips
|
|
|
|
|
|
|
|
static __always_inline void
|
|
|
|
libc_fesetround_mips (int round)
|
|
|
|
{
|
|
|
|
fpu_control_t cw;
|
|
|
|
|
|
|
|
/* Get current state. */
|
|
|
|
_FPU_GETCW (cw);
|
|
|
|
|
|
|
|
/* Set rounding bits. */
|
|
|
|
cw &= ~_FPU_RC_MASK;
|
|
|
|
cw |= round;
|
|
|
|
|
|
|
|
/* Set new state. */
|
|
|
|
_FPU_SETCW (cw);
|
|
|
|
}
|
|
|
|
# define libc_fesetround libc_fesetround_mips
|
|
|
|
# define libc_fesetroundf libc_fesetround_mips
|
|
|
|
# define libc_fesetroundl libc_fesetround_mips
|
|
|
|
|
|
|
|
static __always_inline void
|
|
|
|
libc_feholdexcept_setround_mips (fenv_t *envp, int round)
|
|
|
|
{
|
|
|
|
fpu_control_t cw;
|
|
|
|
|
|
|
|
/* Save the current state. */
|
|
|
|
_FPU_GETCW (cw);
|
|
|
|
envp->__fp_control_register = cw;
|
|
|
|
|
|
|
|
/* Clear all exception enable bits and flags. */
|
2014-02-26 22:19:20 +00:00
|
|
|
cw &= ~(_FPU_MASK_ALL);
|
2013-09-23 16:58:30 +00:00
|
|
|
|
|
|
|
/* Set rounding bits. */
|
|
|
|
cw &= ~_FPU_RC_MASK;
|
|
|
|
cw |= round;
|
|
|
|
|
|
|
|
/* Set new state. */
|
|
|
|
_FPU_SETCW (cw);
|
|
|
|
}
|
|
|
|
# define libc_feholdexcept_setround libc_feholdexcept_setround_mips
|
|
|
|
# define libc_feholdexcept_setroundf libc_feholdexcept_setround_mips
|
|
|
|
# define libc_feholdexcept_setroundl libc_feholdexcept_setround_mips
|
|
|
|
|
2014-02-26 22:19:20 +00:00
|
|
|
# define libc_feholdsetround libc_feholdexcept_setround_mips
|
|
|
|
# define libc_feholdsetroundf libc_feholdexcept_setround_mips
|
|
|
|
# define libc_feholdsetroundl libc_feholdexcept_setround_mips
|
|
|
|
|
2013-09-23 16:58:30 +00:00
|
|
|
static __always_inline void
|
|
|
|
libc_fesetenv_mips (fenv_t *envp)
|
|
|
|
{
|
2015-09-24 23:12:49 +00:00
|
|
|
fpu_control_t cw __attribute__ ((unused));
|
2013-09-23 16:58:30 +00:00
|
|
|
|
|
|
|
/* Read current state to flush fpu pipeline. */
|
|
|
|
_FPU_GETCW (cw);
|
|
|
|
|
|
|
|
_FPU_SETCW (envp->__fp_control_register);
|
|
|
|
}
|
|
|
|
# define libc_fesetenv libc_fesetenv_mips
|
|
|
|
# define libc_fesetenvf libc_fesetenv_mips
|
|
|
|
# define libc_fesetenvl libc_fesetenv_mips
|
|
|
|
|
2014-02-26 22:19:20 +00:00
|
|
|
static __always_inline int
|
|
|
|
libc_feupdateenv_test_mips (fenv_t *envp, int excepts)
|
2013-09-23 16:58:30 +00:00
|
|
|
{
|
2014-02-26 22:19:20 +00:00
|
|
|
/* int ret = fetestexcept (excepts); feupdateenv (envp); return ret; */
|
|
|
|
int cw, temp;
|
2013-09-23 16:58:30 +00:00
|
|
|
|
2014-02-26 22:19:20 +00:00
|
|
|
/* Get current control word. */
|
|
|
|
_FPU_GETCW (cw);
|
2013-09-23 16:58:30 +00:00
|
|
|
|
|
|
|
/* Set flag bits (which are accumulative), and *also* set the
|
|
|
|
cause bits. The setting of the cause bits is what actually causes
|
|
|
|
the hardware to generate the exception, if the corresponding enable
|
|
|
|
bit is set as well. */
|
2014-02-26 22:19:20 +00:00
|
|
|
temp = cw & FE_ALL_EXCEPT;
|
2013-09-23 16:58:30 +00:00
|
|
|
temp |= envp->__fp_control_register | (temp << CAUSE_SHIFT);
|
|
|
|
|
|
|
|
/* Set new state. */
|
|
|
|
_FPU_SETCW (temp);
|
2014-02-26 22:19:20 +00:00
|
|
|
|
|
|
|
return cw & excepts & FE_ALL_EXCEPT;
|
|
|
|
}
|
|
|
|
# define libc_feupdateenv_test libc_feupdateenv_test_mips
|
|
|
|
# define libc_feupdateenv_testf libc_feupdateenv_test_mips
|
|
|
|
# define libc_feupdateenv_testl libc_feupdateenv_test_mips
|
|
|
|
|
|
|
|
static __always_inline void
|
|
|
|
libc_feupdateenv_mips (fenv_t *envp)
|
|
|
|
{
|
|
|
|
libc_feupdateenv_test_mips (envp, 0);
|
2013-09-23 16:58:30 +00:00
|
|
|
}
|
|
|
|
# define libc_feupdateenv libc_feupdateenv_mips
|
|
|
|
# define libc_feupdateenvf libc_feupdateenv_mips
|
|
|
|
# define libc_feupdateenvl libc_feupdateenv_mips
|
|
|
|
|
2014-02-26 22:19:20 +00:00
|
|
|
# define libc_feresetround libc_feupdateenv_mips
|
|
|
|
# define libc_feresetroundf libc_feupdateenv_mips
|
|
|
|
# define libc_feresetroundl libc_feupdateenv_mips
|
|
|
|
|
|
|
|
static __always_inline int
|
|
|
|
libc_fetestexcept_mips (int excepts)
|
|
|
|
{
|
|
|
|
int cw;
|
|
|
|
|
|
|
|
/* Get current control word. */
|
|
|
|
_FPU_GETCW (cw);
|
|
|
|
|
|
|
|
return cw & excepts & FE_ALL_EXCEPT;
|
|
|
|
}
|
|
|
|
# define libc_fetestexcept libc_fetestexcept_mips
|
|
|
|
# define libc_fetestexceptf libc_fetestexcept_mips
|
|
|
|
# define libc_fetestexceptl libc_fetestexcept_mips
|
|
|
|
|
|
|
|
/* Enable support for rounding mode context. */
|
|
|
|
# define HAVE_RM_CTX 1
|
|
|
|
|
|
|
|
static __always_inline void
|
|
|
|
libc_feholdexcept_setround_mips_ctx (struct rm_ctx *ctx, int round)
|
|
|
|
{
|
|
|
|
fpu_control_t old, new;
|
|
|
|
|
|
|
|
/* Save the current state. */
|
|
|
|
_FPU_GETCW (old);
|
|
|
|
ctx->env.__fp_control_register = old;
|
|
|
|
|
|
|
|
/* Clear all exception enable bits and flags. */
|
|
|
|
new = old & ~(_FPU_MASK_ALL);
|
|
|
|
|
|
|
|
/* Set rounding bits. */
|
|
|
|
new = (new & ~_FPU_RC_MASK) | round;
|
|
|
|
|
|
|
|
if (__glibc_unlikely (new != old))
|
|
|
|
{
|
|
|
|
_FPU_SETCW (new);
|
|
|
|
ctx->updated_status = true;
|
|
|
|
}
|
|
|
|
else
|
|
|
|
ctx->updated_status = false;
|
|
|
|
}
|
|
|
|
# define libc_feholdexcept_setround_ctx libc_feholdexcept_setround_mips_ctx
|
|
|
|
# define libc_feholdexcept_setroundf_ctx libc_feholdexcept_setround_mips_ctx
|
|
|
|
# define libc_feholdexcept_setroundl_ctx libc_feholdexcept_setround_mips_ctx
|
|
|
|
|
|
|
|
static __always_inline void
|
|
|
|
libc_fesetenv_mips_ctx (struct rm_ctx *ctx)
|
|
|
|
{
|
|
|
|
libc_fesetenv_mips (&ctx->env);
|
|
|
|
}
|
|
|
|
# define libc_fesetenv_ctx libc_fesetenv_mips_ctx
|
|
|
|
# define libc_fesetenvf_ctx libc_fesetenv_mips_ctx
|
|
|
|
# define libc_fesetenvl_ctx libc_fesetenv_mips_ctx
|
|
|
|
|
|
|
|
static __always_inline void
|
|
|
|
libc_feupdateenv_mips_ctx (struct rm_ctx *ctx)
|
|
|
|
{
|
|
|
|
if (__glibc_unlikely (ctx->updated_status))
|
|
|
|
libc_feupdateenv_test_mips (&ctx->env, 0);
|
|
|
|
}
|
|
|
|
# define libc_feupdateenv_ctx libc_feupdateenv_mips_ctx
|
|
|
|
# define libc_feupdateenvf_ctx libc_feupdateenv_mips_ctx
|
|
|
|
# define libc_feupdateenvl_ctx libc_feupdateenv_mips_ctx
|
2014-03-11 22:30:40 +00:00
|
|
|
# define libc_feresetround_ctx libc_feupdateenv_mips_ctx
|
|
|
|
# define libc_feresetroundf_ctx libc_feupdateenv_mips_ctx
|
|
|
|
# define libc_feresetroundl_ctx libc_feupdateenv_mips_ctx
|
2014-02-26 22:19:20 +00:00
|
|
|
|
|
|
|
static __always_inline void
|
|
|
|
libc_feholdsetround_mips_ctx (struct rm_ctx *ctx, int round)
|
|
|
|
{
|
|
|
|
fpu_control_t old, new;
|
|
|
|
|
|
|
|
/* Save the current state. */
|
|
|
|
_FPU_GETCW (old);
|
|
|
|
ctx->env.__fp_control_register = old;
|
|
|
|
|
|
|
|
/* Set rounding bits. */
|
|
|
|
new = (old & ~_FPU_RC_MASK) | round;
|
|
|
|
|
|
|
|
if (__glibc_unlikely (new != old))
|
|
|
|
{
|
|
|
|
_FPU_SETCW (new);
|
|
|
|
ctx->updated_status = true;
|
|
|
|
}
|
|
|
|
else
|
|
|
|
ctx->updated_status = false;
|
|
|
|
}
|
|
|
|
# define libc_feholdsetround_ctx libc_feholdsetround_mips_ctx
|
|
|
|
# define libc_feholdsetroundf_ctx libc_feholdsetround_mips_ctx
|
|
|
|
# define libc_feholdsetroundl_ctx libc_feholdsetround_mips_ctx
|
|
|
|
|
2013-09-23 16:58:30 +00:00
|
|
|
#endif
|
|
|
|
|
2013-04-02 11:51:02 +00:00
|
|
|
#include_next <math_private.h>
|
|
|
|
|
|
|
|
#endif
|