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