|
Packit |
5c3484 |
dnl ARM v6 mpn_divexact_1
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Contributed to the GNU project by Torbjörn Granlund.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Copyright 2012, 2013 Free Software Foundation, Inc.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl This file is part of the GNU MP Library.
|
|
Packit |
5c3484 |
dnl
|
|
Packit |
5c3484 |
dnl The GNU MP Library is free software; you can redistribute it and/or modify
|
|
Packit |
5c3484 |
dnl it under the terms of either:
|
|
Packit |
5c3484 |
dnl
|
|
Packit |
5c3484 |
dnl * the GNU Lesser General Public License as published by the Free
|
|
Packit |
5c3484 |
dnl Software Foundation; either version 3 of the License, or (at your
|
|
Packit |
5c3484 |
dnl option) any later version.
|
|
Packit |
5c3484 |
dnl
|
|
Packit |
5c3484 |
dnl or
|
|
Packit |
5c3484 |
dnl
|
|
Packit |
5c3484 |
dnl * the GNU General Public License as published by the Free Software
|
|
Packit |
5c3484 |
dnl Foundation; either version 2 of the License, or (at your option) any
|
|
Packit |
5c3484 |
dnl later version.
|
|
Packit |
5c3484 |
dnl
|
|
Packit |
5c3484 |
dnl or both in parallel, as here.
|
|
Packit |
5c3484 |
dnl
|
|
Packit |
5c3484 |
dnl The GNU MP Library is distributed in the hope that it will be useful, but
|
|
Packit |
5c3484 |
dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
|
|
Packit |
5c3484 |
dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
|
Packit |
5c3484 |
dnl for more details.
|
|
Packit |
5c3484 |
dnl
|
|
Packit |
5c3484 |
dnl You should have received copies of the GNU General Public License and the
|
|
Packit |
5c3484 |
dnl GNU Lesser General Public License along with the GNU MP Library. If not,
|
|
Packit |
5c3484 |
dnl see https://www.gnu.org/licenses/.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
include(`../config.m4')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C cycles/limb cycles/limb
|
|
Packit |
5c3484 |
C norm unorm modexact_1c_odd
|
|
Packit |
5c3484 |
C StrongARM - -
|
|
Packit |
5c3484 |
C XScale - -
|
|
Packit |
5c3484 |
C Cortex-A7 ? ?
|
|
Packit |
5c3484 |
C Cortex-A8 ? ?
|
|
Packit |
5c3484 |
C Cortex-A9 9 10 9
|
|
Packit |
5c3484 |
C Cortex-A15 7 7 7
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C Architecture requirements:
|
|
Packit |
5c3484 |
C v5 -
|
|
Packit |
5c3484 |
C v5t clz
|
|
Packit |
5c3484 |
C v5te -
|
|
Packit |
5c3484 |
C v6 umaal
|
|
Packit |
5c3484 |
C v6t2 -
|
|
Packit |
5c3484 |
C v7a -
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`rp', `r0')
|
|
Packit |
5c3484 |
define(`up', `r1')
|
|
Packit |
5c3484 |
define(`n', `r2')
|
|
Packit |
5c3484 |
define(`d', `r3')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`cy', `r7')
|
|
Packit |
5c3484 |
define(`cnt', `r6')
|
|
Packit |
5c3484 |
define(`tnc', `r10')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ASM_START()
|
|
Packit |
5c3484 |
PROLOGUE(mpn_divexact_1)
|
|
Packit |
5c3484 |
push {r4,r5,r6,r7,r8,r9}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
tst d, #1
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
rsb r4, d, #0
|
|
Packit |
5c3484 |
and r4, r4, d
|
|
Packit |
5c3484 |
clz r4, r4
|
|
Packit |
5c3484 |
rsb cnt, r4, #31 C count_trailing_zeros
|
|
Packit |
5c3484 |
mov d, d, lsr cnt
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C binvert limb
|
|
Packit |
5c3484 |
LEA( r4, binvert_limb_table)
|
|
Packit |
5c3484 |
and r12, d, #254
|
|
Packit |
5c3484 |
ldrb r4, [r4, r12, lsr #1]
|
|
Packit |
5c3484 |
mul r12, r4, r4
|
|
Packit |
5c3484 |
mul r12, d, r12
|
|
Packit |
5c3484 |
rsb r12, r12, r4, lsl #1
|
|
Packit |
5c3484 |
mul r4, r12, r12
|
|
Packit |
5c3484 |
mul r4, d, r4
|
|
Packit |
5c3484 |
rsb r4, r4, r12, lsl #1 C r4 = inverse
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ldr r5, [up], #4 C up[0]
|
|
Packit |
5c3484 |
mov cy, #0
|
|
Packit |
5c3484 |
rsb r8, r4, #0 C r8 = -inverse
|
|
Packit |
5c3484 |
beq L(unnorm)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(norm):
|
|
Packit |
5c3484 |
subs n, n, #1
|
|
Packit |
5c3484 |
mul r5, r5, r4
|
|
Packit |
5c3484 |
beq L(end)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(top): ldr r9, [up], #4
|
|
Packit |
5c3484 |
mov r12, #0
|
|
Packit |
5c3484 |
str r5, [rp], #4
|
|
Packit |
5c3484 |
umaal r12, cy, r5, d
|
|
Packit |
5c3484 |
mul r5, r9, r4
|
|
Packit |
5c3484 |
mla r5, cy, r8, r5
|
|
Packit |
5c3484 |
subs n, n, #1
|
|
Packit |
5c3484 |
bne L(top)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(end): str r5, [rp]
|
|
Packit |
5c3484 |
pop {r4,r5,r6,r7,r8,r9}
|
|
Packit |
5c3484 |
bx r14
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(unnorm):
|
|
Packit |
5c3484 |
push {r10,r11}
|
|
Packit |
5c3484 |
rsb tnc, cnt, #32
|
|
Packit |
5c3484 |
mov r11, r5, lsr cnt
|
|
Packit |
5c3484 |
subs n, n, #1
|
|
Packit |
5c3484 |
beq L(edx)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ldr r12, [up], #4
|
|
Packit |
5c3484 |
orr r9, r11, r12, lsl tnc
|
|
Packit |
5c3484 |
mov r11, r12, lsr cnt
|
|
Packit |
5c3484 |
mul r5, r9, r4
|
|
Packit |
5c3484 |
subs n, n, #1
|
|
Packit |
5c3484 |
beq L(edu)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(tpu): ldr r12, [up], #4
|
|
Packit |
5c3484 |
orr r9, r11, r12, lsl tnc
|
|
Packit |
5c3484 |
mov r11, r12, lsr cnt
|
|
Packit |
5c3484 |
mov r12, #0
|
|
Packit |
5c3484 |
str r5, [rp], #4
|
|
Packit |
5c3484 |
umaal r12, cy, r5, d
|
|
Packit |
5c3484 |
mul r5, r9, r4
|
|
Packit |
5c3484 |
mla r5, cy, r8, r5
|
|
Packit |
5c3484 |
subs n, n, #1
|
|
Packit |
5c3484 |
bne L(tpu)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(edu): str r5, [rp], #4
|
|
Packit |
5c3484 |
mov r12, #0
|
|
Packit |
5c3484 |
umaal r12, cy, r5, d
|
|
Packit |
5c3484 |
mul r5, r11, r4
|
|
Packit |
5c3484 |
mla r5, cy, r8, r5
|
|
Packit |
5c3484 |
str r5, [rp]
|
|
Packit |
5c3484 |
pop {r10,r11}
|
|
Packit |
5c3484 |
pop {r4,r5,r6,r7,r8,r9}
|
|
Packit |
5c3484 |
bx r14
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(edx): mul r5, r11, r4
|
|
Packit |
5c3484 |
str r5, [rp]
|
|
Packit |
5c3484 |
pop {r10,r11}
|
|
Packit |
5c3484 |
pop {r4,r5,r6,r7,r8,r9}
|
|
Packit |
5c3484 |
bx r14
|
|
Packit |
5c3484 |
EPILOGUE()
|