|
Packit |
5c3484 |
dnl AMD64 mpn_sqr_basecase optimised for AMD bobcat.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Copyright 2003-2005, 2007, 2008, 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 4.5
|
|
Packit |
5c3484 |
C AMD K10 4.5
|
|
Packit |
5c3484 |
C AMD bd1 4.75
|
|
Packit |
5c3484 |
C AMD bobcat 5
|
|
Packit |
5c3484 |
C Intel P4 17.7
|
|
Packit |
5c3484 |
C Intel core2 5.5
|
|
Packit |
5c3484 |
C Intel NHM 5.43
|
|
Packit |
5c3484 |
C Intel SBR 3.92
|
|
Packit |
5c3484 |
C Intel atom 23
|
|
Packit |
5c3484 |
C VIA nano 5.63
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C This sqr_basecase is based on mul_1 and addmul_1, since these both run at the
|
|
Packit |
5c3484 |
C multiply insn bandwidth, without any apparent loop branch exit pipeline
|
|
Packit |
5c3484 |
C replays experienced on K8. The structure is unusual: it falls into mul_1 in
|
|
Packit |
5c3484 |
C the same way for all n, then it splits into 4 different wind-down blocks and
|
|
Packit |
5c3484 |
C 4 separate addmul_1 loops.
|
|
Packit |
5c3484 |
C
|
|
Packit |
5c3484 |
C We have not tried using the same addmul_1 loops with a switch into feed-in
|
|
Packit |
5c3484 |
C code, as we do in other basecase implementations. Doing that could save
|
|
Packit |
5c3484 |
C substantial code volume, but would also probably add some overhead.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C TODO
|
|
Packit |
5c3484 |
C * Tune un < 4 code.
|
|
Packit |
5c3484 |
C * Perhaps implement a larger final corner (it is now 2 x 1).
|
|
Packit |
5c3484 |
C * Lots of space could be saved by replacing the "switch" code by gradual
|
|
Packit |
5c3484 |
C jumps out from mul_1 winddown code, perhaps with no added overhead.
|
|
Packit |
5c3484 |
C * Are the ALIGN(16) really necessary? They add about 25 bytes of padding.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ABI_SUPPORT(DOS64)
|
|
Packit |
5c3484 |
ABI_SUPPORT(STD64)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C Standard parameters
|
|
Packit |
5c3484 |
define(`rp', `%rdi')
|
|
Packit |
5c3484 |
define(`up', `%rsi')
|
|
Packit |
5c3484 |
define(`un_param', `%rdx')
|
|
Packit |
5c3484 |
C Standard allocations
|
|
Packit |
5c3484 |
define(`un', `%rbx')
|
|
Packit |
5c3484 |
define(`w0', `%r8')
|
|
Packit |
5c3484 |
define(`w1', `%r9')
|
|
Packit |
5c3484 |
define(`w2', `%r10')
|
|
Packit |
5c3484 |
define(`w3', `%r11')
|
|
Packit |
5c3484 |
define(`n', `%rbp')
|
|
Packit |
5c3484 |
define(`v0', `%rcx')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C Temp macro for allowing control over indexing.
|
|
Packit |
5c3484 |
C Define to return $1 for more conservative ptr handling.
|
|
Packit |
5c3484 |
define(`X',`$2')
|
|
Packit |
5c3484 |
dnl define(`X',`$1')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ASM_START()
|
|
Packit |
5c3484 |
TEXT
|
|
Packit |
5c3484 |
ALIGN(64)
|
|
Packit |
5c3484 |
PROLOGUE(mpn_sqr_basecase)
|
|
Packit |
5c3484 |
FUNC_ENTRY(3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov (up), %rax
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
cmp $2, R32(un_param)
|
|
Packit |
5c3484 |
jae L(ge2)
|
|
Packit |
5c3484 |
|
|
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(ge2): mov (up), v0
|
|
Packit |
5c3484 |
jnz L(g2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
mov 8(up), %rax
|
|
Packit |
5c3484 |
mov %rdx, w0
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 8(up), v0
|
|
Packit |
5c3484 |
mov (up), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w0
|
|
Packit |
5c3484 |
mov w0, 8(rp)
|
|
Packit |
5c3484 |
mov %rdx, w0 C CAUTION: r8 realloc
|
|
Packit |
5c3484 |
adc $0, w0
|
|
Packit |
5c3484 |
mov 8(up), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add w1, w0
|
|
Packit |
5c3484 |
adc $0, %rdx
|
|
Packit |
5c3484 |
add w0, %rax
|
|
Packit |
5c3484 |
adc $0, %rdx
|
|
Packit |
5c3484 |
mov %rax, 16(rp)
|
|
Packit |
5c3484 |
mov %rdx, 24(rp)
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(g2): cmp $3, R32(un_param)
|
|
Packit |
5c3484 |
ja L(g3)
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
mov %rdx, 8(rp)
|
|
Packit |
5c3484 |
mov 8(up), %rax
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
mov %rax, 16(rp)
|
|
Packit |
5c3484 |
mov %rdx, 24(rp)
|
|
Packit |
5c3484 |
mov 16(up), %rax
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
mov %rax, 32(rp)
|
|
Packit |
5c3484 |
mov %rdx, 40(rp)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov (up), v0
|
|
Packit |
5c3484 |
mov 8(up), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
mov 16(up), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
xor R32(w2), R32(w2)
|
|
Packit |
5c3484 |
add %rax, w1
|
|
Packit |
5c3484 |
adc %rdx, w2
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov 8(up), v0
|
|
Packit |
5c3484 |
mov 16(up), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
xor R32(w3), R32(w3)
|
|
Packit |
5c3484 |
add %rax, w2
|
|
Packit |
5c3484 |
adc %rdx, w3
|
|
Packit |
5c3484 |
add w0, w0
|
|
Packit |
5c3484 |
adc w1, w1
|
|
Packit |
5c3484 |
adc w2, w2
|
|
Packit |
5c3484 |
adc w3, w3
|
|
Packit |
5c3484 |
mov $0, R32(v0)
|
|
Packit |
5c3484 |
adc v0, v0
|
|
Packit |
5c3484 |
add w0, 8(rp)
|
|
Packit |
5c3484 |
adc w1, 16(rp)
|
|
Packit |
5c3484 |
adc w2, 24(rp)
|
|
Packit |
5c3484 |
adc w3, 32(rp)
|
|
Packit |
5c3484 |
adc v0, 40(rp)
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(g3): push %rbx
|
|
Packit |
5c3484 |
push %rbp
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov 8(up), %rax
|
|
Packit |
5c3484 |
lea -24(rp,un_param,8), rp
|
|
Packit |
5c3484 |
lea -24(up,un_param,8), up
|
|
Packit |
5c3484 |
neg un_param
|
|
Packit |
5c3484 |
push un_param C for sqr_diag_addlsh1
|
|
Packit |
5c3484 |
lea (un_param), un
|
|
Packit |
5c3484 |
lea 3(un_param), n
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
jmp L(L3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(top): mov w0, -16(rp,n,8)
|
|
Packit |
5c3484 |
add w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
mov w2, -8(rp,n,8)
|
|
Packit |
5c3484 |
add w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 8(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
mov w0, (rp,n,8)
|
|
Packit |
5c3484 |
add w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
L(L3): mov 16(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
mov w2, 8(rp,n,8)
|
|
Packit |
5c3484 |
add w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 24(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add $4, n
|
|
Packit |
5c3484 |
js L(top)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov w0, -16(rp,n,8)
|
|
Packit |
5c3484 |
add w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
test n, n
|
|
Packit |
5c3484 |
jz L(r2)
|
|
Packit |
5c3484 |
cmp $2, R32(n)
|
|
Packit |
5c3484 |
ja L(r3)
|
|
Packit |
5c3484 |
jz L(r0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(r1): mov X((up,n,8),8(up)), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
mov w2, X(-8(rp,n,8),(rp))
|
|
Packit |
5c3484 |
add w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov X(8(up,n,8),16(up)), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
mov w0, X((rp,n,8),8(rp))
|
|
Packit |
5c3484 |
add w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov w2, X(8(rp,n,8),16(rp))
|
|
Packit |
5c3484 |
mov w3, X(16(rp,n,8),24(rp))
|
|
Packit |
5c3484 |
add $5, un
|
|
Packit |
5c3484 |
jmp L(to0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(r2): mov X((up,n,8),(up)), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
mov w2, X(-8(rp,n,8),-8(rp))
|
|
Packit |
5c3484 |
add w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov X(8(up,n,8),8(up)), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
mov w0, X((rp,n,8),(rp))
|
|
Packit |
5c3484 |
add w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov X(16(up,n,8),16(up)), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
mov w2, X(8(rp,n,8),8(rp))
|
|
Packit |
5c3484 |
add w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov w0, X(16(rp,n,8),16(rp))
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov w1, X(24(rp,n,8),24(rp))
|
|
Packit |
5c3484 |
add $6, un
|
|
Packit |
5c3484 |
jmp L(to1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(r3): mov w2, X(-8(rp,n,8),16(rp))
|
|
Packit |
5c3484 |
mov w3, X((rp,n,8),24(rp))
|
|
Packit |
5c3484 |
add $3, un
|
|
Packit |
5c3484 |
jmp L(to2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(r0): mov X((up,n,8),16(up)), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
mov w2, X(-8(rp,n,8),8(rp))
|
|
Packit |
5c3484 |
add w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov w0, X((rp,n,8),16(rp))
|
|
Packit |
5c3484 |
mov w1, X(8(rp,n,8),24(rp))
|
|
Packit |
5c3484 |
add $4, un
|
|
Packit |
5c3484 |
C jmp L(to3)
|
|
Packit |
5c3484 |
C fall through into main loop
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(outer):
|
|
Packit |
5c3484 |
mov un, n
|
|
Packit |
5c3484 |
mov (up,un,8), v0
|
|
Packit |
5c3484 |
mov 8(up,un,8), %rax
|
|
Packit |
5c3484 |
lea 8(rp), rp
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
jmp L(al3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(ta3): add w0, -16(rp,n,8)
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
add w2, -8(rp,n,8)
|
|
Packit |
5c3484 |
adc w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 8(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add w0, (rp,n,8)
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
L(al3): mov 16(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
add w2, 8(rp,n,8)
|
|
Packit |
5c3484 |
adc w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 24(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add $4, n
|
|
Packit |
5c3484 |
js L(ta3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
add w0, X(-16(rp,n,8),8(rp))
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add w2, X(-8(rp,n,8),16(rp))
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov w3, X((rp,n,8),24(rp))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(to2): mov un, n
|
|
Packit |
5c3484 |
cmp $-4, R32(un)
|
|
Packit |
5c3484 |
jnc L(end)
|
|
Packit |
5c3484 |
add $4, un
|
|
Packit |
5c3484 |
mov 8(up,n,8), v0
|
|
Packit |
5c3484 |
mov 16(up,n,8), %rax
|
|
Packit |
5c3484 |
lea 8(rp), rp
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
jmp L(al2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(ta2): add w0, -16(rp,n,8)
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
add w2, -8(rp,n,8)
|
|
Packit |
5c3484 |
adc w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 8(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add w0, (rp,n,8)
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov 16(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
add w2, 8(rp,n,8)
|
|
Packit |
5c3484 |
adc w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
L(al2): mov 24(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add $4, n
|
|
Packit |
5c3484 |
js L(ta2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
add w0, X(-16(rp,n,8),8(rp))
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add w2, X(-8(rp,n,8),16(rp))
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov w3, X((rp,n,8),24(rp))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(to1): mov un, n
|
|
Packit |
5c3484 |
mov -16(up,un,8), v0
|
|
Packit |
5c3484 |
mov -8(up,un,8), %rax
|
|
Packit |
5c3484 |
lea 8(rp), rp
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
jmp L(al1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(ta1): add w0, -16(rp,n,8)
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
L(al1): mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
add w2, -8(rp,n,8)
|
|
Packit |
5c3484 |
adc w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 8(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add w0, (rp,n,8)
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov 16(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
add w2, 8(rp,n,8)
|
|
Packit |
5c3484 |
adc w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 24(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add $4, n
|
|
Packit |
5c3484 |
js L(ta1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
add w0, X(-16(rp,n,8),8(rp))
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add w2, X(-8(rp,n,8),16(rp))
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov w3, X((rp,n,8),24(rp))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(to0): mov un, n
|
|
Packit |
5c3484 |
mov -8(up,un,8), v0
|
|
Packit |
5c3484 |
mov (up,un,8), %rax
|
|
Packit |
5c3484 |
lea 8(rp), rp
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
jmp L(al0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(ta0): add w0, -16(rp,n,8)
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov (up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
add w2, -8(rp,n,8)
|
|
Packit |
5c3484 |
adc w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
L(al0): mov 8(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add w0, (rp,n,8)
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov 16(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
add w2, 8(rp,n,8)
|
|
Packit |
5c3484 |
adc w3, w0
|
|
Packit |
5c3484 |
adc $0, w1
|
|
Packit |
5c3484 |
mov 24(up,n,8), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add $4, n
|
|
Packit |
5c3484 |
js L(ta0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
add w0, X(-16(rp,n,8),8(rp))
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add w2, X(-8(rp,n,8),16(rp))
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov w3, X((rp,n,8),24(rp))
|
|
Packit |
5c3484 |
jmp L(outer)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(end): mov X(8(up,un,8),(up)), v0
|
|
Packit |
5c3484 |
mov X(16(up,un,8),8(up)), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w0
|
|
Packit |
5c3484 |
mov %rdx, w1
|
|
Packit |
5c3484 |
mov X(24(up,un,8),16(up)), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
mov %rax, w2
|
|
Packit |
5c3484 |
mov %rdx, w3
|
|
Packit |
5c3484 |
add w0, X(24(rp,un,8),16(rp))
|
|
Packit |
5c3484 |
adc w1, w2
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
add w2, X(32(rp,un,8),24(rp))
|
|
Packit |
5c3484 |
adc $0, w3
|
|
Packit |
5c3484 |
mov X(16(up,un,8),8(up)), v0
|
|
Packit |
5c3484 |
mov X(24(up,un,8),16(up)), %rax
|
|
Packit |
5c3484 |
mul v0
|
|
Packit |
5c3484 |
add %rax, w3
|
|
Packit |
5c3484 |
mov w3, X(40(rp,un,8),32(rp))
|
|
Packit |
5c3484 |
adc $0, %rdx
|
|
Packit |
5c3484 |
mov %rdx, X(48(rp,un,8),40(rp))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C sqr_diag_addlsh1
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
lea 16(up), up
|
|
Packit |
5c3484 |
lea 40(rp), rp
|
|
Packit |
5c3484 |
pop n
|
|
Packit |
5c3484 |
lea 2(n,n), n
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov (up,n,4), %rax
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
xor R32(w2), R32(w2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov 8(rp,n,8), w0
|
|
Packit |
5c3484 |
mov %rax, (rp,n,8)
|
|
Packit |
5c3484 |
jmp L(lm)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(8)
|
|
Packit |
5c3484 |
L(tsd): add %rbx, w0
|
|
Packit |
5c3484 |
adc %rax, w1
|
|
Packit |
5c3484 |
mov w0, -8(rp,n,8)
|
|
Packit |
5c3484 |
mov 8(rp,n,8), w0
|
|
Packit |
5c3484 |
mov w1, (rp,n,8)
|
|
Packit |
5c3484 |
L(lm): mov 16(rp,n,8), w1
|
|
Packit |
5c3484 |
adc w0, w0
|
|
Packit |
5c3484 |
adc w1, w1
|
|
Packit |
5c3484 |
lea (%rdx,w2), %rbx
|
|
Packit |
5c3484 |
mov 8(up,n,4), %rax
|
|
Packit |
5c3484 |
setc R8(w2)
|
|
Packit |
5c3484 |
mul %rax
|
|
Packit |
5c3484 |
add $2, n
|
|
Packit |
5c3484 |
js L(tsd)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(esd): add %rbx, w0
|
|
Packit |
5c3484 |
adc %rax, w1
|
|
Packit |
5c3484 |
mov w0, X(-8(rp,n,8),-8(rp))
|
|
Packit |
5c3484 |
mov w1, X((rp,n,8),(rp))
|
|
Packit |
5c3484 |
adc w2, %rdx
|
|
Packit |
5c3484 |
mov %rdx, X(8(rp,n,8),8(rp))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
pop %rbp
|
|
Packit |
5c3484 |
pop %rbx
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
EPILOGUE()
|