fenv.h revision 1.1.2.2
1 1.1.2.2 pgoyette /*- 2 1.1.2.2 pgoyette * Copyright (c) 2004-2005 David Schultz <das (at) FreeBSD.ORG> 3 1.1.2.2 pgoyette * All rights reserved. 4 1.1.2.2 pgoyette * 5 1.1.2.2 pgoyette * Redistribution and use in source and binary forms, with or without 6 1.1.2.2 pgoyette * modification, are permitted provided that the following conditions 7 1.1.2.2 pgoyette * are met: 8 1.1.2.2 pgoyette * 1. Redistributions of source code must retain the above copyright 9 1.1.2.2 pgoyette * notice, this list of conditions and the following disclaimer. 10 1.1.2.2 pgoyette * 2. Redistributions in binary form must reproduce the above copyright 11 1.1.2.2 pgoyette * notice, this list of conditions and the following disclaimer in the 12 1.1.2.2 pgoyette * documentation and/or other materials provided with the distribution. 13 1.1.2.2 pgoyette * 14 1.1.2.2 pgoyette * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND 15 1.1.2.2 pgoyette * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 16 1.1.2.2 pgoyette * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 17 1.1.2.2 pgoyette * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE 18 1.1.2.2 pgoyette * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 19 1.1.2.2 pgoyette * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS 20 1.1.2.2 pgoyette * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 21 1.1.2.2 pgoyette * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 22 1.1.2.2 pgoyette * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY 23 1.1.2.2 pgoyette * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 24 1.1.2.2 pgoyette * SUCH DAMAGE. 25 1.1.2.2 pgoyette * 26 1.1.2.2 pgoyette * $FreeBSD: releng/10.1/lib/msun/ia64/fenv.h 226218 2011-10-10 15:43:09Z das $ 27 1.1.2.2 pgoyette */ 28 1.1.2.2 pgoyette 29 1.1.2.2 pgoyette #ifndef _IA64_FENV_H_ 30 1.1.2.2 pgoyette #define _IA64_FENV_H_ 31 1.1.2.2 pgoyette 32 1.1.2.2 pgoyette #include <sys/stdint.h> 33 1.1.2.2 pgoyette 34 1.1.2.2 pgoyette #ifndef __fenv_static 35 1.1.2.2 pgoyette #define __fenv_static static 36 1.1.2.2 pgoyette #endif 37 1.1.2.2 pgoyette 38 1.1.2.2 pgoyette typedef __uint64_t fenv_t; 39 1.1.2.2 pgoyette typedef __uint16_t fexcept_t; 40 1.1.2.2 pgoyette 41 1.1.2.2 pgoyette /* Exception flags */ 42 1.1.2.2 pgoyette #define FE_INVALID 0x01 43 1.1.2.2 pgoyette #define FE_DENORMAL 0x02 44 1.1.2.2 pgoyette #define FE_DIVBYZERO 0x04 45 1.1.2.2 pgoyette #define FE_OVERFLOW 0x08 46 1.1.2.2 pgoyette #define FE_UNDERFLOW 0x10 47 1.1.2.2 pgoyette #define FE_INEXACT 0x20 48 1.1.2.2 pgoyette #define FE_ALL_EXCEPT (FE_DIVBYZERO | FE_DENORMAL | FE_INEXACT | \ 49 1.1.2.2 pgoyette FE_INVALID | FE_OVERFLOW | FE_UNDERFLOW) 50 1.1.2.2 pgoyette 51 1.1.2.2 pgoyette /* Rounding modes */ 52 1.1.2.2 pgoyette #define FE_TONEAREST 0x0000 53 1.1.2.2 pgoyette #define FE_DOWNWARD 0x0400 54 1.1.2.2 pgoyette #define FE_UPWARD 0x0800 55 1.1.2.2 pgoyette #define FE_TOWARDZERO 0x0c00 56 1.1.2.2 pgoyette #define _ROUND_MASK (FE_TONEAREST | FE_DOWNWARD | \ 57 1.1.2.2 pgoyette FE_UPWARD | FE_TOWARDZERO) 58 1.1.2.2 pgoyette 59 1.1.2.2 pgoyette __BEGIN_DECLS 60 1.1.2.2 pgoyette 61 1.1.2.2 pgoyette /* Default floating-point environment */ 62 1.1.2.2 pgoyette extern const fenv_t __fe_dfl_env; 63 1.1.2.2 pgoyette #define FE_DFL_ENV (&__fe_dfl_env) 64 1.1.2.2 pgoyette 65 1.1.2.2 pgoyette #define _FPUSW_SHIFT 13 66 1.1.2.2 pgoyette 67 1.1.2.2 pgoyette #define __stfpsr(__r) __asm __volatile("mov %0=ar.fpsr" : "=r" (*(__r))) 68 1.1.2.2 pgoyette #define __ldfpsr(__r) __asm __volatile("mov ar.fpsr=%0;;" : : "r" (__r)) 69 1.1.2.2 pgoyette 70 1.1.2.2 pgoyette __fenv_static inline int 71 1.1.2.2 pgoyette feclearexcept(int __excepts) 72 1.1.2.2 pgoyette { 73 1.1.2.2 pgoyette fenv_t __fpsr; 74 1.1.2.2 pgoyette 75 1.1.2.2 pgoyette __stfpsr(&__fpsr); 76 1.1.2.2 pgoyette __fpsr &= ~((fenv_t)__excepts << _FPUSW_SHIFT); 77 1.1.2.2 pgoyette __ldfpsr(__fpsr); 78 1.1.2.2 pgoyette return (0); 79 1.1.2.2 pgoyette } 80 1.1.2.2 pgoyette 81 1.1.2.2 pgoyette __fenv_static inline int 82 1.1.2.2 pgoyette fegetexceptflag(fexcept_t *__flagp, int __excepts) 83 1.1.2.2 pgoyette { 84 1.1.2.2 pgoyette fenv_t __fpsr; 85 1.1.2.2 pgoyette 86 1.1.2.2 pgoyette __stfpsr(&__fpsr); 87 1.1.2.2 pgoyette *__flagp = (fexcept_t)(__fpsr >> _FPUSW_SHIFT) & __excepts; 88 1.1.2.2 pgoyette return (0); 89 1.1.2.2 pgoyette } 90 1.1.2.2 pgoyette 91 1.1.2.2 pgoyette __fenv_static inline int 92 1.1.2.2 pgoyette fesetexceptflag(const fexcept_t *__flagp, int __excepts) 93 1.1.2.2 pgoyette { 94 1.1.2.2 pgoyette fenv_t __fpsr; 95 1.1.2.2 pgoyette 96 1.1.2.2 pgoyette __stfpsr(&__fpsr); 97 1.1.2.2 pgoyette __fpsr &= ~((fenv_t)__excepts << _FPUSW_SHIFT); 98 1.1.2.2 pgoyette __fpsr |= (fenv_t)(__excepts & *__flagp) << _FPUSW_SHIFT; 99 1.1.2.2 pgoyette __ldfpsr(__fpsr); 100 1.1.2.2 pgoyette return (0); 101 1.1.2.2 pgoyette } 102 1.1.2.2 pgoyette 103 1.1.2.2 pgoyette /* 104 1.1.2.2 pgoyette * It is worthwhile to use the inline version of this function iff it 105 1.1.2.2 pgoyette * is called with arguments that are compile-time constants (due to 106 1.1.2.2 pgoyette * dead code elimination). Unfortunately, gcc isn't smart enough to 107 1.1.2.2 pgoyette * figure this out automatically, and there's no way to tell it. 108 1.1.2.2 pgoyette * We assume that constant arguments will be the common case. 109 1.1.2.2 pgoyette */ 110 1.1.2.2 pgoyette __fenv_static inline int 111 1.1.2.2 pgoyette feraiseexcept(int __excepts) 112 1.1.2.2 pgoyette { 113 1.1.2.2 pgoyette volatile double d; 114 1.1.2.2 pgoyette 115 1.1.2.2 pgoyette /* 116 1.1.2.2 pgoyette * With a compiler that supports the FENV_ACCESS pragma 117 1.1.2.2 pgoyette * properly, simple expressions like '0.0 / 0.0' should 118 1.1.2.2 pgoyette * be sufficient to generate traps. Unfortunately, we 119 1.1.2.2 pgoyette * need to bring a volatile variable into the equation 120 1.1.2.2 pgoyette * to prevent incorrect optimizations. 121 1.1.2.2 pgoyette */ 122 1.1.2.2 pgoyette if (__excepts & FE_INVALID) { 123 1.1.2.2 pgoyette d = 0.0; 124 1.1.2.2 pgoyette d = 0.0 / d; 125 1.1.2.2 pgoyette } 126 1.1.2.2 pgoyette if (__excepts & FE_DIVBYZERO) { 127 1.1.2.2 pgoyette d = 0.0; 128 1.1.2.2 pgoyette d = 1.0 / d; 129 1.1.2.2 pgoyette } 130 1.1.2.2 pgoyette if (__excepts & FE_OVERFLOW) { 131 1.1.2.2 pgoyette d = 0x1.ffp1023; 132 1.1.2.2 pgoyette d *= 2.0; 133 1.1.2.2 pgoyette } 134 1.1.2.2 pgoyette if (__excepts & FE_UNDERFLOW) { 135 1.1.2.2 pgoyette d = 0x1p-1022; 136 1.1.2.2 pgoyette d /= 0x1p1023; 137 1.1.2.2 pgoyette } 138 1.1.2.2 pgoyette if (__excepts & FE_INEXACT) { 139 1.1.2.2 pgoyette d = 0x1p-1022; 140 1.1.2.2 pgoyette d += 1.0; 141 1.1.2.2 pgoyette } 142 1.1.2.2 pgoyette return (0); 143 1.1.2.2 pgoyette } 144 1.1.2.2 pgoyette 145 1.1.2.2 pgoyette __fenv_static inline int 146 1.1.2.2 pgoyette fetestexcept(int __excepts) 147 1.1.2.2 pgoyette { 148 1.1.2.2 pgoyette fenv_t __fpsr; 149 1.1.2.2 pgoyette 150 1.1.2.2 pgoyette __stfpsr(&__fpsr); 151 1.1.2.2 pgoyette return ((__fpsr >> _FPUSW_SHIFT) & __excepts); 152 1.1.2.2 pgoyette } 153 1.1.2.2 pgoyette 154 1.1.2.2 pgoyette 155 1.1.2.2 pgoyette __fenv_static inline int 156 1.1.2.2 pgoyette fegetround(void) 157 1.1.2.2 pgoyette { 158 1.1.2.2 pgoyette fenv_t __fpsr; 159 1.1.2.2 pgoyette 160 1.1.2.2 pgoyette __stfpsr(&__fpsr); 161 1.1.2.2 pgoyette return (__fpsr & _ROUND_MASK); 162 1.1.2.2 pgoyette } 163 1.1.2.2 pgoyette 164 1.1.2.2 pgoyette __fenv_static inline int 165 1.1.2.2 pgoyette fesetround(int __round) 166 1.1.2.2 pgoyette { 167 1.1.2.2 pgoyette fenv_t __fpsr; 168 1.1.2.2 pgoyette 169 1.1.2.2 pgoyette if (__round & ~_ROUND_MASK) 170 1.1.2.2 pgoyette return (-1); 171 1.1.2.2 pgoyette __stfpsr(&__fpsr); 172 1.1.2.2 pgoyette __fpsr &= ~_ROUND_MASK; 173 1.1.2.2 pgoyette __fpsr |= __round; 174 1.1.2.2 pgoyette __ldfpsr(__fpsr); 175 1.1.2.2 pgoyette return (0); 176 1.1.2.2 pgoyette } 177 1.1.2.2 pgoyette 178 1.1.2.2 pgoyette __fenv_static inline int 179 1.1.2.2 pgoyette fegetenv(fenv_t *__envp) 180 1.1.2.2 pgoyette { 181 1.1.2.2 pgoyette 182 1.1.2.2 pgoyette __stfpsr(__envp); 183 1.1.2.2 pgoyette return (0); 184 1.1.2.2 pgoyette } 185 1.1.2.2 pgoyette 186 1.1.2.2 pgoyette __fenv_static inline int 187 1.1.2.2 pgoyette feholdexcept(fenv_t *__envp) 188 1.1.2.2 pgoyette { 189 1.1.2.2 pgoyette fenv_t __fpsr; 190 1.1.2.2 pgoyette 191 1.1.2.2 pgoyette __stfpsr(&__fpsr); 192 1.1.2.2 pgoyette *__envp = __fpsr; 193 1.1.2.2 pgoyette __fpsr &= ~((fenv_t)FE_ALL_EXCEPT << _FPUSW_SHIFT); 194 1.1.2.2 pgoyette __fpsr |= FE_ALL_EXCEPT; 195 1.1.2.2 pgoyette __ldfpsr(__fpsr); 196 1.1.2.2 pgoyette return (0); 197 1.1.2.2 pgoyette } 198 1.1.2.2 pgoyette 199 1.1.2.2 pgoyette __fenv_static inline int 200 1.1.2.2 pgoyette fesetenv(const fenv_t *__envp) 201 1.1.2.2 pgoyette { 202 1.1.2.2 pgoyette 203 1.1.2.2 pgoyette __ldfpsr(*__envp); 204 1.1.2.2 pgoyette return (0); 205 1.1.2.2 pgoyette } 206 1.1.2.2 pgoyette 207 1.1.2.2 pgoyette int feupdateenv(const fenv_t *__envp); 208 1.1.2.2 pgoyette 209 1.1.2.2 pgoyette #if defined(_NETBSD_SOURCE) || defined(_GNU_SOURCE) 210 1.1.2.2 pgoyette 211 1.1.2.2 pgoyette /* We currently provide no external definitions of the functions below. */ 212 1.1.2.2 pgoyette 213 1.1.2.2 pgoyette static inline int 214 1.1.2.2 pgoyette feenableexcept(int __mask) 215 1.1.2.2 pgoyette { 216 1.1.2.2 pgoyette fenv_t __newfpsr, __oldfpsr; 217 1.1.2.2 pgoyette 218 1.1.2.2 pgoyette __stfpsr(&__oldfpsr); 219 1.1.2.2 pgoyette __newfpsr = __oldfpsr & ~(__mask & FE_ALL_EXCEPT); 220 1.1.2.2 pgoyette __ldfpsr(__newfpsr); 221 1.1.2.2 pgoyette return (~__oldfpsr & FE_ALL_EXCEPT); 222 1.1.2.2 pgoyette } 223 1.1.2.2 pgoyette 224 1.1.2.2 pgoyette static inline int 225 1.1.2.2 pgoyette fedisableexcept(int __mask) 226 1.1.2.2 pgoyette { 227 1.1.2.2 pgoyette fenv_t __newfpsr, __oldfpsr; 228 1.1.2.2 pgoyette 229 1.1.2.2 pgoyette __stfpsr(&__oldfpsr); 230 1.1.2.2 pgoyette __newfpsr = __oldfpsr | (__mask & FE_ALL_EXCEPT); 231 1.1.2.2 pgoyette __ldfpsr(__newfpsr); 232 1.1.2.2 pgoyette return (~__oldfpsr & FE_ALL_EXCEPT); 233 1.1.2.2 pgoyette } 234 1.1.2.2 pgoyette 235 1.1.2.2 pgoyette static inline int 236 1.1.2.2 pgoyette fegetexcept(void) 237 1.1.2.2 pgoyette { 238 1.1.2.2 pgoyette fenv_t __fpsr; 239 1.1.2.2 pgoyette 240 1.1.2.2 pgoyette __stfpsr(&__fpsr); 241 1.1.2.2 pgoyette return (~__fpsr & FE_ALL_EXCEPT); 242 1.1.2.2 pgoyette } 243 1.1.2.2 pgoyette 244 1.1.2.2 pgoyette #endif /* _NETBSD_SOURCE || _GNU_SOURCE */ 245 1.1.2.2 pgoyette 246 1.1.2.2 pgoyette __END_DECLS 247 1.1.2.2 pgoyette 248 1.1.2.2 pgoyette #endif /* !_IA64_FENV_H_ */ 249
Indexes created Fri Oct 17 09:09:57 GMT 2025