|
Packit |
5c3484 |
dnl AMD64 mpn_mullo_basecase optimised for Intel Sandy bridge and Ivy bridge.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Contributed to the GNU project by Torbjörn Granlund.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Copyright 2008, 2009, 2011-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 mul_2 addmul_2
|
|
Packit |
5c3484 |
C AMD K8,K9
|
|
Packit |
5c3484 |
C AMD K10
|
|
Packit |
5c3484 |
C AMD bull
|
|
Packit |
5c3484 |
C AMD pile
|
|
Packit |
5c3484 |
C AMD steam
|
|
Packit |
5c3484 |
C AMD bobcat
|
|
Packit |
5c3484 |
C AMD jaguar
|
|
Packit |
5c3484 |
C Intel P4
|
|
Packit |
5c3484 |
C Intel core
|
|
Packit |
5c3484 |
C Intel NHM
|
|
Packit |
5c3484 |
C Intel SBR 2.5 2.95
|
|
Packit |
5c3484 |
C Intel IBR 2.3 2.68
|
|
Packit |
5c3484 |
C Intel HWL 2.0 2.5
|
|
Packit |
5c3484 |
C Intel BWL
|
|
Packit |
5c3484 |
C Intel atom
|
|
Packit |
5c3484 |
C VIA nano
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C The inner loops of this code are the result of running a code generation and
|
|
Packit |
5c3484 |
C optimisation tool suite written by David Harvey and Torbjörn Granlund.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C TODO
|
|
Packit |
5c3484 |
C * Implement proper cor2, replacing current cor0.
|
|
Packit |
5c3484 |
C * Offset n by 2 in order to avoid the outer loop cmp. (And sqr_basecase?)
|
|
Packit |
5c3484 |
C * Micro-optimise.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C When playing with pointers, set this to $2 to fall back to conservative
|
|
Packit |
5c3484 |
C indexing in wind-down code.
|
|
Packit |
5c3484 |
define(`I',`$1')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`rp', `%rdi')
|
|
Packit |
5c3484 |
define(`up', `%rsi')
|
|
Packit |
5c3484 |
define(`vp_param', `%rdx')
|
|
Packit |
5c3484 |
define(`n', `%rcx')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`vp', `%r8')
|
|
Packit |
5c3484 |
define(`X0', `%r14')
|
|
Packit |
5c3484 |
define(`X1', `%r15')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`w0', `%r10')
|
|
Packit |
5c3484 |
define(`w1', `%r11')
|
|
Packit |
5c3484 |
define(`w2', `%r12')
|
|
Packit |
5c3484 |
define(`w3', `%r13')
|
|
Packit |
5c3484 |
define(`i', `%rbp')
|
|
Packit |
5c3484 |
define(`v0', `%r9')
|
|
Packit |
5c3484 |
define(`v1', `%rbx')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C rax rbx rcx rdx rdi rsi rbp r8 r9 r10 r11 r12 r13 r14 r15
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ABI_SUPPORT(DOS64)
|
|
Packit |
5c3484 |
ABI_SUPPORT(STD64)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ASM_START()
|
|
Packit |
5c3484 |
TEXT
|
|
Packit |
5c3484 |
ALIGN(32)
|
|
Packit |
5c3484 |
PROLOGUE(mpn_mullo_basecase)
|
|
Packit |
5c3484 |
FUNC_ENTRY(4)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov (up), %rax
|
|
Packit |
5c3484 |
mov vp_param, vp
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
cmp $4, n
|
|
Packit |
5c3484 |
jb L(small)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov (vp_param), v0
|
|
Packit |
5c3484 |
push %rbx
|
|
Packit |
5c3484 |
lea (rp,n,8), rp C point rp at R[un]
|
|
Packit |
5c3484 |
push %rbp
|
|
Packit |
5c3484 |
lea (up,n,8), up C point up right after U's end
|
|
Packit |
5c3484 |
push %r12
|
|
Packit |
5c3484 |
neg n
|
|
Packit |
5c3484 |
push %r13
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov 8(vp), v1
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
test $1, R8(n)
|
|
Packit |
5c3484 |
jnz L(m2b1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(m2b0):lea (n), i
|
|
Packit |
5c3484 |
xor w0, w0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
jmp L(m2l0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(m2b1):lea 1(n), i
|
|
Packit |
5c3484 |
xor w1, w1
|
|
Packit |
5c3484 |
xor w2, w2
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
jmp L(m2l1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(32)
|
|
Packit |
5c3484 |
L(m2tp):mul v0
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
L(m2l1):mov -8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w1, w0
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
mov w0, -8(rp,i,8)
|
|
Packit |
5c3484 |
mov %rdx, w0
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
add w3, w2
|
|
Packit |
5c3484 |
L(m2l0):mov (up,i,8), %rax
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
mov w2, (rp,i,8)
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w2 C FIXME: dead in last iteration
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
adc $0, w2 C FIXME: dead in last iteration
|
|
Packit |
5c3484 |
add $2, i
|
|
Packit |
5c3484 |
jnc L(m2tp)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(m2ed):imul v0, %rax
|
|
Packit |
5c3484 |
add w0, %rax
|
|
Packit |
5c3484 |
add w1, %rax
|
|
Packit |
5c3484 |
mov %rax, I(-8(rp),-8(rp,i,8))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
add $2, n
|
|
Packit |
5c3484 |
lea 16(vp), vp
|
|
Packit |
5c3484 |
lea -16(up), up
|
|
Packit |
5c3484 |
cmp $-2, n
|
|
Packit |
5c3484 |
jge L(cor1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
push %r14
|
|
Packit |
5c3484 |
push %r15
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(outer):
|
|
Packit |
5c3484 |
mov (vp), v0
|
|
Packit |
5c3484 |
mov 8(vp), v1
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
test $1, R8(n)
|
|
Packit |
5c3484 |
jnz L(a1x1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(a1x0):mov (rp,n,8), X1
|
|
Packit |
5c3484 |
xor w2, w2
|
|
Packit |
5c3484 |
xor w1, w1
|
|
Packit |
5c3484 |
test $2, R8(n)
|
|
Packit |
5c3484 |
jnz L(a110)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(a100):lea 1(n), i
|
|
Packit |
5c3484 |
jmp L(lo0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(a110):lea 3(n), i
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add %rax, X1
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mov 8(rp,n,8), X0
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
jmp L(lo2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(a1x1):mov (rp,n,8), X0
|
|
Packit |
5c3484 |
xor w0, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
test $2, R8(n)
|
|
Packit |
5c3484 |
jz L(a111)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(a101):lea 2(n), i
|
|
Packit |
5c3484 |
add %rax, X0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
mov 8(rp,n,8), X1
|
|
Packit |
5c3484 |
jmp L(lo1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(a111):lea (n), i
|
|
Packit |
5c3484 |
xor w3, w3
|
|
Packit |
5c3484 |
jmp L(lo3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(32)
|
|
Packit |
5c3484 |
L(top):
|
|
Packit |
5c3484 |
L(lo2): mul v1
|
|
Packit |
5c3484 |
mov %rdx, w0
|
|
Packit |
5c3484 |
add %rax, X0
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add w2, X0
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
mov -16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, X0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov -16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
mov X1, -24(rp,i,8)
|
|
Packit |
5c3484 |
mov -8(rp,i,8), X1
|
|
Packit |
5c3484 |
add w3, X0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
L(lo1): mov %rdx, w2
|
|
Packit |
5c3484 |
mov X0, -16(rp,i,8)
|
|
Packit |
5c3484 |
add %rax, X1
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
mov -8(up,i,8), %rax
|
|
Packit |
5c3484 |
add w0, X1
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
L(lo0): add %rax, X1
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov -8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
mov (rp,i,8), X0
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov %rdx, w0
|
|
Packit |
5c3484 |
add %rax, X0
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w2, X0
|
|
Packit |
5c3484 |
mov X1, -8(rp,i,8)
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
L(lo3): add %rax, X0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
add w3, X0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
mov 8(rp,i,8), X1
|
|
Packit |
5c3484 |
add %rax, X1
|
|
Packit |
5c3484 |
mov %rdx, w2
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mov X0, (rp,i,8)
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w0, X1
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
add %rax, X1
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mov 16(rp,i,8), X0
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add $4, i
|
|
Packit |
5c3484 |
jnc L(top)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(end): imul v1, %rax
|
|
Packit |
5c3484 |
add %rax, X0
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add w2, X0
|
|
Packit |
5c3484 |
mov I(-8(up),-16(up,i,8)), %rax
|
|
Packit |
5c3484 |
imul v0, %rax
|
|
Packit |
5c3484 |
add X0, %rax
|
|
Packit |
5c3484 |
mov X1, I(-16(rp),-24(rp,i,8))
|
|
Packit |
5c3484 |
add w3, %rax
|
|
Packit |
5c3484 |
mov %rax, I(-8(rp),-16(rp,i,8))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
add $2, n
|
|
Packit |
5c3484 |
lea 16(vp), vp
|
|
Packit |
5c3484 |
lea -16(up), up
|
|
Packit |
5c3484 |
cmp $-2, n
|
|
Packit |
5c3484 |
jl L(outer)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
pop %r15
|
|
Packit |
5c3484 |
pop %r14
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
jnz L(cor0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(cor1):mov (vp), v0
|
|
Packit |
5c3484 |
mov 8(vp), v1
|
|
Packit |
5c3484 |
mov -16(up), %rax
|
|
Packit |
5c3484 |
mul v0 C u0 x v2
|
|
Packit |
5c3484 |
add -16(rp), %rax C FIXME: rp[0] still available in reg?
|
|
Packit |
5c3484 |
adc -8(rp), %rdx C FIXME: rp[1] still available in reg?
|
|
Packit |
5c3484 |
mov -8(up), %r10
|
|
Packit |
5c3484 |
imul v0, %r10
|
|
Packit |
5c3484 |
mov -16(up), %r11
|
|
Packit |
5c3484 |
imul v1, %r11
|
|
Packit |
5c3484 |
mov %rax, -16(rp)
|
|
Packit |
5c3484 |
add %r10, %r11
|
|
Packit |
5c3484 |
add %rdx, %r11
|
|
Packit |
5c3484 |
mov %r11, -8(rp)
|
|
Packit |
5c3484 |
pop %r13
|
|
Packit |
5c3484 |
pop %r12
|
|
Packit |
5c3484 |
pop %rbp
|
|
Packit |
5c3484 |
pop %rbx
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(cor0):mov (vp), %r11
|
|
Packit |
5c3484 |
imul -8(up), %r11
|
|
Packit |
5c3484 |
add %rax, %r11
|
|
Packit |
5c3484 |
mov %r11, -8(rp)
|
|
Packit |
5c3484 |
pop %r13
|
|
Packit |
5c3484 |
pop %r12
|
|
Packit |
5c3484 |
pop %rbp
|
|
Packit |
5c3484 |
pop %rbx
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(small):
|
|
Packit |
5c3484 |
cmp $2, n
|
|
Packit |
5c3484 |
jae L(gt1)
|
|
Packit |
5c3484 |
L(n1): imul (vp_param), %rax
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
L(gt1): ja L(gt2)
|
|
Packit |
5c3484 |
L(n2): mov (vp_param), %r9
|
|
Packit |
5c3484 |
mul %r9
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
mov 8(up), %rax
|
|
Packit |
5c3484 |
imul %r9, %rax
|
|
Packit |
5c3484 |
add %rax, %rdx
|
|
Packit |
5c3484 |
mov 8(vp), %r9
|
|
Packit |
5c3484 |
mov (up), %rcx
|
|
Packit |
5c3484 |
imul %r9, %rcx
|
|
Packit |
5c3484 |
add %rcx, %rdx
|
|
Packit |
5c3484 |
mov %rdx, 8(rp)
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
L(gt2):
|
|
Packit |
5c3484 |
L(n3): mov (vp_param), %r9
|
|
Packit |
5c3484 |
mul %r9 C u0 x v0
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
mov %rdx, %r10
|
|
Packit |
5c3484 |
mov 8(up), %rax
|
|
Packit |
5c3484 |
mul %r9 C u1 x v0
|
|
Packit |
5c3484 |
imul 16(up), %r9 C u2 x v0
|
|
Packit |
5c3484 |
add %rax, %r10
|
|
Packit |
5c3484 |
adc %rdx, %r9
|
|
Packit |
5c3484 |
mov 8(vp), %r11
|
|
Packit |
5c3484 |
mov (up), %rax
|
|
Packit |
5c3484 |
mul %r11 C u0 x v1
|
|
Packit |
5c3484 |
add %rax, %r10
|
|
Packit |
5c3484 |
adc %rdx, %r9
|
|
Packit |
5c3484 |
imul 8(up), %r11 C u1 x v1
|
|
Packit |
5c3484 |
add %r11, %r9
|
|
Packit |
5c3484 |
mov %r10, 8(rp)
|
|
Packit |
5c3484 |
mov 16(vp), %r10
|
|
Packit |
5c3484 |
mov (up), %rax
|
|
Packit |
5c3484 |
imul %rax, %r10 C u0 x v2
|
|
Packit |
5c3484 |
add %r10, %r9
|
|
Packit |
5c3484 |
mov %r9, 16(rp)
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
EPILOGUE()
|
|
rpm-build |
c3cd4f |
CF_PROT
|