1 /* Definitions of libc internal inline math functions implemented 2 by the m68881/2. 3 Copyright (C) 1991,92,93,94,96,97,98,99 Free Software Foundation, Inc. 4 This file is part of the GNU C Library. 5 6 The GNU C Library is free software; you can redistribute it and/or 7 modify it under the terms of the GNU Lesser General Public 8 License as published by the Free Software Foundation; either 9 version 2.1 of the License, or (at your option) any later version. 10 11 The GNU C Library is distributed in the hope that it will be useful, 12 but WITHOUT ANY WARRANTY; without even the implied warranty of 13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 14 Lesser General Public License for more details. 15 16 You should have received a copy of the GNU Lesser General Public 17 License along with the GNU C Library; if not, write to the Free 18 Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 19 02111-1307 USA. */ 20 21 /* This file contains the definitions of the inline math functions that 22 are only used internally inside libm, not visible to the user. */ 23 24 __inline_mathop (__ieee754_acos, acos) 25 __inline_mathop (__ieee754_asin, asin) 26 __inline_mathop (__ieee754_cosh, cosh) 27 __inline_mathop (__ieee754_sinh, sinh) 28 __inline_mathop (__ieee754_exp, etox) 29 __inline_mathop (__ieee754_exp2, twotox) 30 __inline_mathop (__ieee754_exp10, tentox) 31 __inline_mathop (__ieee754_log10, log10) 32 __inline_mathop (__ieee754_log2, log2) 33 __inline_mathop (__ieee754_log, logn) 34 __inline_mathop (__ieee754_sqrt, sqrt) 35 __inline_mathop (__ieee754_atanh, atanh) 36 37 __m81_defun (double, __ieee754_remainder, (double __x, double __y)) 38 { 39 double __result; 40 __asm ("frem%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); 41 return __result; 42 } 43 44 __m81_defun (float, __ieee754_remainderf, (float __x, float __y)) 45 { 46 float __result; 47 __asm ("frem%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); 48 return __result; 49 } 50 51 __m81_defun (long double, 52 __ieee754_remainderl, (long double __x, long double __y)) 53 { 54 long double __result; 55 __asm ("frem%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); 56 return __result; 57 } 58 59 __m81_defun (double, __ieee754_fmod, (double __x, double __y)) 60 { 61 double __result; 62 __asm ("fmod%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); 63 return __result; 64 } 65 66 __m81_defun (float, __ieee754_fmodf, (float __x, float __y)) 67 { 68 float __result; 69 __asm ("fmod%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); 70 return __result; 71 } 72 73 __m81_defun (long double, 74 __ieee754_fmodl, (long double __x, long double __y)) 75 { 76 long double __result; 77 __asm ("fmod%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); 78 return __result; 79 } 80 81 /* Get the m68881 condition codes, to quickly check multiple conditions. */ 82 static __inline__ unsigned long 83 __m81_test (long double __val) 84 { 85 unsigned long __fpsr; 86 __asm ("ftst%.x %1; fmove%.l %/fpsr,%0" : "=dm" (__fpsr) : "f" (__val)); 87 return __fpsr; 88 } 89 90 /* Bit values returned by __m81_test. */ 91 #define __M81_COND_NAN (1 << 24) 92 #define __M81_COND_INF (2 << 24) 93 #define __M81_COND_ZERO (4 << 24) 94 #define __M81_COND_NEG (8 << 24) 95