• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
1 # Alpha 21064 __udiv_qrnnd
2 # Copyright (C) 1992, 1994, 1995, 2000 Free Software Foundation, Inc.
3
4 # This file is part of GCC.
5
6 # The GNU MP Library is free software; you can redistribute it and/or modify
7 # it under the terms of the GNU General Public License as published by
8 # the Free Software Foundation; either version 2 of the License, or (at your
9 # option) any later version.
10
11 # In addition to the permissions in the GNU General Public License, the
12 # Free Software Foundation gives you unlimited permission to link the
13 # compiled version of this file with other programs, and to distribute
14 # those programs without any restriction coming from the use of this
15 # file.  (The General Public License restrictions do apply in other
16 # respects; for example, they cover modification of the file, and
17 # distribution when not linked into another program.)
18
19 # This file is distributed in the hope that it will be useful, but
20 # WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
21 # or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public
22 # License for more details.
23
24 # You should have received a copy of the GNU General Public License
25 # along with GCC; see the file COPYING.  If not, write to the
26 # Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
27 # MA 02111-1307, USA.
28
29        .set noreorder
30        .set noat
31
32	.text
33
34	.globl __udiv_qrnnd
35	.ent __udiv_qrnnd
36__udiv_qrnnd:
37	.frame $30,0,$26,0
38	.prologue 0
39
40#define cnt	$2
41#define tmp	$3
42#define rem_ptr	$16
43#define n1	$17
44#define n0	$18
45#define d	$19
46#define qb	$20
47#define AT	$at
48
49	ldiq	cnt,16
50	blt	d,$largedivisor
51
52$loop1:	cmplt	n0,0,tmp
53	addq	n1,n1,n1
54	bis	n1,tmp,n1
55	addq	n0,n0,n0
56	cmpule	d,n1,qb
57	subq	n1,d,tmp
58	cmovne	qb,tmp,n1
59	bis	n0,qb,n0
60	cmplt	n0,0,tmp
61	addq	n1,n1,n1
62	bis	n1,tmp,n1
63	addq	n0,n0,n0
64	cmpule	d,n1,qb
65	subq	n1,d,tmp
66	cmovne	qb,tmp,n1
67	bis	n0,qb,n0
68	cmplt	n0,0,tmp
69	addq	n1,n1,n1
70	bis	n1,tmp,n1
71	addq	n0,n0,n0
72	cmpule	d,n1,qb
73	subq	n1,d,tmp
74	cmovne	qb,tmp,n1
75	bis	n0,qb,n0
76	cmplt	n0,0,tmp
77	addq	n1,n1,n1
78	bis	n1,tmp,n1
79	addq	n0,n0,n0
80	cmpule	d,n1,qb
81	subq	n1,d,tmp
82	cmovne	qb,tmp,n1
83	bis	n0,qb,n0
84	subq	cnt,1,cnt
85	bgt	cnt,$loop1
86	stq	n1,0(rem_ptr)
87	bis	$31,n0,$0
88	ret	$31,($26),1
89
90$largedivisor:
91	and	n0,1,$4
92
93	srl	n0,1,n0
94	sll	n1,63,tmp
95	or	tmp,n0,n0
96	srl	n1,1,n1
97
98	and	d,1,$6
99	srl	d,1,$5
100	addq	$5,$6,$5
101
102$loop2:	cmplt	n0,0,tmp
103	addq	n1,n1,n1
104	bis	n1,tmp,n1
105	addq	n0,n0,n0
106	cmpule	$5,n1,qb
107	subq	n1,$5,tmp
108	cmovne	qb,tmp,n1
109	bis	n0,qb,n0
110	cmplt	n0,0,tmp
111	addq	n1,n1,n1
112	bis	n1,tmp,n1
113	addq	n0,n0,n0
114	cmpule	$5,n1,qb
115	subq	n1,$5,tmp
116	cmovne	qb,tmp,n1
117	bis	n0,qb,n0
118	cmplt	n0,0,tmp
119	addq	n1,n1,n1
120	bis	n1,tmp,n1
121	addq	n0,n0,n0
122	cmpule	$5,n1,qb
123	subq	n1,$5,tmp
124	cmovne	qb,tmp,n1
125	bis	n0,qb,n0
126	cmplt	n0,0,tmp
127	addq	n1,n1,n1
128	bis	n1,tmp,n1
129	addq	n0,n0,n0
130	cmpule	$5,n1,qb
131	subq	n1,$5,tmp
132	cmovne	qb,tmp,n1
133	bis	n0,qb,n0
134	subq	cnt,1,cnt
135	bgt	cnt,$loop2
136
137	addq	n1,n1,n1
138	addq	$4,n1,n1
139	bne	$6,$Odd
140	stq	n1,0(rem_ptr)
141	bis	$31,n0,$0
142	ret	$31,($26),1
143
144$Odd:
145	/* q' in n0. r' in n1 */
146	addq	n1,n0,n1
147
148	cmpult	n1,n0,tmp	# tmp := carry from addq
149	subq	n1,d,AT
150	addq	n0,tmp,n0
151	cmovne	tmp,AT,n1
152
153	cmpult	n1,d,tmp
154	addq	n0,1,AT
155	cmoveq	tmp,AT,n0
156	subq	n1,d,AT
157	cmoveq	tmp,AT,n1
158
159	stq	n1,0(rem_ptr)
160	bis	$31,n0,$0
161	ret	$31,($26),1
162
163	.end	__udiv_qrnnd
164