fixunstfdi_ieee754.c revision 1.2 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