![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to udiv_qrnnd.S CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / mpn / alpha|
Return to udiv_qrnnd.S CVS log | Up to [local] / OpenXM_contrib / gmp / mpn / alpha |
1.1 maekawa 1: # Alpha 21064 __udiv_qrnnd
2:
1.1.1.2 ! maekawa 3: # Copyright (C) 1992, 1994, 1995, 1997, 2000 Free Software Foundation, Inc.
1.1 maekawa 4:
5: # This file is part of the GNU MP Library.
6:
7: # The GNU MP Library is free software; you can redistribute it and/or modify
1.1.1.2 ! maekawa 8: # it under the terms of the GNU Lesser General Public License as published by
! 9: # the Free Software Foundation; either version 2.1 of the License, or (at your
1.1 maekawa 10: # option) any later version.
11:
12: # The GNU MP Library is distributed in the hope that it will be useful, but
13: # WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
1.1.1.2 ! maekawa 14: # or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
1.1 maekawa 15: # License for more details.
16:
1.1.1.2 ! maekawa 17: # You should have received a copy of the GNU Lesser General Public License
1.1 maekawa 18: # along with the GNU MP Library; see the file COPYING.LIB. If not, write to
19: # the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
20: # MA 02111-1307, USA.
21:
22:
1.1.1.2 ! maekawa 23: .set noreorder
! 24: .set noat
1.1 maekawa 25: .text
1.1.1.2 ! maekawa 26: .align 3
! 27: .globl __gmpn_udiv_qrnnd
! 28: .ent __gmpn_udiv_qrnnd
! 29: __gmpn_udiv_qrnnd:
! 30: .frame $30,0,$26,0
! 31: .prologue 0
1.1 maekawa 32: #define cnt $2
33: #define tmp $3
34: #define rem_ptr $16
35: #define n1 $17
36: #define n0 $18
37: #define d $19
38: #define qb $20
39:
40: ldiq cnt,16
41: blt d,.Largedivisor
42:
43: .Loop1: cmplt n0,0,tmp
44: addq n1,n1,n1
45: bis n1,tmp,n1
46: addq n0,n0,n0
47: cmpule d,n1,qb
48: subq n1,d,tmp
49: cmovne qb,tmp,n1
50: bis n0,qb,n0
51: cmplt n0,0,tmp
52: addq n1,n1,n1
53: bis n1,tmp,n1
54: addq n0,n0,n0
55: cmpule d,n1,qb
56: subq n1,d,tmp
57: cmovne qb,tmp,n1
58: bis n0,qb,n0
59: cmplt n0,0,tmp
60: addq n1,n1,n1
61: bis n1,tmp,n1
62: addq n0,n0,n0
63: cmpule d,n1,qb
64: subq n1,d,tmp
65: cmovne qb,tmp,n1
66: bis n0,qb,n0
67: cmplt n0,0,tmp
68: addq n1,n1,n1
69: bis n1,tmp,n1
70: addq n0,n0,n0
71: cmpule d,n1,qb
72: subq n1,d,tmp
73: cmovne qb,tmp,n1
74: bis n0,qb,n0
75: subq cnt,1,cnt
76: bgt cnt,.Loop1
77: stq n1,0(rem_ptr)
78: bis $31,n0,$0
79: ret $31,($26),1
80:
81: .Largedivisor:
82: and n0,1,$4
83:
84: srl n0,1,n0
85: sll n1,63,tmp
86: or tmp,n0,n0
87: srl n1,1,n1
88:
89: and d,1,$6
90: srl d,1,$5
91: addq $5,$6,$5
92:
93: .Loop2: cmplt n0,0,tmp
94: addq n1,n1,n1
95: bis n1,tmp,n1
96: addq n0,n0,n0
97: cmpule $5,n1,qb
98: subq n1,$5,tmp
99: cmovne qb,tmp,n1
100: bis n0,qb,n0
101: cmplt n0,0,tmp
102: addq n1,n1,n1
103: bis n1,tmp,n1
104: addq n0,n0,n0
105: cmpule $5,n1,qb
106: subq n1,$5,tmp
107: cmovne qb,tmp,n1
108: bis n0,qb,n0
109: cmplt n0,0,tmp
110: addq n1,n1,n1
111: bis n1,tmp,n1
112: addq n0,n0,n0
113: cmpule $5,n1,qb
114: subq n1,$5,tmp
115: cmovne qb,tmp,n1
116: bis n0,qb,n0
117: cmplt n0,0,tmp
118: addq n1,n1,n1
119: bis n1,tmp,n1
120: addq n0,n0,n0
121: cmpule $5,n1,qb
122: subq n1,$5,tmp
123: cmovne qb,tmp,n1
124: bis n0,qb,n0
125: subq cnt,1,cnt
126: bgt cnt,.Loop2
127:
128: addq n1,n1,n1
129: addq $4,n1,n1
130: bne $6,.LOdd
131: stq n1,0(rem_ptr)
132: bis $31,n0,$0
133: ret $31,($26),1
134:
135: .LOdd:
136: /* q' in n0. r' in n1 */
137: addq n1,n0,n1
138: cmpult n1,n0,tmp # tmp := carry from addq
139: beq tmp,.LLp6
140: addq n0,1,n0
141: subq n1,d,n1
142: .LLp6: cmpult n1,d,tmp
143: bne tmp,.LLp7
144: addq n0,1,n0
145: subq n1,d,n1
146: .LLp7:
147: stq n1,0(rem_ptr)
148: bis $31,n0,$0
149: ret $31,($26),1
150:
1.1.1.2 ! maekawa 151: .end __gmpn_udiv_qrnnd