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