1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
|
dnl PPC-64 mpn_divrem_2 -- Divide an mpn number by a normalized 2-limb number.
dnl Copyright 2007, 2008 Free Software Foundation, Inc.
dnl This file is part of the GNU MP Library.
dnl
dnl The GNU MP Library is free software; you can redistribute it and/or modify
dnl it under the terms of either:
dnl
dnl * the GNU Lesser General Public License as published by the Free
dnl Software Foundation; either version 3 of the License, or (at your
dnl option) any later version.
dnl
dnl or
dnl
dnl * the GNU General Public License as published by the Free Software
dnl Foundation; either version 2 of the License, or (at your option) any
dnl later version.
dnl
dnl or both in parallel, as here.
dnl
dnl The GNU MP Library is distributed in the hope that it will be useful, but
dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
dnl for more details.
dnl
dnl You should have received copies of the GNU General Public License and the
dnl GNU Lesser General Public License along with the GNU MP Library. If not,
dnl see https://www.gnu.org/licenses/.
include(`../config.m4')
C cycles/limb
C norm frac
C POWER3/PPC630
C POWER4/PPC970 ? ?
C POWER5 37 ?
C POWER6 62 ?
C POWER6 30.5 ?
C INPUT PARAMETERS
C qp = r3
C fn = r4
C up = r5
C un = r6
C dp = r7
ifdef(`DARWIN',,`
define(`r2',`r31')') C FIXME!
ASM_START()
EXTERN_FUNC(mpn_invert_limb)
PROLOGUE(mpn_divrem_2,toc)
mflr r0
std r23, -72(r1)
std r24, -64(r1)
std r25, -56(r1)
std r26, -48(r1)
std r27, -40(r1)
std r28, -32(r1)
std r29, -24(r1)
std r30, -16(r1)
std r31, -8(r1)
std r0, 16(r1)
stdu r1, -192(r1)
mr r24, r3
mr r25, r4
sldi r0, r6, 3
add r26, r5, r0
addi r26, r26, -24
ld r30, 8(r7)
ld r28, 0(r7)
ld r29, 16(r26)
ld r31, 8(r26)
ifelse(0,1,`
li r23, 0
cmpld cr7, r29, r30
blt cr7, L(8)
bgt cr7, L(9)
cmpld cr0, r31, r28
blt cr0, L(8)
L(9): subfc r31, r28, r31
subfe r29, r30, r29
li r23, 1
',`
li r23, 0
cmpld cr7, r29, r30
blt cr7, L(8)
mfcr r0
rlwinm r0, r0, 30, 1
subfc r9, r28, r31
addze. r0, r0
nop
beq cr0, L(8)
subfc r31, r28, r31
subfe r29, r30, r29
li r23, 1
')
L(8):
add r27, r25, r6
addic. r27, r27, -3
blt cr0, L(18)
mr r3, r30
CALL( mpn_invert_limb)
mulld r10, r3, r30
mulhdu r0, r3, r28
addc r8, r10, r28
subfe r11, r1, r1
addc r10, r8, r0
addze. r11, r11
blt cr0, L(91)
L(40):
subfc r10, r30, r10
addme. r11, r11
addi r3, r3, -1
bge cr0, L(40)
L(91):
addi r5, r27, 1
mtctr r5
sldi r0, r27, 3
add r24, r24, r0
ALIGN(16)
L(loop):
mulhdu r8, r29, r3
mulld r6, r29, r3
addc r6, r6, r31
adde r8, r8, r29
cmpd cr7, r27, r25
mulld r0, r30, r8
mulhdu r11, r28, r8
mulld r10, r28, r8
subf r31, r0, r31
li r7, 0
blt cr7, L(60)
ld r7, 0(r26)
addi r26, r26, -8
nop
L(60): subfc r7, r28, r7
subfe r31, r30, r31
subfc r7, r10, r7
subfe r4, r11, r31
subfc r9, r6, r4
subfe r9, r1, r1
andc r6, r28, r9
andc r0, r30, r9
addc r31, r7, r6
adde r29, r4, r0
subf r8, r9, r8
cmpld cr7, r29, r30
bge- cr7, L(fix)
L(bck): std r8, 0(r24)
addi r24, r24, -8
addi r27, r27, -1
bdnz L(loop)
L(18):
std r31, 8(r26)
std r29, 16(r26)
mr r3, r23
addi r1, r1, 192
ld r0, 16(r1)
mtlr r0
ld r23, -72(r1)
ld r24, -64(r1)
ld r25, -56(r1)
ld r26, -48(r1)
ld r27, -40(r1)
ld r28, -32(r1)
ld r29, -24(r1)
ld r30, -16(r1)
ld r31, -8(r1)
blr
L(fix):
mfcr r0
rlwinm r0, r0, 30, 1
subfc r9, r28, r31
addze. r0, r0
beq cr0, L(bck)
subfc r31, r28, r31
subfe r29, r30, r29
addi r8, r8, 1
b L(bck)
EPILOGUE()
|
