|
Packit |
5c3484 |
dnl AMD64 mpn_rsh1add_n -- rp[] = (up[] + vp[]) >> 1
|
|
Packit |
5c3484 |
dnl AMD64 mpn_rsh1sub_n -- rp[] = (up[] - vp[]) >> 1
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Copyright 2003, 2005, 2009, 2011, 2012 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 AMD K8,K9 2.14 (mpn_add_n + mpn_rshift need 4.125)
|
|
Packit |
5c3484 |
C AMD K10 2.14 (mpn_add_n + mpn_rshift need 4.125)
|
|
Packit |
5c3484 |
C Intel P4 12.75
|
|
Packit |
5c3484 |
C Intel core2 3.75
|
|
Packit |
5c3484 |
C Intel NMH 4.4
|
|
Packit |
5c3484 |
C Intel SBR ?
|
|
Packit |
5c3484 |
C Intel atom ?
|
|
Packit |
5c3484 |
C VIA nano 3.25
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C TODO
|
|
Packit |
5c3484 |
C * Rewrite to use indexed addressing, like addlsh1.asm and sublsh1.asm.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C INPUT PARAMETERS
|
|
Packit |
5c3484 |
define(`rp', `%rdi')
|
|
Packit |
5c3484 |
define(`up', `%rsi')
|
|
Packit |
5c3484 |
define(`vp', `%rdx')
|
|
Packit |
5c3484 |
define(`n',` %rcx')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ifdef(`OPERATION_rsh1add_n', `
|
|
Packit |
5c3484 |
define(ADDSUB, add)
|
|
Packit |
5c3484 |
define(ADCSBB, adc)
|
|
Packit |
5c3484 |
define(func_n, mpn_rsh1add_n)
|
|
Packit |
5c3484 |
define(func_nc, mpn_rsh1add_nc)')
|
|
Packit |
5c3484 |
ifdef(`OPERATION_rsh1sub_n', `
|
|
Packit |
5c3484 |
define(ADDSUB, sub)
|
|
Packit |
5c3484 |
define(ADCSBB, sbb)
|
|
Packit |
5c3484 |
define(func_n, mpn_rsh1sub_n)
|
|
Packit |
5c3484 |
define(func_nc, mpn_rsh1sub_nc)')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
MULFUNC_PROLOGUE(mpn_rsh1add_n mpn_rsh1add_nc mpn_rsh1sub_n mpn_rsh1sub_nc)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ABI_SUPPORT(DOS64)
|
|
Packit |
5c3484 |
ABI_SUPPORT(STD64)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ASM_START()
|
|
Packit |
5c3484 |
TEXT
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
PROLOGUE(func_nc)
|
|
Packit |
5c3484 |
FUNC_ENTRY(4)
|
|
Packit |
5c3484 |
IFDOS(` mov 56(%rsp), %r8 ')
|
|
Packit |
5c3484 |
push %rbx
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
xor R32(%rax), R32(%rax)
|
|
Packit |
5c3484 |
neg %r8 C set C flag from parameter
|
|
Packit |
5c3484 |
mov (up), %rbx
|
|
Packit |
5c3484 |
ADCSBB (vp), %rbx
|
|
Packit |
5c3484 |
jmp L(ent)
|
|
Packit |
5c3484 |
EPILOGUE()
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
PROLOGUE(func_n)
|
|
Packit |
5c3484 |
FUNC_ENTRY(4)
|
|
Packit |
5c3484 |
push %rbx
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
xor R32(%rax), R32(%rax)
|
|
Packit |
5c3484 |
mov (up), %rbx
|
|
Packit |
5c3484 |
ADDSUB (vp), %rbx
|
|
Packit |
5c3484 |
L(ent):
|
|
Packit |
5c3484 |
rcr %rbx C rotate, save acy
|
|
Packit |
5c3484 |
adc R32(%rax), R32(%rax) C return value
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov R32(n), R32(%r11)
|
|
Packit |
5c3484 |
and $3, R32(%r11)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
cmp $1, R32(%r11)
|
|
Packit |
5c3484 |
je L(do) C jump if n = 1 5 9 ...
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(n1): cmp $2, R32(%r11)
|
|
Packit |
5c3484 |
jne L(n2) C jump unless n = 2 6 10 ...
|
|
Packit |
5c3484 |
add %rbx, %rbx C rotate carry limb, restore acy
|
|
Packit |
5c3484 |
mov 8(up), %r10
|
|
Packit |
5c3484 |
ADCSBB 8(vp), %r10
|
|
Packit |
5c3484 |
lea 8(up), up
|
|
Packit |
5c3484 |
lea 8(vp), vp
|
|
Packit |
5c3484 |
lea 8(rp), rp
|
|
Packit |
5c3484 |
rcr %r10
|
|
Packit |
5c3484 |
rcr %rbx
|
|
Packit |
5c3484 |
mov %rbx, -8(rp)
|
|
Packit |
5c3484 |
jmp L(cj1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(n2): cmp $3, R32(%r11)
|
|
Packit |
5c3484 |
jne L(n3) C jump unless n = 3 7 11 ...
|
|
Packit |
5c3484 |
add %rbx, %rbx C rotate carry limb, restore acy
|
|
Packit |
5c3484 |
mov 8(up), %r9
|
|
Packit |
5c3484 |
mov 16(up), %r10
|
|
Packit |
5c3484 |
ADCSBB 8(vp), %r9
|
|
Packit |
5c3484 |
ADCSBB 16(vp), %r10
|
|
Packit |
5c3484 |
lea 16(up), up
|
|
Packit |
5c3484 |
lea 16(vp), vp
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
rcr %r10
|
|
Packit |
5c3484 |
rcr %r9
|
|
Packit |
5c3484 |
rcr %rbx
|
|
Packit |
5c3484 |
mov %rbx, -16(rp)
|
|
Packit |
5c3484 |
jmp L(cj2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(n3): dec n C come here for n = 4 8 12 ...
|
|
Packit |
5c3484 |
add %rbx, %rbx C rotate carry limb, restore acy
|
|
Packit |
5c3484 |
mov 8(up), %r8
|
|
Packit |
5c3484 |
mov 16(up), %r9
|
|
Packit |
5c3484 |
ADCSBB 8(vp), %r8
|
|
Packit |
5c3484 |
ADCSBB 16(vp), %r9
|
|
Packit |
5c3484 |
mov 24(up), %r10
|
|
Packit |
5c3484 |
ADCSBB 24(vp), %r10
|
|
Packit |
5c3484 |
lea 24(up), up
|
|
Packit |
5c3484 |
lea 24(vp), vp
|
|
Packit |
5c3484 |
lea 24(rp), rp
|
|
Packit |
5c3484 |
rcr %r10
|
|
Packit |
5c3484 |
rcr %r9
|
|
Packit |
5c3484 |
rcr %r8
|
|
Packit |
5c3484 |
rcr %rbx
|
|
Packit |
5c3484 |
mov %rbx, -24(rp)
|
|
Packit |
5c3484 |
mov %r8, -16(rp)
|
|
Packit |
5c3484 |
L(cj2): mov %r9, -8(rp)
|
|
Packit |
5c3484 |
L(cj1): mov %r10, %rbx
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(do):
|
|
Packit |
5c3484 |
shr $2, n C 4
|
|
Packit |
5c3484 |
je L(end) C 2
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(top): add %rbx, %rbx C rotate carry limb, restore acy
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov 8(up), %r8
|
|
Packit |
5c3484 |
mov 16(up), %r9
|
|
Packit |
5c3484 |
ADCSBB 8(vp), %r8
|
|
Packit |
5c3484 |
ADCSBB 16(vp), %r9
|
|
Packit |
5c3484 |
mov 24(up), %r10
|
|
Packit |
5c3484 |
mov 32(up), %r11
|
|
Packit |
5c3484 |
ADCSBB 24(vp), %r10
|
|
Packit |
5c3484 |
ADCSBB 32(vp), %r11
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 32(up), up
|
|
Packit |
5c3484 |
lea 32(vp), vp
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
rcr %r11 C rotate, save acy
|
|
Packit |
5c3484 |
rcr %r10
|
|
Packit |
5c3484 |
rcr %r9
|
|
Packit |
5c3484 |
rcr %r8
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
rcr %rbx
|
|
Packit |
5c3484 |
mov %rbx, (rp)
|
|
Packit |
5c3484 |
mov %r8, 8(rp)
|
|
Packit |
5c3484 |
mov %r9, 16(rp)
|
|
Packit |
5c3484 |
mov %r10, 24(rp)
|
|
Packit |
5c3484 |
mov %r11, %rbx
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 32(rp), rp
|
|
Packit |
5c3484 |
dec n
|
|
Packit |
5c3484 |
jne L(top)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(end): mov %rbx, (rp)
|
|
Packit |
5c3484 |
pop %rbx
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
EPILOGUE()
|