| 1 | /* -*- buffer-read-only: t -*- vi: set ro: */
|
|---|
| 2 | /* DO NOT EDIT! GENERATED AUTOMATICALLY! */
|
|---|
| 3 | /* Manipulating the FPU control word.
|
|---|
| 4 | Copyright (C) 2007-2011 Free Software Foundation, Inc.
|
|---|
| 5 | Written by Bruno Haible <bruno@clisp.org>, 2007.
|
|---|
| 6 |
|
|---|
| 7 | This program is free software: you can redistribute it and/or modify
|
|---|
| 8 | it under the terms of the GNU General Public License as published by
|
|---|
| 9 | the Free Software Foundation; either version 3 of the License, or
|
|---|
| 10 | (at your option) any later version.
|
|---|
| 11 |
|
|---|
| 12 | This program is distributed in the hope that it will be useful,
|
|---|
| 13 | but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|---|
| 14 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|---|
| 15 | GNU General Public License for more details.
|
|---|
| 16 |
|
|---|
| 17 | You should have received a copy of the GNU General Public License
|
|---|
| 18 | along with this program. If not, see <http://www.gnu.org/licenses/>. */
|
|---|
| 19 |
|
|---|
| 20 | #ifndef _FPUCW_H
|
|---|
| 21 | #define _FPUCW_H
|
|---|
| 22 |
|
|---|
| 23 | /* The i386 floating point hardware (the 387 compatible FPU, not the modern
|
|---|
| 24 | SSE/SSE2 hardware) has a controllable rounding precision. It is specified
|
|---|
| 25 | through the 'PC' bits in the FPU control word ('fctrl' register). (See
|
|---|
| 26 | the GNU libc i386 <fpu_control.h> header for details.)
|
|---|
| 27 |
|
|---|
| 28 | On some platforms, such as Linux or Solaris, the default precision setting
|
|---|
| 29 | is set to "extended precision". This means that 'long double' instructions
|
|---|
| 30 | operate correctly, but 'double' computations often produce slightly
|
|---|
| 31 | different results as on strictly IEEE 754 conforming systems.
|
|---|
| 32 |
|
|---|
| 33 | On some platforms, such as NetBSD, the default precision is set to
|
|---|
| 34 | "double precision". This means that 'long double' instructions will operate
|
|---|
| 35 | only as 'double', i.e. lead wrong results.
|
|---|
| 36 |
|
|---|
| 37 | The FPU control word is under control of the application, i.e. it is
|
|---|
| 38 | not required to be set either way by the ABI. (In fact, the i386 ABI
|
|---|
| 39 | http://refspecs.freestandards.org/elf/abi386-4.pdf page 3-12 = page 38
|
|---|
| 40 | is not clear about it. But in any case, gcc treats the control word
|
|---|
| 41 | like a "preserved" register: it emits code that assumes that the control
|
|---|
| 42 | word is preserved across calls, and it restores the control word at the
|
|---|
| 43 | end of functions that modify it.)
|
|---|
| 44 |
|
|---|
| 45 | See Vincent LefÚvre's page http://www.vinc17.org/research/extended.en.html
|
|---|
| 46 | for a good explanation.
|
|---|
| 47 | See http://www.uwsg.iu.edu/hypermail/linux/kernel/0103.0/0453.html for
|
|---|
| 48 | some argumentation which setting should be the default. */
|
|---|
| 49 |
|
|---|
| 50 | /* This header file provides the following facilities:
|
|---|
| 51 | fpucw_t integral type holding the value of 'fctrl'
|
|---|
| 52 | FPU_PC_MASK bit mask denoting the precision control
|
|---|
| 53 | FPU_PC_DOUBLE precision control for 53 bits mantissa
|
|---|
| 54 | FPU_PC_EXTENDED precision control for 64 bits mantissa
|
|---|
| 55 | GET_FPUCW () yields the current FPU control word
|
|---|
| 56 | SET_FPUCW (word) sets the FPU control word
|
|---|
| 57 | DECL_LONG_DOUBLE_ROUNDING variable declaration for
|
|---|
| 58 | BEGIN/END_LONG_DOUBLE_ROUNDING
|
|---|
| 59 | BEGIN_LONG_DOUBLE_ROUNDING () starts a sequence of instructions with
|
|---|
| 60 | 'long double' safe operation precision
|
|---|
| 61 | END_LONG_DOUBLE_ROUNDING () ends a sequence of instructions with
|
|---|
| 62 | 'long double' safe operation precision
|
|---|
| 63 | */
|
|---|
| 64 |
|
|---|
| 65 | /* Inline assembler like this works only with GNU C. */
|
|---|
| 66 | #if (defined __i386__ || defined __x86_64__) && defined __GNUC__
|
|---|
| 67 |
|
|---|
| 68 | typedef unsigned short fpucw_t; /* glibc calls this fpu_control_t */
|
|---|
| 69 |
|
|---|
| 70 | # define FPU_PC_MASK 0x0300
|
|---|
| 71 | # define FPU_PC_DOUBLE 0x200 /* glibc calls this _FPU_DOUBLE */
|
|---|
| 72 | # define FPU_PC_EXTENDED 0x300 /* glibc calls this _FPU_EXTENDED */
|
|---|
| 73 |
|
|---|
| 74 | # define GET_FPUCW() \
|
|---|
| 75 | ({ fpucw_t _cw; \
|
|---|
| 76 | __asm__ __volatile__ ("fnstcw %0" : "=m" (*&_cw)); \
|
|---|
| 77 | _cw; \
|
|---|
| 78 | })
|
|---|
| 79 | # define SET_FPUCW(word) \
|
|---|
| 80 | (void)({ fpucw_t _ncw = (word); \
|
|---|
| 81 | __asm__ __volatile__ ("fldcw %0" : : "m" (*&_ncw)); \
|
|---|
| 82 | })
|
|---|
| 83 |
|
|---|
| 84 | # define DECL_LONG_DOUBLE_ROUNDING \
|
|---|
| 85 | fpucw_t oldcw;
|
|---|
| 86 | # define BEGIN_LONG_DOUBLE_ROUNDING() \
|
|---|
| 87 | (void)(oldcw = GET_FPUCW (), \
|
|---|
| 88 | SET_FPUCW ((oldcw & ~FPU_PC_MASK) | FPU_PC_EXTENDED))
|
|---|
| 89 | # define END_LONG_DOUBLE_ROUNDING() \
|
|---|
| 90 | SET_FPUCW (oldcw)
|
|---|
| 91 |
|
|---|
| 92 | #else
|
|---|
| 93 |
|
|---|
| 94 | typedef unsigned int fpucw_t;
|
|---|
| 95 |
|
|---|
| 96 | # define FPU_PC_MASK 0
|
|---|
| 97 | # define FPU_PC_DOUBLE 0
|
|---|
| 98 | # define FPU_PC_EXTENDED 0
|
|---|
| 99 |
|
|---|
| 100 | # define GET_FPUCW() 0
|
|---|
| 101 | # define SET_FPUCW(word) (void)(word)
|
|---|
| 102 |
|
|---|
| 103 | # define DECL_LONG_DOUBLE_ROUNDING
|
|---|
| 104 | # define BEGIN_LONG_DOUBLE_ROUNDING()
|
|---|
| 105 | # define END_LONG_DOUBLE_ROUNDING()
|
|---|
| 106 |
|
|---|
| 107 | #endif
|
|---|
| 108 |
|
|---|
| 109 | #endif /* _FPUCW_H */
|
|---|
