343be16fcbfe4e9a018eeb2d1a156d16f10b7fd3
[kopensolaris-gnu/glibc.git] / sysdeps / powerpc / fenv_libc.h
1 /* Internal libc stuff for floating point environment routines.
2    Copyright (C) 1997 Free Software Foundation, Inc.
3    This file is part of the GNU C Library.
4
5    The GNU C Library is free software; you can redistribute it and/or
6    modify it under the terms of the GNU Library General Public License as
7    published by the Free Software Foundation; either version 2 of the
8    License, or (at your option) any later version.
9
10    The GNU C Library is distributed in the hope that it will be useful,
11    but WITHOUT ANY WARRANTY; without even the implied warranty of
12    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
13    Library General Public License for more details.
14
15    You should have received a copy of the GNU Library General Public
16    License along with the GNU C Library; see the file COPYING.LIB.  If not,
17    write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
18    Boston, MA 02111-1307, USA.  */
19
20 #ifndef _FENV_LIBC_H
21 #define _FENV_LIBC_H    1
22
23 #include <fenv.h>
24
25 /* The sticky bits in the FPSCR indicating exceptions have occurred.  */
26 #define FPSCR_STICKY_BITS ((FE_ALL_EXCEPT | FE_ALL_INVALID) & ~FE_INVALID)
27
28 /* Equivalent to fegetenv, but returns a fenv_t instead of taking a
29    pointer.  */
30 #define fegetenv_register() \
31         ({ fenv_t env; asm volatile ("mffs %0" : "=f" (env)); env; })
32
33 /* Equivalent to fesetenv, but takes a fenv_t instead of a pointer.  */
34 #define fesetenv_register(env) \
35         ({ double d = (env); asm volatile ("mtfsf 0xff,%0" : : "f" (d)); })
36
37 /* This very handy macro:
38    - Sets the rounding mode to 'round to nearest';
39    - Sets the processor into IEEE mode; and
40    - Prevents exceptions from being raised for inexact results.
41    These things happen to be exactly what you need for typical elementary
42    functions.  */
43 #define relax_fenv_state() asm ("mtfsfi 7,0")
44
45 /* Set/clear a particular FPSCR bit (for instance,
46    reset_fpscr_bit(FPSCR_VE);
47    prevents INVALID exceptions from being raised).  */
48 #define set_fpscr_bit(x) asm volatile ("mtfsb1 %0" : : "i"(x))
49 #define reset_fpscr_bit(x) asm volatile ("mtfsb0 %0" : : "i"(x))
50
51 typedef union
52 {
53   fenv_t fenv;
54   unsigned int l[2];
55 } fenv_union_t;
56
57 /* Definitions of all the FPSCR bit numbers */
58 enum {
59   FPSCR_FX = 0,    /* exception summary */
60   FPSCR_FEX,       /* enabled exception summary */
61   FPSCR_VX,        /* invalid operation summary */
62   FPSCR_OX,        /* overflow */
63   FPSCR_UX,        /* underflow */
64   FPSCR_ZX,        /* zero divide */
65   FPSCR_XX,        /* inexact */
66   FPSCR_VXSNAN,    /* invalid operation for SNaN */
67   FPSCR_VXISI,     /* invalid operation for Inf-Inf */
68   FPSCR_VXIDI,     /* invalid operation for Inf/Inf */
69   FPSCR_VXZDZ,     /* invalid operation for 0/0 */
70   FPSCR_VXIMZ,     /* invalid operation for Inf*0 */
71   FPSCR_VXVC,      /* invalid operation for invalid compare */
72   FPSCR_FR,        /* fraction rounded [fraction was incremented by round] */
73   FPSCR_FI,        /* fraction inexact */
74   FPSCR_FPRF_C,    /* result class descriptor */
75   FPSCR_FPRF_FL,   /* result less than (usually, less than 0) */
76   FPSCR_FPRF_FG,   /* result greater than */
77   FPSCR_FPRF_FE,   /* result equal to */
78   FPSCR_FPRF_FU,   /* result unordered */
79   FPSCR_20,        /* reserved */
80   FPSCR_VXSOFT,    /* invalid operation set by software */
81   FPSCR_VXSQRT,    /* invalid operation for square root */
82   FPSCR_VXCVI,     /* invalid operation for invalid integer convert */
83   FPSCR_VE,        /* invalid operation exception enable */
84   FPSCR_OE,        /* overflow exception enable */
85   FPSCR_UE,        /* underflow exception enable */
86   FPSCR_ZE,        /* zero divide exception enable */
87   FPSCR_XE,        /* inexact exception enable */
88   FPSCR_NI         /* non-IEEE mode (typically, no denormalised numbers) */
89   /* the remaining two least-significant bits keep the rounding mode */
90 };
91
92 /* This operation (i) sets the appropriate FPSCR bits for its
93    parameter, (ii) converts SNaN to the corresponding NaN, and (iii)
94    otherwise passes its parameter through unchanged (in particular, -0
95    and +0 stay as they were).  The `obvious' way to do this is optimised
96    out by gcc.  */
97 #define f_wash(x) \
98    ({ double d; asm volatile ("fmul %0,%1,%2" \
99                               : "=f"(d) \
100                               : "f" (x), "f"((float)1.0)); d; })
101 #define f_washf(x) \
102    ({ float f; asm volatile ("fmuls %0,%1,%2" \
103                              : "=f"(f) \
104                              : "f" (x), "f"((float)1.0)); f; })
105  
106 #endif /* fenv_libc.h */