fixunsdfsi.S revision 1.6 1 1.6 joerg /* $NetBSD: fixunsdfsi.S,v 1.6 2013/09/12 15:36:16 joerg Exp $ */
2 1.4 christos
3 1.1 cgd /*
4 1.1 cgd * Copyright (c) 1992, 1993
5 1.1 cgd * The Regents of the University of California. All rights reserved.
6 1.1 cgd *
7 1.1 cgd * This software was developed by the Computer Systems Engineering group
8 1.1 cgd * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
9 1.1 cgd * contributed to Berkeley.
10 1.1 cgd *
11 1.1 cgd * Redistribution and use in source and binary forms, with or without
12 1.1 cgd * modification, are permitted provided that the following conditions
13 1.1 cgd * are met:
14 1.1 cgd * 1. Redistributions of source code must retain the above copyright
15 1.1 cgd * notice, this list of conditions and the following disclaimer.
16 1.1 cgd * 2. Redistributions in binary form must reproduce the above copyright
17 1.1 cgd * notice, this list of conditions and the following disclaimer in the
18 1.1 cgd * documentation and/or other materials provided with the distribution.
19 1.5 agc * 3. Neither the name of the University nor the names of its contributors
20 1.1 cgd * may be used to endorse or promote products derived from this software
21 1.1 cgd * without specific prior written permission.
22 1.1 cgd *
23 1.1 cgd * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
24 1.1 cgd * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
25 1.1 cgd * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
26 1.1 cgd * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
27 1.1 cgd * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
28 1.1 cgd * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
29 1.1 cgd * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
30 1.1 cgd * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
31 1.1 cgd * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
32 1.1 cgd * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
33 1.1 cgd * SUCH DAMAGE.
34 1.1 cgd *
35 1.1 cgd * from: Header: fixunsdfsi.s,v 1.3 91/10/08 00:03:15 torek Exp
36 1.1 cgd */
37 1.1 cgd
38 1.4 christos #include <machine/asm.h>
39 1.1 cgd #if defined(LIBC_SCCS) && !defined(lint)
40 1.4 christos #if 0
41 1.1 cgd .asciz "@(#)fixunsdfsi.s 8.1 (Berkeley) 6/4/93"
42 1.4 christos #else
43 1.6 joerg RCSID("$NetBSD: fixunsdfsi.S,v 1.6 2013/09/12 15:36:16 joerg Exp $")
44 1.4 christos #endif
45 1.1 cgd #endif /* LIBC_SCCS and not lint */
46 1.1 cgd
47 1.1 cgd /*
48 1.1 cgd * Convert double to unsigned integer (for gcc).
49 1.1 cgd *
50 1.1 cgd * I have made the output for NaN agree with the Sun compiler, not
51 1.1 cgd * that it really matters, by using `fbul,a'.
52 1.1 cgd */
53 1.1 cgd
54 1.1 cgd
55 1.1 cgd .align 8
56 1.1 cgd Lbig:
57 1.1 cgd .word 0x41e00000 ! .double 0r2147483648.0e+00
58 1.1 cgd .word 0 ! (who me, not trust the assembler?)
59 1.1 cgd
60 1.1 cgd ENTRY(__fixunsdfsi)
61 1.1 cgd sub %sp, 8, %sp
62 1.1 cgd std %o0, [%sp + 64] ! get argument into fpu reg
63 1.1 cgd ldd [%sp + 64], %f0
64 1.6 joerg #ifdef __PIC__
65 1.2 pk PICCY_SET(Lbig, %g1, %o0)
66 1.2 pk ldd [%g1], %f2
67 1.2 pk #else
68 1.1 cgd sethi %hi(Lbig), %g1
69 1.1 cgd ldd [%g1 + %lo(Lbig)], %f2
70 1.2 pk #endif
71 1.1 cgd fcmped %f0, %f2 ! d < 2^31, or NaN, or -Inf?
72 1.1 cgd nop ! (fpop2 delay)
73 1.1 cgd fbul,a 1f ! if so, use fdtoi to convert to int
74 1.1 cgd fdtoi %f0, %f0 ! (this includes negatives!)
75 1.1 cgd
76 1.1 cgd ! d does not fit in an int, so subtract 2^31, convert,
77 1.1 cgd ! and add 2^31 again (sigh). Just hope the intermediate
78 1.1 cgd ! fits (if not, the result is undefined anyway).
79 1.1 cgd
80 1.1 cgd fsubd %f0, %f2, %f0 ! d -= 2^31
81 1.1 cgd fdtoi %f0, %f0 ! convert to int
82 1.1 cgd st %f0, [%sp + 64] ! move into return reg
83 1.1 cgd ld [%sp + 64], %o0
84 1.1 cgd sethi %hi(0x80000000), %o1
85 1.1 cgd add %o0, %o1, %o0 ! add 2^31
86 1.1 cgd retl
87 1.1 cgd add %sp, 8, %sp
88 1.1 cgd
89 1.1 cgd 1:
90 1.1 cgd st %f0, [%sp + 64] ! return result
91 1.1 cgd ld [%sp + 64], %o0
92 1.1 cgd retl
93 1.1 cgd add %sp, 8, %sp
94