|
Packit |
5c3484 |
dnl AMD64 mpn_sqr_basecase optimised for Intel Haswell.
|
|
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 sqr_diag_addlsh1
|
|
Packit |
5c3484 |
C AMD K8,K9 n/a n/a n/a
|
|
Packit |
5c3484 |
C AMD K10 n/a n/a n/a
|
|
Packit |
5c3484 |
C AMD bull n/a n/a n/a
|
|
Packit |
5c3484 |
C AMD pile n/a n/a n/a
|
|
Packit |
5c3484 |
C AMD steam ? ? ?
|
|
Packit |
5c3484 |
C AMD bobcat n/a n/a n/a
|
|
Packit |
5c3484 |
C AMD jaguar ? ? ?
|
|
Packit |
5c3484 |
C Intel P4 n/a n/a n/a
|
|
Packit |
5c3484 |
C Intel core n/a n/a n/a
|
|
Packit |
5c3484 |
C Intel NHM n/a n/a n/a
|
|
Packit |
5c3484 |
C Intel SBR n/a n/a n/a
|
|
Packit |
5c3484 |
C Intel IBR n/a n/a n/a
|
|
Packit |
5c3484 |
C Intel HWL 1.86 2.15 ~2.5
|
|
Packit |
5c3484 |
C Intel BWL ? ? ?
|
|
Packit |
5c3484 |
C Intel atom n/a n/a n/a
|
|
Packit |
5c3484 |
C VIA nano n/a n/a n/a
|
|
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, except
|
|
Packit |
5c3484 |
C that the sqr_diag_addlsh1 loop was manually written.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C TODO
|
|
Packit |
5c3484 |
C * Replace current unoptimised sqr_diag_addlsh1 loop; 1.75 c/l might be
|
|
Packit |
5c3484 |
C possible.
|
|
Packit |
5c3484 |
C * Consider splitting outer loop into 2, one for n = 1 (mod 2) and one for
|
|
Packit |
5c3484 |
C n = 0 (mod 2). These loops could fall into specific "corner" code.
|
|
Packit |
5c3484 |
C * Consider splitting outer loop into 4.
|
|
Packit |
5c3484 |
C * Streamline pointer updates.
|
|
Packit |
5c3484 |
C * Perhaps suppress a few more xor insns in feed-in code.
|
|
Packit |
5c3484 |
C * Make sure we write no dead registers in feed-in code.
|
|
Packit |
5c3484 |
C * We might use 32-bit size ops, since n >= 2^32 is non-terminating. Watch
|
|
Packit |
5c3484 |
C out for negative sizes being zero-extended, though.
|
|
Packit |
5c3484 |
C * Provide straight-line code for n = 4; then look for simplifications in
|
|
Packit |
5c3484 |
C main code.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`rp', `%rdi')
|
|
Packit |
5c3484 |
define(`up', `%rsi')
|
|
Packit |
5c3484 |
define(`un_param',`%rdx')
|
|
Packit |
5c3484 |
|
|
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_sqr_basecase)
|
|
Packit |
5c3484 |
FUNC_ENTRY(3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
cmp $2, un_param
|
|
Packit |
5c3484 |
jae L(gt1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov (up), %rdx
|
|
Packit |
5c3484 |
mulx( %rdx, %rax, %rdx)
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
mov %rdx, 8(rp)
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(gt1): jne L(gt2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov (up), %rdx
|
|
Packit |
5c3484 |
mov 8(up), %rcx
|
|
Packit |
5c3484 |
mulx( %rcx, %r9, %r10) C v0 * v1 W 1 2
|
|
Packit |
5c3484 |
mulx( %rdx, %rax, %r8) C v0 * v0 W 0 1
|
|
Packit |
5c3484 |
mov %rcx, %rdx
|
|
Packit |
5c3484 |
mulx( %rdx, %r11, %rdx) C v1 * v1 W 2 3
|
|
Packit |
5c3484 |
add %r9, %r9 C W 1
|
|
Packit |
5c3484 |
adc %r10, %r10 C W 2
|
|
Packit |
5c3484 |
adc $0, %rdx C W 3
|
|
Packit |
5c3484 |
add %r9, %r8 C W 1
|
|
Packit |
5c3484 |
adc %r11, %r10 C W 2
|
|
Packit |
5c3484 |
adc $0, %rdx C W 3
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
mov %r8, 8(rp)
|
|
Packit |
5c3484 |
mov %r10, 16(rp)
|
|
Packit |
5c3484 |
mov %rdx, 24(rp)
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(gt2): cmp $4, un_param
|
|
Packit |
5c3484 |
jae L(gt3)
|
|
Packit |
5c3484 |
define(`v0', `%r8')
|
|
Packit |
5c3484 |
define(`v1', `%r9')
|
|
Packit |
5c3484 |
define(`w0', `%r10')
|
|
Packit |
5c3484 |
define(`w2', `%r11')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov (up), v0
|
|
Packit |
5c3484 |
mov 8(up), %rdx
|
|
Packit |
5c3484 |
mov %rdx, v1
|
|
Packit |
5c3484 |
mulx( v0, w2, %rax)
|
|
Packit |
5c3484 |
mov 16(up), %rdx
|
|
Packit |
5c3484 |
mulx( v0, w0, %rcx)
|
|
Packit |
5c3484 |
mov w2, %r8
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
adc $0, %rcx
|
|
Packit |
5c3484 |
mulx( v1, %rdx, %rax)
|
|
Packit |
5c3484 |
add %rcx, %rdx
|
|
Packit |
5c3484 |
mov %rdx, 24(rp)
|
|
Packit |
5c3484 |
adc $0, %rax
|
|
Packit |
5c3484 |
mov %rax, 32(rp)
|
|
Packit |
5c3484 |
xor R32(%rcx), R32(%rcx)
|
|
Packit |
5c3484 |
mov (up), %rdx
|
|
Packit |
5c3484 |
mulx( %rdx, %rax, w2)
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
add %r8, %r8
|
|
Packit |
5c3484 |
adc w0, w0
|
|
Packit |
5c3484 |
setc R8(%rcx)
|
|
Packit |
5c3484 |
mov 8(up), %rdx
|
|
Packit |
5c3484 |
mulx( %rdx, %rax, %rdx)
|
|
Packit |
5c3484 |
add w2, %r8
|
|
Packit |
5c3484 |
adc %rax, w0
|
|
Packit |
5c3484 |
mov %r8, 8(rp)
|
|
Packit |
5c3484 |
mov w0, 16(rp)
|
|
Packit |
5c3484 |
mov 24(rp), %r8
|
|
Packit |
5c3484 |
mov 32(rp), w0
|
|
Packit |
5c3484 |
lea (%rdx,%rcx), w2
|
|
Packit |
5c3484 |
adc %r8, %r8
|
|
Packit |
5c3484 |
adc w0, w0
|
|
Packit |
5c3484 |
setc R8(%rcx)
|
|
Packit |
5c3484 |
mov 16(up), %rdx
|
|
Packit |
5c3484 |
mulx( %rdx, %rax, %rdx)
|
|
Packit |
5c3484 |
add w2, %r8
|
|
Packit |
5c3484 |
adc %rax, w0
|
|
Packit |
5c3484 |
mov %r8, 24(rp)
|
|
Packit |
5c3484 |
mov w0, 32(rp)
|
|
Packit |
5c3484 |
adc %rcx, %rdx
|
|
Packit |
5c3484 |
mov %rdx, 40(rp)
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(gt3):
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`v0', `%r8')
|
|
Packit |
5c3484 |
define(`v1', `%r9')
|
|
Packit |
5c3484 |
define(`w0', `%r10')
|
|
Packit |
5c3484 |
define(`w1', `%r11')
|
|
Packit |
5c3484 |
define(`w2', `%rbx')
|
|
Packit |
5c3484 |
define(`w3', `%rbp')
|
|
Packit |
5c3484 |
define(`un', `%r12')
|
|
Packit |
5c3484 |
define(`n', `%rcx')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`X0', `%r13')
|
|
Packit |
5c3484 |
define(`X1', `%r14')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(do_mul_2):
|
|
Packit |
5c3484 |
push %rbx
|
|
Packit |
5c3484 |
push %rbp
|
|
Packit |
5c3484 |
push %r12
|
|
Packit |
5c3484 |
push %r13
|
|
Packit |
5c3484 |
push %r14
|
|
Packit |
5c3484 |
mov $0, R32(un)
|
|
Packit |
5c3484 |
sub un_param, un C free up rdx
|
|
Packit |
5c3484 |
push un
|
|
Packit |
5c3484 |
mov (up), v0
|
|
Packit |
5c3484 |
mov 8(up), %rdx
|
|
Packit |
5c3484 |
lea 2(un), n
|
|
Packit |
5c3484 |
sar $2, n C FIXME: suppress, change loop?
|
|
Packit |
5c3484 |
inc un C decrement |un|
|
|
Packit |
5c3484 |
mov %rdx, v1
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
test $1, R8(un)
|
|
Packit |
5c3484 |
jnz L(mx1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(mx0): mulx( v0, w2, w1)
|
|
Packit |
5c3484 |
mov 16(up), %rdx
|
|
Packit |
5c3484 |
mov w2, 8(rp)
|
|
Packit |
5c3484 |
xor w2, w2
|
|
Packit |
5c3484 |
mulx( v0, w0, w3)
|
|
Packit |
5c3484 |
test $2, R8(un)
|
|
Packit |
5c3484 |
jz L(m00)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(m10): lea -8(rp), rp
|
|
Packit |
5c3484 |
lea -8(up), up
|
|
Packit |
5c3484 |
jmp L(mlo2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(m00): lea 8(up), up
|
|
Packit |
5c3484 |
lea 8(rp), rp
|
|
Packit |
5c3484 |
jmp L(mlo0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(mx1): mulx( v0, w0, w3)
|
|
Packit |
5c3484 |
mov 16(up), %rdx
|
|
Packit |
5c3484 |
mov w0, 8(rp)
|
|
Packit |
5c3484 |
xor w0, w0
|
|
Packit |
5c3484 |
mulx( v0, w2, w1)
|
|
Packit |
5c3484 |
test $2, R8(un)
|
|
Packit |
5c3484 |
jz L(mlo3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(m01): lea 16(rp), rp
|
|
Packit |
5c3484 |
lea 16(up), up
|
|
Packit |
5c3484 |
jmp L(mlo1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(32)
|
|
Packit |
5c3484 |
L(mtop):mulx( v1, %rax, w0)
|
|
Packit |
5c3484 |
add %rax, w2 C 0
|
|
Packit |
5c3484 |
mov (up), %rdx
|
|
Packit |
5c3484 |
mulx( v0, %rax, w1)
|
|
Packit |
5c3484 |
adc $0, w0 C 1
|
|
Packit |
5c3484 |
add %rax, w2 C 0
|
|
Packit |
5c3484 |
L(mlo1):adc $0, w1 C 1
|
|
Packit |
5c3484 |
add w3, w2 C 0
|
|
Packit |
5c3484 |
mov w2, (rp) C 0
|
|
Packit |
5c3484 |
adc $0, w1 C 1
|
|
Packit |
5c3484 |
mulx( v1, %rax, w2)
|
|
Packit |
5c3484 |
add %rax, w0 C 1
|
|
Packit |
5c3484 |
mov 8(up), %rdx
|
|
Packit |
5c3484 |
adc $0, w2 C 2
|
|
Packit |
5c3484 |
mulx( v0, %rax, w3)
|
|
Packit |
5c3484 |
add %rax, w0 C 1
|
|
Packit |
5c3484 |
adc $0, w3 C 2
|
|
Packit |
5c3484 |
L(mlo0):add w1, w0 C 1
|
|
Packit |
5c3484 |
mov w0, 8(rp) C 1
|
|
Packit |
5c3484 |
adc $0, w3 C 2
|
|
Packit |
5c3484 |
mulx( v1, %rax, w0)
|
|
Packit |
5c3484 |
add %rax, w2 C 2
|
|
Packit |
5c3484 |
mov 16(up), %rdx
|
|
Packit |
5c3484 |
mulx( v0, %rax, w1)
|
|
Packit |
5c3484 |
adc $0, w0 C 3
|
|
Packit |
5c3484 |
add %rax, w2 C 2
|
|
Packit |
5c3484 |
adc $0, w1 C 3
|
|
Packit |
5c3484 |
L(mlo3):add w3, w2 C 2
|
|
Packit |
5c3484 |
mov w2, 16(rp) C 2
|
|
Packit |
5c3484 |
adc $0, w1 C 3
|
|
Packit |
5c3484 |
mulx( v1, %rax, w2)
|
|
Packit |
5c3484 |
add %rax, w0 C 3
|
|
Packit |
5c3484 |
mov 24(up), %rdx
|
|
Packit |
5c3484 |
adc $0, w2 C 4
|
|
Packit |
5c3484 |
mulx( v0, %rax, w3)
|
|
Packit |
5c3484 |
add %rax, w0 C 3
|
|
Packit |
5c3484 |
adc $0, w3 C 4
|
|
Packit |
5c3484 |
L(mlo2):add w1, w0 C 3
|
|
Packit |
5c3484 |
lea 32(up), up
|
|
Packit |
5c3484 |
mov w0, 24(rp) C 3
|
|
Packit |
5c3484 |
adc $0, w3 C 4
|
|
Packit |
5c3484 |
inc n
|
|
Packit |
5c3484 |
lea 32(rp), rp
|
|
Packit |
5c3484 |
jnz L(mtop)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(mend):mulx( v1, %rdx, %rax)
|
|
Packit |
5c3484 |
add %rdx, w2
|
|
Packit |
5c3484 |
adc $0, %rax
|
|
Packit |
5c3484 |
add w3, w2
|
|
Packit |
5c3484 |
mov w2, (rp)
|
|
Packit |
5c3484 |
adc $0, %rax
|
|
Packit |
5c3484 |
mov %rax, 8(rp)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 16(up), up
|
|
Packit |
5c3484 |
lea -16(rp), rp
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(do_addmul_2):
|
|
Packit |
5c3484 |
L(outer):
|
|
Packit |
5c3484 |
lea (up,un,8), up C put back up to 2 positions above last time
|
|
Packit |
5c3484 |
lea 48(rp,un,8), rp C put back rp to 4 positions above last time
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov -8(up), v0 C shared between addmul_2 and corner
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
add $2, un C decrease |un|
|
|
Packit |
5c3484 |
cmp $-2, un
|
|
Packit |
5c3484 |
jge L(corner)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov (up), v1
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 1(un), n
|
|
Packit |
5c3484 |
sar $2, n C FIXME: suppress, change loop?
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov v1, %rdx
|
|
Packit |
5c3484 |
test $1, R8(un)
|
|
Packit |
5c3484 |
jnz L(bx1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(bx0): mov (rp), X0
|
|
Packit |
5c3484 |
mov 8(rp), X1
|
|
Packit |
5c3484 |
mulx( v0, %rax, w1)
|
|
Packit |
5c3484 |
add %rax, X0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov X0, (rp)
|
|
Packit |
5c3484 |
xor w2, w2
|
|
Packit |
5c3484 |
test $2, R8(un)
|
|
Packit |
5c3484 |
jnz L(b10)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b00): mov 8(up), %rdx
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
lea 16(up), up
|
|
Packit |
5c3484 |
jmp L(lo0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b10): mov 8(up), %rdx
|
|
Packit |
5c3484 |
mov 16(rp), X0
|
|
Packit |
5c3484 |
lea 32(up), up
|
|
Packit |
5c3484 |
inc n
|
|
Packit |
5c3484 |
mulx( v0, %rax, w3)
|
|
Packit |
5c3484 |
jz L(ex)
|
|
Packit |
5c3484 |
jmp L(lo2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(bx1): mov (rp), X1
|
|
Packit |
5c3484 |
mov 8(rp), X0
|
|
Packit |
5c3484 |
mulx( v0, %rax, w3)
|
|
Packit |
5c3484 |
mov 8(up), %rdx
|
|
Packit |
5c3484 |
add %rax, X1
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
xor w0, w0
|
|
Packit |
5c3484 |
mov X1, (rp)
|
|
Packit |
5c3484 |
mulx( v0, %rax, w1)
|
|
Packit |
5c3484 |
test $2, R8(un)
|
|
Packit |
5c3484 |
jz L(b11)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b01): mov 16(rp), X1
|
|
Packit |
5c3484 |
lea 24(rp), rp
|
|
Packit |
5c3484 |
lea 24(up), up
|
|
Packit |
5c3484 |
jmp L(lo1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b11): lea 8(rp), rp
|
|
Packit |
5c3484 |
lea 8(up), up
|
|
Packit |
5c3484 |
jmp L(lo3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(32)
|
|
Packit |
5c3484 |
L(top): mulx( v0, %rax, w3)
|
|
Packit |
5c3484 |
add w0, X1
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
L(lo2): add %rax, X1
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mulx( v1, %rax, w0)
|
|
Packit |
5c3484 |
add %rax, X0
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
lea 32(rp), rp
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
mov -16(up), %rdx
|
|
Packit |
5c3484 |
mov X1, -24(rp)
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add w2, X0
|
|
Packit |
5c3484 |
mov -8(rp), X1
|
|
Packit |
5c3484 |
mulx( v0, %rax, w1)
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
L(lo1): add %rax, X0
|
|
Packit |
5c3484 |
mulx( v1, %rax, w2)
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
add w3, X0
|
|
Packit |
5c3484 |
mov X0, -16(rp)
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
add %rax, X1
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
add w0, X1
|
|
Packit |
5c3484 |
mov -8(up), %rdx
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
L(lo0): mulx( v0, %rax, w3)
|
|
Packit |
5c3484 |
add %rax, X1
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov (rp), X0
|
|
Packit |
5c3484 |
mulx( v1, %rax, w0)
|
|
Packit |
5c3484 |
add %rax, X0
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
mov X1, -8(rp)
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov (up), %rdx
|
|
Packit |
5c3484 |
add w2, X0
|
|
Packit |
5c3484 |
mulx( v0, %rax, w1)
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
L(lo3): add %rax, X0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mulx( v1, %rax, w2)
|
|
Packit |
5c3484 |
add w3, X0
|
|
Packit |
5c3484 |
mov 8(rp), X1
|
|
Packit |
5c3484 |
mov X0, (rp)
|
|
Packit |
5c3484 |
mov 16(rp), X0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
add %rax, X1
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
mov 8(up), %rdx
|
|
Packit |
5c3484 |
lea 32(up), up
|
|
Packit |
5c3484 |
inc n
|
|
Packit |
5c3484 |
jnz L(top)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(end): mulx( v0, %rax, w3)
|
|
Packit |
5c3484 |
add w0, X1
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
L(ex): add %rax, X1
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mulx( v1, %rdx, %rax)
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
mov X1, 8(rp)
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add w2, %rdx
|
|
Packit |
5c3484 |
adc $0, %rax
|
|
Packit |
5c3484 |
add %rdx, w3
|
|
Packit |
5c3484 |
mov w3, 16(rp)
|
|
Packit |
5c3484 |
adc $0, %rax
|
|
Packit |
5c3484 |
mov %rax, 24(rp)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
jmp L(outer) C loop until a small corner remains
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(corner):
|
|
Packit |
5c3484 |
pop un
|
|
Packit |
5c3484 |
mov (up), %rdx
|
|
Packit |
5c3484 |
jg L(small_corner)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov %rdx, v1
|
|
Packit |
5c3484 |
mov (rp), X0
|
|
Packit |
5c3484 |
mov %rax, X1 C Tricky rax reuse of last iteration
|
|
Packit |
5c3484 |
mulx( v0, %rax, w1)
|
|
Packit |
5c3484 |
add %rax, X0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov X0, (rp)
|
|
Packit |
5c3484 |
mov 8(up), %rdx
|
|
Packit |
5c3484 |
mulx( v0, %rax, w3)
|
|
Packit |
5c3484 |
add %rax, X1
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mulx( v1, %rdx, %rax)
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
mov X1, 8(rp)
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add w3, %rdx
|
|
Packit |
5c3484 |
mov %rdx, 16(rp)
|
|
Packit |
5c3484 |
adc $0, %rax
|
|
Packit |
5c3484 |
mov %rax, 24(rp)
|
|
Packit |
5c3484 |
lea 32(rp), rp
|
|
Packit |
5c3484 |
lea 16(up), up
|
|
Packit |
5c3484 |
jmp L(com)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(small_corner):
|
|
Packit |
5c3484 |
mulx( v0, X1, w3)
|
|
Packit |
5c3484 |
add %rax, X1 C Tricky rax reuse of last iteration
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov X1, (rp)
|
|
Packit |
5c3484 |
mov w3, 8(rp)
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
lea 8(up), up
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(com):
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(sqr_diag_addlsh1):
|
|
Packit |
5c3484 |
lea 8(up,un,8), up C put back up at its very beginning
|
|
Packit |
5c3484 |
lea (rp,un,8), rp
|
|
Packit |
5c3484 |
lea (rp,un,8), rp C put back rp at its very beginning
|
|
Packit |
5c3484 |
inc un
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov -8(up), %rdx
|
|
Packit |
5c3484 |
xor R32(%rbx), R32(%rbx) C clear CF as side effect
|
|
Packit |
5c3484 |
mulx( %rdx, %rax, %r10)
|
|
Packit |
5c3484 |
mov %rax, 8(rp)
|
|
Packit |
5c3484 |
mov 16(rp), %r8
|
|
Packit |
5c3484 |
mov 24(rp), %r9
|
|
Packit |
5c3484 |
jmp L(dm)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(dtop):mov 32(rp), %r8
|
|
Packit |
5c3484 |
mov 40(rp), %r9
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
lea (%rdx,%rbx), %r10
|
|
Packit |
5c3484 |
L(dm): adc %r8, %r8
|
|
Packit |
5c3484 |
adc %r9, %r9
|
|
Packit |
5c3484 |
setc R8(%rbx)
|
|
Packit |
5c3484 |
mov (up), %rdx
|
|
Packit |
5c3484 |
lea 8(up), up
|
|
Packit |
5c3484 |
mulx( %rdx, %rax, %rdx)
|
|
Packit |
5c3484 |
add %r10, %r8
|
|
Packit |
5c3484 |
adc %rax, %r9
|
|
Packit |
5c3484 |
mov %r8, 16(rp)
|
|
Packit |
5c3484 |
mov %r9, 24(rp)
|
|
Packit |
5c3484 |
inc un
|
|
Packit |
5c3484 |
jnz L(dtop)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(dend):adc %rbx, %rdx
|
|
Packit |
5c3484 |
mov %rdx, 32(rp)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
pop %r14
|
|
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 |
EPILOGUE()
|