|
Packit |
5c3484 |
dnl ARM64 mpn_mod_34lsub1 -- remainder modulo 2^48-1.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Copyright 2012-2014 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
|
|
Packit |
5c3484 |
C Cortex-A53 ?
|
|
Packit |
5c3484 |
C Cortex-A57 ?
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`ap', x0)
|
|
Packit |
5c3484 |
define(`n', x1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
changecom(@&*$)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C mp_limb_t mpn_mod_34lsub1 (mp_srcptr up, mp_size_t n)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C TODO
|
|
Packit |
5c3484 |
C * An alternative inner loop which could run at 0.722 c/l:
|
|
Packit |
5c3484 |
C adds x8, x8, x2
|
|
Packit |
5c3484 |
C adcs x9, x9, x3
|
|
Packit |
5c3484 |
C ldp x2, x3, [ap, #-32]
|
|
Packit |
5c3484 |
C adcs x10, x10, x4
|
|
Packit |
5c3484 |
C adc x12, x12, xzr
|
|
Packit |
5c3484 |
C adds x8, x8, x5
|
|
Packit |
5c3484 |
C ldp x4, x5, [ap, #-16]
|
|
Packit |
5c3484 |
C sub n, n, #6
|
|
Packit |
5c3484 |
C adcs x9, x9, x6
|
|
Packit |
5c3484 |
C adcs x10, x10, x7
|
|
Packit |
5c3484 |
C ldp x6, x7, [ap], #48
|
|
Packit |
5c3484 |
C adc x12, x12, xzr
|
|
Packit |
5c3484 |
C tbz n, #63, L(top)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ASM_START()
|
|
Packit |
5c3484 |
TEXT
|
|
Packit |
5c3484 |
ALIGN(32)
|
|
Packit |
5c3484 |
PROLOGUE(mpn_mod_34lsub1)
|
|
Packit |
5c3484 |
subs n, n, #3
|
|
Packit |
5c3484 |
mov x8, #0
|
|
Packit |
5c3484 |
b.lt L(le2) C n <= 2
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ldp x2, x3, [ap, #0]
|
|
Packit |
5c3484 |
ldr x4, [ap, #16]
|
|
Packit |
5c3484 |
add ap, ap, #24
|
|
Packit |
5c3484 |
subs n, n, #3
|
|
Packit |
5c3484 |
b.lt L(sum) C n <= 5
|
|
Packit |
5c3484 |
cmn x0, #0 C clear carry
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(top): ldp x5, x6, [ap, #0]
|
|
Packit |
5c3484 |
ldr x7, [ap, #16]
|
|
Packit |
5c3484 |
add ap, ap, #24
|
|
Packit |
5c3484 |
sub n, n, #3
|
|
Packit |
5c3484 |
adcs x2, x2, x5
|
|
Packit |
5c3484 |
adcs x3, x3, x6
|
|
Packit |
5c3484 |
adcs x4, x4, x7
|
|
Packit |
5c3484 |
tbz n, #63, L(top)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
adc x8, xzr, xzr C x8 <= 1
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(sum): cmn n, #2
|
|
Packit |
5c3484 |
mov x5, #0
|
|
Packit |
5c3484 |
b.lo 1f
|
|
Packit |
5c3484 |
ldr x5, [ap], #8
|
|
Packit |
5c3484 |
1: mov x6, #0
|
|
Packit |
5c3484 |
b.ls 1f
|
|
Packit |
5c3484 |
ldr x6, [ap], #8
|
|
Packit |
5c3484 |
1: adds x2, x2, x5
|
|
Packit |
5c3484 |
adcs x3, x3, x6
|
|
Packit |
5c3484 |
adcs x4, x4, xzr
|
|
Packit |
5c3484 |
adc x8, x8, xzr C x8 <= 2
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(sum2):
|
|
Packit |
5c3484 |
and x0, x2, #0xffffffffffff
|
|
Packit |
5c3484 |
add x0, x0, x2, lsr #48
|
|
Packit |
5c3484 |
add x0, x0, x8
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lsl x8, x3, #16
|
|
Packit |
5c3484 |
and x1, x8, #0xffffffffffff
|
|
Packit |
5c3484 |
add x0, x0, x1
|
|
Packit |
5c3484 |
add x0, x0, x3, lsr #32
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lsl x8, x4, #32
|
|
Packit |
5c3484 |
and x1, x8, #0xffffffffffff
|
|
Packit |
5c3484 |
add x0, x0, x1
|
|
Packit |
5c3484 |
add x0, x0, x4, lsr #16
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(le2): cmn n, #1
|
|
Packit |
5c3484 |
b.ne L(1)
|
|
Packit |
5c3484 |
ldp x2, x3, [ap]
|
|
Packit |
5c3484 |
mov x4, #0
|
|
Packit |
5c3484 |
b L(sum2)
|
|
Packit |
5c3484 |
L(1): ldr x2, [ap]
|
|
Packit |
5c3484 |
and x0, x2, #0xffffffffffff
|
|
Packit |
5c3484 |
add x0, x0, x2, lsr #48
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
EPILOGUE()
|