1 /*- 2 * SPDX-License-Identifier: BSD-2-Clause 3 * 4 * Copyright (c) 2003 Mike Barcroft <mike@FreeBSD.org> 5 * Copyright (c) 2002 David Schultz <das@FreeBSD.ORG> 6 * All rights reserved. 7 * 8 * Redistribution and use in source and binary forms, with or without 9 * modification, are permitted provided that the following conditions 10 * are met: 11 * 1. Redistributions of source code must retain the above copyright 12 * notice, this list of conditions and the following disclaimer. 13 * 2. Redistributions in binary form must reproduce the above copyright 14 * notice, this list of conditions and the following disclaimer in the 15 * documentation and/or other materials provided with the distribution. 16 * 17 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND 18 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 19 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 20 * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE 21 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 22 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS 23 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 24 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 25 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY 26 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 27 * SUCH DAMAGE. 28 */ 29 30 #ifndef _FPMATH_H_ 31 #define _FPMATH_H_ 32 33 #include <endian.h> 34 #ifndef LD128_ENABLE 35 #include "_fpmath.h" 36 #endif 37 38 #ifndef _IEEE_WORD_ORDER 39 #define _IEEE_WORD_ORDER __BYTE_ORDER 40 #endif 41 42 #ifdef LD128_ENABLE 43 union IEEEl2bits { 44 long double e; 45 struct { 46 unsigned long manl :64; 47 unsigned long manh :48; 48 unsigned int exp :15; 49 unsigned int sign :1; 50 } bits; 51 /* TODO andrew: Check the packing here */ 52 struct { 53 unsigned long manl :64; 54 unsigned long manh :48; 55 unsigned int expsign :16; 56 } xbits; 57 }; 58 59 #define LDBL_NBIT 0 60 #define LDBL_IMPLICIT_NBIT 61 #define mask_nbit_l(u) ((void)0) 62 63 #define LDBL_MANH_SIZE 48 64 #define LDBL_MANL_SIZE 64 65 66 #define LDBL_TO_ARRAY32(u, a) do { \ 67 (a)[0] = (uint32_t)(u).bits.manl; \ 68 (a)[1] = (uint32_t)((u).bits.manl >> 32); \ 69 (a)[2] = (uint32_t)(u).bits.manh; \ 70 (a)[3] = (uint32_t)((u).bits.manh >> 32); \ 71 } while(0) 72 #else 73 union IEEEf2bits { 74 float f; 75 struct { 76 #if __BYTE_ORDER == __LITTLE_ENDIAN 77 unsigned int man :23; 78 unsigned int exp :8; 79 unsigned int sign :1; 80 #else /* __BIG_ENDIAN */ 81 unsigned int sign :1; 82 unsigned int exp :8; 83 unsigned int man :23; 84 #endif 85 } bits; 86 }; 87 88 #define DBL_MANH_SIZE 20 89 #define DBL_MANL_SIZE 32 90 91 union IEEEd2bits { 92 double d; 93 struct { 94 #if __BYTE_ORDER == __LITTLE_ENDIAN 95 #if _IEEE_WORD_ORDER == __LITTLE_ENDIAN 96 unsigned int manl :32; 97 #endif 98 unsigned int manh :20; 99 unsigned int exp :11; 100 unsigned int sign :1; 101 #if _IEEE_WORD_ORDER == __BIG_ENDIAN 102 unsigned int manl :32; 103 #endif 104 #else /* __BIG_ENDIAN */ 105 unsigned int sign :1; 106 unsigned int exp :11; 107 unsigned int manh :20; 108 unsigned int manl :32; 109 #endif 110 } bits; 111 }; 112 #endif 113 #endif /* !_FPMATH_H */