1 /*- 2 * Copyright (c) 2003 Mike Barcroft <mike@FreeBSD.org> 3 * Copyright (c) 2002 David Schultz <das@FreeBSD.ORG> 4 * All rights reserved. 5 * 6 * Redistribution and use in source and binary forms, with or without 7 * modification, are permitted provided that the following conditions 8 * are met: 9 * 1. Redistributions of source code must retain the above copyright 10 * notice, this list of conditions and the following disclaimer. 11 * 2. Redistributions in binary form must reproduce the above copyright 12 * notice, this list of conditions and the following disclaimer in the 13 * documentation and/or other materials provided with the distribution. 14 * 15 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND 16 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 17 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 18 * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE 19 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 20 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS 21 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 22 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 23 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY 24 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 25 * SUCH DAMAGE. 26 * 27 * $FreeBSD: src/lib/libc/include/fpmath.h,v 1.3 2005/02/06 03:23:31 das Exp $ 28 */ 29 30 #include <endian.h> 31 #if defined(__arm__) 32 #include "arm/_fpmath.h" 33 #elif defined(__aarch64__) 34 #include "arm64/_fpmath.h" 35 #elif defined(__i386__) 36 #include "i386/_fpmath.h" 37 #elif defined(__x86_64__) 38 #include "amd64/_fpmath.h" 39 #elif defined(__mips__) 40 #include "mips/_fpmath.h" 41 #elif !defined(__le32__) && !defined(__le64__) 42 #error Unknown arch! 43 #endif 44 45 union IEEEf2bits { 46 float f; 47 struct { 48 #if __BYTE_ORDER == __LITTLE_ENDIAN 49 unsigned int man :23; 50 unsigned int exp :8; 51 unsigned int sign :1; 52 #else /* _BIG_ENDIAN */ 53 unsigned int sign :1; 54 unsigned int exp :8; 55 unsigned int man :23; 56 #endif 57 } bits; 58 }; 59 60 #define DBL_MANH_SIZE 20 61 #define DBL_MANL_SIZE 32 62 63 union IEEEd2bits { 64 double d; 65 struct { 66 /* #ifdef __ARMEB__ */ 67 #if (__BYTE_ORDER == __BIG_ENDIAN) || (defined(__arm__) && !defined(__VFP_FP__)) 68 unsigned int manh :20; 69 unsigned int exp :11; 70 unsigned int sign :1; 71 unsigned int manl :32; 72 #elif __BYTE_ORDER == __LITTLE_ENDIAN 73 unsigned int manl :32; 74 unsigned int manh :20; 75 unsigned int exp :11; 76 unsigned int sign :1; 77 #elif __BYTE_ORDER == __BIG_ENDIAN 78 unsigned int sign :1; 79 unsigned int exp :11; 80 unsigned int manh :20; 81 unsigned int manl :32; 82 #endif 83 } bits; 84 }; 85 86 /* 87 * The BSD "long double" functions are broken when sizeof(long double) == sizeof(double). 88 * Android works around those cases by replacing the broken functions with our own trivial stubs 89 * that call the regular "double" function. 90 */ 91 #define __fpclassifyl __broken__fpclassify 92 #define __isfinitel __broken__isfinitel 93 #define __isinfl __broken__isinfl 94 #define __isnanl __broken__isnanl 95 #define __isnormall __broken__isnormall 96 #define __signbitl __broken_signbitl 97