1 1.2 rin /* $NetBSD: fixunstfdi_ieee754.c,v 1.2 2024/01/05 09:41:58 rin Exp $ */ 2 1.1 matt /*- 3 1.1 matt * Copyright (c) 2011 The NetBSD Foundation, Inc. 4 1.1 matt * All rights reserved. 5 1.1 matt * 6 1.1 matt * This code is derived from software contributed to The NetBSD Foundation 7 1.1 matt * by Matt Thomas of 3am Software Foundry. 8 1.1 matt * 9 1.1 matt * Redistribution and use in source and binary forms, with or without 10 1.1 matt * modification, are permitted provided that the following conditions 11 1.1 matt * are met: 12 1.1 matt * 1. Redistributions of source code must retain the above copyright 13 1.1 matt * notice, this list of conditions and the following disclaimer. 14 1.1 matt * 2. Redistributions in binary form must reproduce the above copyright 15 1.1 matt * notice, this list of conditions and the following disclaimer in the 16 1.1 matt * documentation and/or other materials provided with the distribution. 17 1.1 matt * 18 1.1 matt * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS 19 1.1 matt * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED 20 1.1 matt * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR 21 1.1 matt * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS 22 1.1 matt * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 23 1.1 matt * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF 24 1.1 matt * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS 25 1.1 matt * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN 26 1.1 matt * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) 27 1.1 matt * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 28 1.1 matt * POSSIBILITY OF SUCH DAMAGE. 29 1.1 matt */ 30 1.1 matt 31 1.1 matt #include <sys/cdefs.h> 32 1.1 matt 33 1.1 matt #if defined(LIBC_SCCS) && !defined(lint) 34 1.2 rin __RCSID("$NetBSD: fixunstfdi_ieee754.c,v 1.2 2024/01/05 09:41:58 rin Exp $"); 35 1.1 matt #endif /* LIBC_SCCS and not lint */ 36 1.1 matt 37 1.1 matt #include <stdbool.h> 38 1.2 rin #include <stddef.h> 39 1.1 matt #include <stdint.h> 40 1.1 matt #include <float.h> 41 1.1 matt #include <machine/ieee.h> 42 1.1 matt 43 1.1 matt #if defined(__x86_64__) || defined(__i486__) 44 1.1 matt #define FIXUNS __fixunsxfdi 45 1.1 matt #else 46 1.1 matt #define FIXUNS __fixunstfdi 47 1.1 matt #endif 48 1.1 matt 49 1.1 matt uint64_t __fixunsgen64(int, bool, size_t, size_t, const uint32_t *); 50 1.1 matt 51 1.1 matt uint64_t FIXUNS(long double); 52 1.1 matt 53 1.1 matt /* 54 1.1 matt * Convert long double to uint64_t. All operations are done module 2^64. 55 1.1 matt */ 56 1.1 matt uint64_t 57 1.1 matt FIXUNS(long double x) 58 1.1 matt { 59 1.1 matt const union ieee_ext_u extu = { .extu_ld = x }; 60 1.1 matt uint32_t frac[(EXT_FRACBITS + 31)/32 + 2]; 61 1.1 matt 62 1.1 matt frac[0] = 0; 63 1.1 matt frac[1] = 0; 64 1.1 matt 65 1.1 matt EXT_TO_ARRAY32(extu, &frac[2]); 66 1.1 matt 67 1.1 matt return __fixunsgen64( 68 1.1 matt extu.extu_ext.ext_exp - EXT_EXP_BIAS, 69 1.1 matt extu.extu_ext.ext_sign != 0, 70 1.1 matt LDBL_MANT_DIG, 71 1.1 matt EXT_FRACHBITS, 72 1.1 matt &frac[__arraycount(frac)-1]); 73 1.1 matt } 74