|
Packit |
5c3484 |
dnl X86-64 mpn_sqr_basecase optimised for Intel Nehalem/Westmere.
|
|
Packit |
5c3484 |
dnl It also seems good for Conroe/Wolfdale.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Contributed to the GNU project by Torbjörn Granlund.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Copyright 2008, 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
|
|
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 4.9 4.18-4.25 3.87
|
|
Packit |
5c3484 |
C Intel NHM 3.8 4.06-4.2 3.5
|
|
Packit |
5c3484 |
C Intel SBR
|
|
Packit |
5c3484 |
C Intel IBR
|
|
Packit |
5c3484 |
C Intel HWL
|
|
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 Code structure:
|
|
Packit |
5c3484 |
C
|
|
Packit |
5c3484 |
C
|
|
Packit |
5c3484 |
C m_2(0m4) m_2(2m4) m_2(1m4) m_2(3m4)
|
|
Packit |
5c3484 |
C | | | |
|
|
Packit |
5c3484 |
C | | | |
|
|
Packit |
5c3484 |
C | | | |
|
|
Packit |
5c3484 |
C \|/ \|/ \|/ \|/
|
|
Packit |
5c3484 |
C ____________ ____________
|
|
Packit |
5c3484 |
C / \ / \
|
|
Packit |
5c3484 |
C \|/ \ \|/ \
|
|
Packit |
5c3484 |
C am_2(3m4) am_2(1m4) am_2(0m4) am_2(2m4)
|
|
Packit |
5c3484 |
C \ /|\ \ /|\
|
|
Packit |
5c3484 |
C \____________/ \____________/
|
|
Packit |
5c3484 |
C \ /
|
|
Packit |
5c3484 |
C \ /
|
|
Packit |
5c3484 |
C \ /
|
|
Packit |
5c3484 |
C tail(0m2) tail(1m2)
|
|
Packit |
5c3484 |
C \ /
|
|
Packit |
5c3484 |
C \ /
|
|
Packit |
5c3484 |
C sqr_diag_addlsh1
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C TODO
|
|
Packit |
5c3484 |
C * Tune. None done so far.
|
|
Packit |
5c3484 |
C * Currently 2761 bytes, making it smaller would be nice.
|
|
Packit |
5c3484 |
C * Consider using a jumptab-based entry sequence. One might even use a mask-
|
|
Packit |
5c3484 |
C less sequence, if the table is large enough to support tuneup's needs.
|
|
Packit |
5c3484 |
C The code would be, using non-PIC code,
|
|
Packit |
5c3484 |
C lea tab(%rip),%rax; jmp *(n,%rax)
|
|
Packit |
5c3484 |
C or,
|
|
Packit |
5c3484 |
C lea tab(%rip),%rax; lea (%rip),%rbx; add (n,%rax),%rbx; jmp *%rbx
|
|
Packit |
5c3484 |
C using PIC code. The table entries would be Ln1,Ln2,Ln3,Lm0,Lm1,Lm2,Lm3,..
|
|
Packit |
5c3484 |
C with the last four entries repeated a safe number of times.
|
|
Packit |
5c3484 |
C * Consider expanding feed-in code in order to avoid zeroing registers.
|
|
Packit |
5c3484 |
C * Zero consistently with xor.
|
|
Packit |
5c3484 |
C * Check if using "lea (reg),reg" should be done in more places; we have some
|
|
Packit |
5c3484 |
C explicit "mov %rax,reg" now.
|
|
Packit |
5c3484 |
C * Try zeroing with xor in m2 loops.
|
|
Packit |
5c3484 |
C * Try re-rolling the m2 loops to avoid the current 9 insn code duplication
|
|
Packit |
5c3484 |
C between loop header and wind-down code.
|
|
Packit |
5c3484 |
C * Consider adc reg,reg instead of adc $0,reg in m2 loops. This save a byte.
|
|
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 |
C Define this to $1 to use late loop index variable as zero, $2 to use an
|
|
Packit |
5c3484 |
C explicit $0.
|
|
Packit |
5c3484 |
define(`Z',`$1')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`rp', `%rdi')
|
|
Packit |
5c3484 |
define(`up', `%rsi')
|
|
Packit |
5c3484 |
define(`n_param', `%rdx')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`n', `%r8')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`v0', `%r10')
|
|
Packit |
5c3484 |
define(`v1', `%r11')
|
|
Packit |
5c3484 |
define(`w0', `%rbx')
|
|
Packit |
5c3484 |
define(`w1', `%rcx')
|
|
Packit |
5c3484 |
define(`w2', `%rbp')
|
|
Packit |
5c3484 |
define(`w3', `%r9')
|
|
Packit |
5c3484 |
define(`i', `%r13')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`X0', `%r12')
|
|
Packit |
5c3484 |
define(`X1', `%r14')
|
|
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 |
define(`ALIGNx', `ALIGN(16)')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
define(`N', 85)
|
|
Packit |
5c3484 |
ifdef(`N',,`define(`N',0)')
|
|
Packit |
5c3484 |
define(`MOV', `ifelse(eval(N & $3),0,`mov $1, $2',`lea ($1), $2')')
|
|
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 $4, n_param
|
|
Packit |
5c3484 |
jl L(small)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
push %rbx
|
|
Packit |
5c3484 |
push %rbp
|
|
Packit |
5c3484 |
push %r12
|
|
Packit |
5c3484 |
push %r13
|
|
Packit |
5c3484 |
push %r14
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov (up), v0
|
|
Packit |
5c3484 |
mov 8(up), %rax
|
|
Packit |
5c3484 |
mov %rax, v1
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov $1, R32(n)
|
|
Packit |
5c3484 |
sub n_param, n C n = -n_param+1
|
|
Packit |
5c3484 |
push n
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea (up,n_param,8), up
|
|
Packit |
5c3484 |
lea (rp,n_param,8), rp
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
test $1, R8(n)
|
|
Packit |
5c3484 |
jnz L(bx1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(bx0): test $2, R8(n)
|
|
Packit |
5c3484 |
mov %rax, (rp,n,8)
|
|
Packit |
5c3484 |
jnz L(b10)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b00): lea (n), i C n = 5, 9, ...
|
|
Packit |
5c3484 |
mov %rdx, w1 C FIXME: Use lea?
|
|
Packit |
5c3484 |
xor R32(w2), R32(w2)
|
|
Packit |
5c3484 |
jmp L(m2e0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b10): lea 2(n), i C n = 7, 11, ...
|
|
Packit |
5c3484 |
mov 8(up,n,8), %rax
|
|
Packit |
5c3484 |
mov %rdx, w3 C FIXME: Use lea?
|
|
Packit |
5c3484 |
xor R32(w0), R32(w0)
|
|
Packit |
5c3484 |
xor R32(w1), R32(w1)
|
|
Packit |
5c3484 |
jmp L(m2e2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(bx1): test $2, R8(n)
|
|
Packit |
5c3484 |
mov %rax, (rp,n,8)
|
|
Packit |
5c3484 |
jz L(b11)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b01): lea 1(n), i C n = 6, 10, ...
|
|
Packit |
5c3484 |
mov %rdx, w0 C FIXME: Use lea?
|
|
Packit |
5c3484 |
xor R32(w1), R32(w1)
|
|
Packit |
5c3484 |
jmp L(m2e1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b11): lea -1(n), i C n = 4, 8, 12, ...
|
|
Packit |
5c3484 |
mov %rdx, w2 C FIXME: Use lea?
|
|
Packit |
5c3484 |
xor R32(w3), R32(w3)
|
|
Packit |
5c3484 |
jmp L(m2e3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGNx
|
|
Packit |
5c3484 |
L(m2top1):
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov -8(up,i,8), %rax
|
|
Packit |
5c3484 |
mov w3, -8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
adc $0, R32(w1)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
adc %rdx, w1
|
|
Packit |
5c3484 |
L(m2e1):mov $0, R32(w2)
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
mov w0, (rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w1
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
adc $0, R32(w2)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w1
|
|
Packit |
5c3484 |
adc %rdx, w2
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov $0, R32(w3)
|
|
Packit |
5c3484 |
add %rax, w1
|
|
Packit |
5c3484 |
adc %rdx, w2
|
|
Packit |
5c3484 |
adc $0, R32(w3)
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
mov w1, 8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w3
|
|
Packit |
5c3484 |
mov $0, R32(w0)
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
adc %rdx, w3
|
|
Packit |
5c3484 |
adc $0, R32(w0)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
mov $0, R32(w1)
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
mov w2, 16(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
add $4, i
|
|
Packit |
5c3484 |
js L(m2top1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov I(-8(up),-8(up,i,8)), %rax
|
|
Packit |
5c3484 |
mov w3, I(-8(rp),-8(rp,i,8))
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
adc R32(w1), R32(w1)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, %rax
|
|
Packit |
5c3484 |
adc w1, %rdx
|
|
Packit |
5c3484 |
mov %rax, I((rp),(rp,i,8))
|
|
Packit |
5c3484 |
mov %rdx, I(8(rp),8(rp,i,8))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
add $2, n C decrease |n|
|
|
Packit |
5c3484 |
jmp L(am2o3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGNx
|
|
Packit |
5c3484 |
L(m2top3):
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov -8(up,i,8), %rax
|
|
Packit |
5c3484 |
mov w3, -8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
adc $0, R32(w1)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
adc %rdx, w1
|
|
Packit |
5c3484 |
mov $0, R32(w2)
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
mov w0, (rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w1
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
adc $0, R32(w2)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w1
|
|
Packit |
5c3484 |
adc %rdx, w2
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov $0, R32(w3)
|
|
Packit |
5c3484 |
add %rax, w1
|
|
Packit |
5c3484 |
adc %rdx, w2
|
|
Packit |
5c3484 |
adc $0, R32(w3)
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
mov w1, 8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w3
|
|
Packit |
5c3484 |
L(m2e3):mov $0, R32(w0)
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
adc %rdx, w3
|
|
Packit |
5c3484 |
adc $0, R32(w0)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
mov $0, R32(w1)
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
mov w2, 16(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
add $4, i
|
|
Packit |
5c3484 |
js L(m2top3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov I(-8(up),-8(up,i,8)), %rax
|
|
Packit |
5c3484 |
mov w3, I(-8(rp),-8(rp,i,8))
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
adc R32(w1), R32(w1)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, %rax
|
|
Packit |
5c3484 |
adc w1, %rdx
|
|
Packit |
5c3484 |
mov %rax, I((rp),(rp,i,8))
|
|
Packit |
5c3484 |
mov %rdx, I(8(rp),8(rp,i,8))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
add $2, n C decrease |n|
|
|
Packit |
5c3484 |
cmp $-1, n
|
|
Packit |
5c3484 |
jz L(cor1) C jumps iff entry n = 4
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(am2o1):
|
|
Packit |
5c3484 |
mov -8(up,n,8), v0
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mov %rax, v1
|
|
Packit |
5c3484 |
lea 1(n), i
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, X1
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 128)
|
|
Packit |
5c3484 |
mov (rp,n,8), w1
|
|
Packit |
5c3484 |
xor R32(w2), R32(w2)
|
|
Packit |
5c3484 |
mov 8(up,n,8), %rax
|
|
Packit |
5c3484 |
xor R32(w3), R32(w3)
|
|
Packit |
5c3484 |
jmp L(lo1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGNx
|
|
Packit |
5c3484 |
L(am2top1):
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, w1
|
|
Packit |
5c3484 |
adc %rax, w2
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
MOV( %rdx, w3, 1)
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
L(lo1): mul v0
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
mov X1, -8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 2)
|
|
Packit |
5c3484 |
adc $0, X1
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
MOV( %rdx, w0, 4)
|
|
Packit |
5c3484 |
mov (rp,i,8), w1
|
|
Packit |
5c3484 |
add w1, w2
|
|
Packit |
5c3484 |
adc %rax, w3
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w2, X0
|
|
Packit |
5c3484 |
adc %rax, X1
|
|
Packit |
5c3484 |
mov X0, (rp,i,8)
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 8)
|
|
Packit |
5c3484 |
adc $0, X0
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mov 8(rp,i,8), w2
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w2, w3
|
|
Packit |
5c3484 |
adc %rax, w0
|
|
Packit |
5c3484 |
MOV( %rdx, w1, 16)
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w3, X1
|
|
Packit |
5c3484 |
mov X1, 8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 32)
|
|
Packit |
5c3484 |
mov 16(rp,i,8), w3
|
|
Packit |
5c3484 |
adc $0, X1
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w3, w0
|
|
Packit |
5c3484 |
MOV( %rdx, w2, 64)
|
|
Packit |
5c3484 |
adc %rax, w1
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w0, X0
|
|
Packit |
5c3484 |
mov X0, 16(rp,i,8)
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 128)
|
|
Packit |
5c3484 |
adc %rax, X1
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
mov 24(rp,i,8), w0
|
|
Packit |
5c3484 |
adc $0, X0
|
|
Packit |
5c3484 |
add $4, i
|
|
Packit |
5c3484 |
jnc L(am2top1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, w1
|
|
Packit |
5c3484 |
adc w2, %rax
|
|
Packit |
5c3484 |
adc Z(i,$0), %rdx
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
adc Z(i,$0), X0
|
|
Packit |
5c3484 |
mov X1, I(-8(rp),-8(rp,i,8))
|
|
Packit |
5c3484 |
add X0, %rax
|
|
Packit |
5c3484 |
mov %rax, I((rp),(rp,i,8))
|
|
Packit |
5c3484 |
adc Z(i,$0), %rdx
|
|
Packit |
5c3484 |
mov %rdx, I(8(rp),8(rp,i,8))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
add $2, n
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(am2o3):
|
|
Packit |
5c3484 |
mov -8(up,n,8), v0
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mov %rax, v1
|
|
Packit |
5c3484 |
lea -1(n), i
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, X1
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 8)
|
|
Packit |
5c3484 |
mov (rp,n,8), w3
|
|
Packit |
5c3484 |
xor R32(w0), R32(w0)
|
|
Packit |
5c3484 |
xor R32(w1), R32(w1)
|
|
Packit |
5c3484 |
mov 8(up,n,8), %rax
|
|
Packit |
5c3484 |
jmp L(lo3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGNx
|
|
Packit |
5c3484 |
L(am2top3):
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, w1
|
|
Packit |
5c3484 |
adc %rax, w2
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
MOV( %rdx, w3, 1)
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
mov X1, -8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 2)
|
|
Packit |
5c3484 |
adc $0, X1
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
MOV( %rdx, w0, 4)
|
|
Packit |
5c3484 |
mov (rp,i,8), w1
|
|
Packit |
5c3484 |
add w1, w2
|
|
Packit |
5c3484 |
adc %rax, w3
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w2, X0
|
|
Packit |
5c3484 |
adc %rax, X1
|
|
Packit |
5c3484 |
mov X0, (rp,i,8)
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 8)
|
|
Packit |
5c3484 |
adc $0, X0
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mov 8(rp,i,8), w2
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w2, w3
|
|
Packit |
5c3484 |
adc %rax, w0
|
|
Packit |
5c3484 |
MOV( %rdx, w1, 16)
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
L(lo3): mul v0
|
|
Packit |
5c3484 |
add w3, X1
|
|
Packit |
5c3484 |
mov X1, 8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 32)
|
|
Packit |
5c3484 |
mov 16(rp,i,8), w3
|
|
Packit |
5c3484 |
adc $0, X1
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w3, w0
|
|
Packit |
5c3484 |
MOV( %rdx, w2, 64)
|
|
Packit |
5c3484 |
adc %rax, w1
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w0, X0
|
|
Packit |
5c3484 |
mov X0, 16(rp,i,8)
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 128)
|
|
Packit |
5c3484 |
adc %rax, X1
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
mov 24(rp,i,8), w0
|
|
Packit |
5c3484 |
adc $0, X0
|
|
Packit |
5c3484 |
add $4, i
|
|
Packit |
5c3484 |
jnc L(am2top3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, w1
|
|
Packit |
5c3484 |
adc w2, %rax
|
|
Packit |
5c3484 |
adc Z(i,$0), %rdx
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
adc Z(i,$0), X0
|
|
Packit |
5c3484 |
mov X1, I(-8(rp),-8(rp,i,8))
|
|
Packit |
5c3484 |
add X0, %rax
|
|
Packit |
5c3484 |
mov %rax, I((rp),(rp,i,8))
|
|
Packit |
5c3484 |
adc Z(i,$0), %rdx
|
|
Packit |
5c3484 |
mov %rdx, I(8(rp),8(rp,i,8))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
add $2, n
|
|
Packit |
5c3484 |
cmp $-1, n
|
|
Packit |
5c3484 |
jnz L(am2o1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(cor1):pop n
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
mov -16(up), v0
|
|
Packit |
5c3484 |
mov -8(up), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w3, %rax
|
|
Packit |
5c3484 |
adc $0, %rdx
|
|
Packit |
5c3484 |
mov %rax, -8(rp)
|
|
Packit |
5c3484 |
mov %rdx, (rp)
|
|
Packit |
5c3484 |
jmp L(sqr_diag_addlsh1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGNx
|
|
Packit |
5c3484 |
L(m2top2):
|
|
Packit |
5c3484 |
L(m2e2):mul v0
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov -8(up,i,8), %rax
|
|
Packit |
5c3484 |
mov w3, -8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
adc $0, R32(w1)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
adc %rdx, w1
|
|
Packit |
5c3484 |
mov $0, R32(w2)
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
mov w0, (rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w1
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
adc $0, R32(w2)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w1
|
|
Packit |
5c3484 |
adc %rdx, w2
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov $0, R32(w3)
|
|
Packit |
5c3484 |
add %rax, w1
|
|
Packit |
5c3484 |
adc %rdx, w2
|
|
Packit |
5c3484 |
adc $0, R32(w3)
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
mov w1, 8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w3
|
|
Packit |
5c3484 |
mov $0, R32(w0)
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
adc %rdx, w3
|
|
Packit |
5c3484 |
adc $0, R32(w0)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
mov $0, R32(w1)
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
mov w2, 16(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
add $4, i
|
|
Packit |
5c3484 |
js L(m2top2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov I(-8(up),-8(up,i,8)), %rax
|
|
Packit |
5c3484 |
mov w3, I(-8(rp),-8(rp,i,8))
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
adc R32(w1), R32(w1)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, %rax
|
|
Packit |
5c3484 |
adc w1, %rdx
|
|
Packit |
5c3484 |
mov %rax, I((rp),(rp,i,8))
|
|
Packit |
5c3484 |
mov %rdx, I(8(rp),8(rp,i,8))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
add $2, n C decrease |n|
|
|
Packit |
5c3484 |
jmp L(am2o0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGNx
|
|
Packit |
5c3484 |
L(m2top0):
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov -8(up,i,8), %rax
|
|
Packit |
5c3484 |
mov w3, -8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
adc $0, R32(w1)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
adc %rdx, w1
|
|
Packit |
5c3484 |
mov $0, R32(w2)
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
mov w0, (rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w1
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
adc $0, R32(w2)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w1
|
|
Packit |
5c3484 |
adc %rdx, w2
|
|
Packit |
5c3484 |
L(m2e0):mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov $0, R32(w3)
|
|
Packit |
5c3484 |
add %rax, w1
|
|
Packit |
5c3484 |
adc %rdx, w2
|
|
Packit |
5c3484 |
adc $0, R32(w3)
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
mov w1, 8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w3
|
|
Packit |
5c3484 |
mov $0, R32(w0)
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
adc %rdx, w3
|
|
Packit |
5c3484 |
adc $0, R32(w0)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
mov $0, R32(w1)
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
mov w2, 16(rp,i,8)
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
add $4, i
|
|
Packit |
5c3484 |
js L(m2top0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov I(-8(up),-8(up,i,8)), %rax
|
|
Packit |
5c3484 |
mov w3, I(-8(rp),-8(rp,i,8))
|
|
Packit |
5c3484 |
adc %rdx, w0
|
|
Packit |
5c3484 |
adc R32(w1), R32(w1)
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, %rax
|
|
Packit |
5c3484 |
adc w1, %rdx
|
|
Packit |
5c3484 |
mov %rax, I((rp),(rp,i,8))
|
|
Packit |
5c3484 |
mov %rdx, I(8(rp),8(rp,i,8))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
add $2, n C decrease |n|
|
|
Packit |
5c3484 |
cmp $-2, n
|
|
Packit |
5c3484 |
jz L(cor2) C jumps iff entry n = 5
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(am2o2):
|
|
Packit |
5c3484 |
mov -8(up,n,8), v0
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mov %rax, v1
|
|
Packit |
5c3484 |
lea -2(n), i
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 32)
|
|
Packit |
5c3484 |
mov (rp,n,8), w0
|
|
Packit |
5c3484 |
xor R32(w1), R32(w1)
|
|
Packit |
5c3484 |
xor R32(w2), R32(w2)
|
|
Packit |
5c3484 |
mov 8(up,n,8), %rax
|
|
Packit |
5c3484 |
jmp L(lo2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGNx
|
|
Packit |
5c3484 |
L(am2top2):
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, w1
|
|
Packit |
5c3484 |
adc %rax, w2
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
MOV( %rdx, w3, 1)
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
mov X1, -8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 2)
|
|
Packit |
5c3484 |
adc $0, X1
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
MOV( %rdx, w0, 4)
|
|
Packit |
5c3484 |
mov (rp,i,8), w1
|
|
Packit |
5c3484 |
add w1, w2
|
|
Packit |
5c3484 |
adc %rax, w3
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w2, X0
|
|
Packit |
5c3484 |
adc %rax, X1
|
|
Packit |
5c3484 |
mov X0, (rp,i,8)
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 8)
|
|
Packit |
5c3484 |
adc $0, X0
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mov 8(rp,i,8), w2
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w2, w3
|
|
Packit |
5c3484 |
adc %rax, w0
|
|
Packit |
5c3484 |
MOV( %rdx, w1, 16)
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w3, X1
|
|
Packit |
5c3484 |
mov X1, 8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 32)
|
|
Packit |
5c3484 |
mov 16(rp,i,8), w3
|
|
Packit |
5c3484 |
adc $0, X1
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w3, w0
|
|
Packit |
5c3484 |
MOV( %rdx, w2, 64)
|
|
Packit |
5c3484 |
adc %rax, w1
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
L(lo2): mul v0
|
|
Packit |
5c3484 |
add w0, X0
|
|
Packit |
5c3484 |
mov X0, 16(rp,i,8)
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 128)
|
|
Packit |
5c3484 |
adc %rax, X1
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
mov 24(rp,i,8), w0
|
|
Packit |
5c3484 |
adc $0, X0
|
|
Packit |
5c3484 |
add $4, i
|
|
Packit |
5c3484 |
jnc L(am2top2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, w1
|
|
Packit |
5c3484 |
adc w2, %rax
|
|
Packit |
5c3484 |
adc Z(i,$0), %rdx
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
adc Z(i,$0), X0
|
|
Packit |
5c3484 |
mov X1, I(-8(rp),-8(rp,i,8))
|
|
Packit |
5c3484 |
add X0, %rax
|
|
Packit |
5c3484 |
mov %rax, I((rp),(rp,i,8))
|
|
Packit |
5c3484 |
adc Z(i,$0), %rdx
|
|
Packit |
5c3484 |
mov %rdx, I(8(rp),8(rp,i,8))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
add $2, n
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(am2o0):
|
|
Packit |
5c3484 |
mov -8(up,n,8), v0
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mov %rax, v1
|
|
Packit |
5c3484 |
lea 0(n), i
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 2)
|
|
Packit |
5c3484 |
xor R32(w0), R32(w0)
|
|
Packit |
5c3484 |
mov (rp,n,8), w2
|
|
Packit |
5c3484 |
xor R32(w3), R32(w3)
|
|
Packit |
5c3484 |
jmp L(lo0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGNx
|
|
Packit |
5c3484 |
L(am2top0):
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, w1
|
|
Packit |
5c3484 |
adc %rax, w2
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
MOV( %rdx, w3, 1)
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
mov X1, -8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 2)
|
|
Packit |
5c3484 |
adc $0, X1
|
|
Packit |
5c3484 |
mov (up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
MOV( %rdx, w0, 4)
|
|
Packit |
5c3484 |
mov (rp,i,8), w1
|
|
Packit |
5c3484 |
add w1, w2
|
|
Packit |
5c3484 |
adc %rax, w3
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
L(lo0): mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w2, X0
|
|
Packit |
5c3484 |
adc %rax, X1
|
|
Packit |
5c3484 |
mov X0, (rp,i,8)
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 8)
|
|
Packit |
5c3484 |
adc $0, X0
|
|
Packit |
5c3484 |
mov 8(up,i,8), %rax
|
|
Packit |
5c3484 |
mov 8(rp,i,8), w2
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w2, w3
|
|
Packit |
5c3484 |
adc %rax, w0
|
|
Packit |
5c3484 |
MOV( %rdx, w1, 16)
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w3, X1
|
|
Packit |
5c3484 |
mov X1, 8(rp,i,8)
|
|
Packit |
5c3484 |
adc %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 32)
|
|
Packit |
5c3484 |
mov 16(rp,i,8), w3
|
|
Packit |
5c3484 |
adc $0, X1
|
|
Packit |
5c3484 |
mov 16(up,i,8), %rax
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w3, w0
|
|
Packit |
5c3484 |
MOV( %rdx, w2, 64)
|
|
Packit |
5c3484 |
adc %rax, w1
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
adc $0, w2
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w0, X0
|
|
Packit |
5c3484 |
mov X0, 16(rp,i,8)
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 128)
|
|
Packit |
5c3484 |
adc %rax, X1
|
|
Packit |
5c3484 |
mov 24(up,i,8), %rax
|
|
Packit |
5c3484 |
mov 24(rp,i,8), w0
|
|
Packit |
5c3484 |
adc $0, X0
|
|
Packit |
5c3484 |
add $4, i
|
|
Packit |
5c3484 |
jnc L(am2top0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, w1
|
|
Packit |
5c3484 |
adc w2, %rax
|
|
Packit |
5c3484 |
adc Z(i,$0), %rdx
|
|
Packit |
5c3484 |
add w1, X1
|
|
Packit |
5c3484 |
adc Z(i,$0), X0
|
|
Packit |
5c3484 |
mov X1, I(-8(rp),-8(rp,i,8))
|
|
Packit |
5c3484 |
add X0, %rax
|
|
Packit |
5c3484 |
mov %rax, I((rp),(rp,i,8))
|
|
Packit |
5c3484 |
adc Z(i,$0), %rdx
|
|
Packit |
5c3484 |
mov %rdx, I(8(rp),8(rp,i,8))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
add $2, n
|
|
Packit |
5c3484 |
cmp $-2, n
|
|
Packit |
5c3484 |
jnz L(am2o2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(cor2):pop n
|
|
Packit |
5c3484 |
mov -24(up), v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w0
|
|
Packit |
5c3484 |
mov -16(up), %rax
|
|
Packit |
5c3484 |
mov %rax, v1
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, X0
|
|
Packit |
5c3484 |
MOV( %rdx, X1, 32)
|
|
Packit |
5c3484 |
mov -8(up), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w2, X0
|
|
Packit |
5c3484 |
mov X0, -16(rp)
|
|
Packit |
5c3484 |
MOV( %rdx, X0, 128)
|
|
Packit |
5c3484 |
adc %rax, X1
|
|
Packit |
5c3484 |
mov -8(up), %rax
|
|
Packit |
5c3484 |
adc $0, X0
|
|
Packit |
5c3484 |
mul v1
|
|
Packit |
5c3484 |
add w0, X1
|
|
Packit |
5c3484 |
adc $0, X0
|
|
Packit |
5c3484 |
mov X1, -8(rp)
|
|
Packit |
5c3484 |
add X0, %rax
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
adc $0, %rdx
|
|
Packit |
5c3484 |
mov %rdx, 8(rp)
|
|
Packit |
5c3484 |
lea 8(rp), rp
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(sqr_diag_addlsh1):
|
|
Packit |
5c3484 |
mov -8(up,n,8), %rax
|
|
Packit |
5c3484 |
shl n
|
|
Packit |
5c3484 |
xor R32(%rbx), R32(%rbx)
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
mov 8(rp,n,8), %r11
|
|
Packit |
5c3484 |
lea (%rdx), %r10
|
|
Packit |
5c3484 |
mov 16(rp,n,8), %r9
|
|
Packit |
5c3484 |
add %r11, %r11
|
|
Packit |
5c3484 |
jmp L(dm)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGNx
|
|
Packit |
5c3484 |
L(dtop):mul %rax
|
|
Packit |
5c3484 |
add %r11, %r10
|
|
Packit |
5c3484 |
mov 8(rp,n,8), %r11
|
|
Packit |
5c3484 |
mov %r10, -8(rp,n,8)
|
|
Packit |
5c3484 |
adc %r9, %rax
|
|
Packit |
5c3484 |
lea (%rdx,%rbx), %r10
|
|
Packit |
5c3484 |
mov 16(rp,n,8), %r9
|
|
Packit |
5c3484 |
adc %r11, %r11
|
|
Packit |
5c3484 |
L(dm): mov %rax, (rp,n,8)
|
|
Packit |
5c3484 |
mov (up,n,4), %rax
|
|
Packit |
5c3484 |
adc %r9, %r9
|
|
Packit |
5c3484 |
setc R8(%rbx)
|
|
Packit |
5c3484 |
add $2, n
|
|
Packit |
5c3484 |
js L(dtop)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
add %r11, %r10
|
|
Packit |
5c3484 |
mov %r10, -8(rp)
|
|
Packit |
5c3484 |
adc %r9, %rax
|
|
Packit |
5c3484 |
lea (%rdx,%rbx), %r10
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
adc $0, %r10
|
|
Packit |
5c3484 |
mov %r10, 8(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 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(small):
|
|
Packit |
5c3484 |
mov (up), %rax
|
|
Packit |
5c3484 |
cmp $2, n_param
|
|
Packit |
5c3484 |
jae L(gt1)
|
|
Packit |
5c3484 |
L(n1):
|
|
Packit |
5c3484 |
mul %rax
|
|
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 |
L(n2): mov %rax, %r8
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
mov 8(up), %r11
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
mov %r11, %rax
|
|
Packit |
5c3484 |
mov %rdx, %r9
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
mov %rax, %r10
|
|
Packit |
5c3484 |
mov %r11, %rax
|
|
Packit |
5c3484 |
mov %rdx, %r11
|
|
Packit |
5c3484 |
mul %r8
|
|
Packit |
5c3484 |
xor %r8, %r8
|
|
Packit |
5c3484 |
add %rax, %r9
|
|
Packit |
5c3484 |
adc %rdx, %r10
|
|
Packit |
5c3484 |
adc %r8, %r11
|
|
Packit |
5c3484 |
add %rax, %r9
|
|
Packit |
5c3484 |
mov %r9, 8(rp)
|
|
Packit |
5c3484 |
adc %rdx, %r10
|
|
Packit |
5c3484 |
mov %r10, 16(rp)
|
|
Packit |
5c3484 |
adc %r8, %r11
|
|
Packit |
5c3484 |
mov %r11, 24(rp)
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(gt2):
|
|
Packit |
5c3484 |
L(n3): mov %rax, %r10
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
mov 8(up), %r11
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
mov %r11, %rax
|
|
Packit |
5c3484 |
mov %rdx, 8(rp)
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
mov 16(up), %rcx
|
|
Packit |
5c3484 |
mov %rax, 16(rp)
|
|
Packit |
5c3484 |
mov %rcx, %rax
|
|
Packit |
5c3484 |
mov %rdx, 24(rp)
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
mov %rax, 32(rp)
|
|
Packit |
5c3484 |
mov %rdx, 40(rp)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov %r11, %rax
|
|
Packit |
5c3484 |
mul %r10
|
|
Packit |
5c3484 |
mov %rax, %r8
|
|
Packit |
5c3484 |
mov %rcx, %rax
|
|
Packit |
5c3484 |
mov %rdx, %r9
|
|
Packit |
5c3484 |
mul %r10
|
|
Packit |
5c3484 |
xor %r10, %r10
|
|
Packit |
5c3484 |
add %rax, %r9
|
|
Packit |
5c3484 |
mov %r11, %rax
|
|
Packit |
5c3484 |
mov %r10, %r11
|
|
Packit |
5c3484 |
adc %rdx, %r10
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul %rcx
|
|
Packit |
5c3484 |
add %rax, %r10
|
|
Packit |
5c3484 |
adc %r11, %rdx
|
|
Packit |
5c3484 |
add %r8, %r8
|
|
Packit |
5c3484 |
adc %r9, %r9
|
|
Packit |
5c3484 |
adc %r10, %r10
|
|
Packit |
5c3484 |
adc %rdx, %rdx
|
|
Packit |
5c3484 |
adc %r11, %r11
|
|
Packit |
5c3484 |
add %r8, 8(rp)
|
|
Packit |
5c3484 |
adc %r9, 16(rp)
|
|
Packit |
5c3484 |
adc %r10, 24(rp)
|
|
Packit |
5c3484 |
adc %rdx, 32(rp)
|
|
Packit |
5c3484 |
adc %r11, 40(rp)
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
EPILOGUE()
|
|
rpm-build |
c3cd4f |
CF_PROT
|