|
Packit |
5c3484 |
dnl AMD64 mpn_copyi optimised for CPUs with fast SSE.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Copyright 2003, 2005, 2007, 2011, 2012, 2015 Free Software Foundation,
|
|
Packit |
5c3484 |
dnl Inc.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Contributed to the GNU project by Torbjörn Granlund.
|
|
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 cycles/limb cycles/limb good
|
|
Packit |
5c3484 |
C aligned unaligned best seen for cpu?
|
|
Packit |
5c3484 |
C AMD K8,K9
|
|
Packit |
5c3484 |
C AMD K10 0.85 1.64 Y/N
|
|
Packit |
5c3484 |
C AMD bull 1.4 1.4 N
|
|
Packit |
5c3484 |
C AMD pile 0.77 0.93 N
|
|
Packit |
5c3484 |
C AMD steam ? ?
|
|
Packit |
5c3484 |
C AMD excavator ? ?
|
|
Packit |
5c3484 |
C AMD bobcat
|
|
Packit |
5c3484 |
C AMD jaguar 0.65 1.02 opt/0.93 Y/N
|
|
Packit |
5c3484 |
C Intel P4 2.3 2.3 Y
|
|
Packit |
5c3484 |
C Intel core 1.0 1.0 0.52/0.64 N
|
|
Packit |
5c3484 |
C Intel NHM 0.5 0.67 Y
|
|
Packit |
5c3484 |
C Intel SBR 0.51 0.75 opt/0.54 Y/N
|
|
Packit |
5c3484 |
C Intel IBR 0.50 0.57 opt/0.54 Y
|
|
Packit |
5c3484 |
C Intel HWL 0.50 0.57 opt/0.51 Y
|
|
Packit |
5c3484 |
C Intel BWL 0.55 0.62 opt/0.55 Y
|
|
Packit |
5c3484 |
C Intel atom
|
|
Packit |
5c3484 |
C Intel SLM 1.02 1.27 opt/1.07 Y/N
|
|
Packit |
5c3484 |
C VIA nano 1.16 5.16 Y/N
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C We try to do as many 16-byte operations as possible. The top-most and
|
|
Packit |
5c3484 |
C bottom-most writes might need 8-byte operations. We can always write using
|
|
Packit |
5c3484 |
C aligned 16-byte operations, we read with both aligned and unaligned 16-byte
|
|
Packit |
5c3484 |
C operations.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C Instead of having separate loops for reading aligned and unaligned, we read
|
|
Packit |
5c3484 |
C using MOVDQU. This seems to work great except for core2; there performance
|
|
Packit |
5c3484 |
C doubles when reading using MOVDQA (for aligned source). It is unclear how to
|
|
Packit |
5c3484 |
C best handle the unaligned case there.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C INPUT PARAMETERS
|
|
Packit |
5c3484 |
define(`rp', `%rdi')
|
|
Packit |
5c3484 |
define(`up', `%rsi')
|
|
Packit |
5c3484 |
define(`n', `%rdx')
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ABI_SUPPORT(DOS64)
|
|
Packit |
5c3484 |
ABI_SUPPORT(STD64)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl define(`movdqu', lddqu)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ASM_START()
|
|
Packit |
5c3484 |
TEXT
|
|
Packit |
5c3484 |
ALIGN(64)
|
|
Packit |
5c3484 |
PROLOGUE(mpn_copyi)
|
|
Packit |
5c3484 |
FUNC_ENTRY(3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
cmp $3, n C NB: bc code below assumes this limit
|
|
Packit |
5c3484 |
jc L(bc)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
test $8, R8(rp) C is rp 16-byte aligned?
|
|
Packit |
5c3484 |
jz L(ali) C jump if rp aligned
|
|
Packit |
5c3484 |
movsq C copy single limb
|
|
Packit |
5c3484 |
dec n
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
sub $16, n
|
|
Packit |
5c3484 |
jc L(sma)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(top): movdqu (up), %xmm0
|
|
Packit |
5c3484 |
movdqu 16(up), %xmm1
|
|
Packit |
5c3484 |
movdqu 32(up), %xmm2
|
|
Packit |
5c3484 |
movdqu 48(up), %xmm3
|
|
Packit |
5c3484 |
movdqu 64(up), %xmm4
|
|
Packit |
5c3484 |
movdqu 80(up), %xmm5
|
|
Packit |
5c3484 |
movdqu 96(up), %xmm6
|
|
Packit |
5c3484 |
movdqu 112(up), %xmm7
|
|
Packit |
5c3484 |
lea 128(up), up
|
|
Packit |
5c3484 |
movdqa %xmm0, (rp)
|
|
Packit |
5c3484 |
movdqa %xmm1, 16(rp)
|
|
Packit |
5c3484 |
movdqa %xmm2, 32(rp)
|
|
Packit |
5c3484 |
movdqa %xmm3, 48(rp)
|
|
Packit |
5c3484 |
movdqa %xmm4, 64(rp)
|
|
Packit |
5c3484 |
movdqa %xmm5, 80(rp)
|
|
Packit |
5c3484 |
movdqa %xmm6, 96(rp)
|
|
Packit |
5c3484 |
movdqa %xmm7, 112(rp)
|
|
Packit |
5c3484 |
lea 128(rp), rp
|
|
Packit |
5c3484 |
L(ali): sub $16, n
|
|
Packit |
5c3484 |
jnc L(top)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(sma): test $8, R8(n)
|
|
Packit |
5c3484 |
jz 1f
|
|
Packit |
5c3484 |
movdqu (up), %xmm0
|
|
Packit |
5c3484 |
movdqu 16(up), %xmm1
|
|
Packit |
5c3484 |
movdqu 32(up), %xmm2
|
|
Packit |
5c3484 |
movdqu 48(up), %xmm3
|
|
Packit |
5c3484 |
lea 64(up), up
|
|
Packit |
5c3484 |
movdqa %xmm0, (rp)
|
|
Packit |
5c3484 |
movdqa %xmm1, 16(rp)
|
|
Packit |
5c3484 |
movdqa %xmm2, 32(rp)
|
|
Packit |
5c3484 |
movdqa %xmm3, 48(rp)
|
|
Packit |
5c3484 |
lea 64(rp), rp
|
|
Packit |
5c3484 |
1:
|
|
Packit |
5c3484 |
test $4, R8(n)
|
|
Packit |
5c3484 |
jz 1f
|
|
Packit |
5c3484 |
movdqu (up), %xmm0
|
|
Packit |
5c3484 |
movdqu 16(up), %xmm1
|
|
Packit |
5c3484 |
lea 32(up), up
|
|
Packit |
5c3484 |
movdqa %xmm0, (rp)
|
|
Packit |
5c3484 |
movdqa %xmm1, 16(rp)
|
|
Packit |
5c3484 |
lea 32(rp), rp
|
|
Packit |
5c3484 |
1:
|
|
Packit |
5c3484 |
test $2, R8(n)
|
|
Packit |
5c3484 |
jz 1f
|
|
Packit |
5c3484 |
movdqu (up), %xmm0
|
|
Packit |
5c3484 |
lea 16(up), up
|
|
Packit |
5c3484 |
movdqa %xmm0, (rp)
|
|
Packit |
5c3484 |
lea 16(rp), rp
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
1:
|
|
Packit |
5c3484 |
L(end): test $1, R8(n)
|
|
Packit |
5c3484 |
jz 1f
|
|
Packit |
5c3484 |
mov (up), %r8
|
|
Packit |
5c3484 |
mov %r8, (rp)
|
|
Packit |
5c3484 |
1:
|
|
Packit |
5c3484 |
FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C Basecase code. Needed for good small operands speed, not for correctness as
|
|
Packit |
5c3484 |
C the above code is currently written. The commented-out lines need to be
|
|
Packit |
5c3484 |
C reinstated if this code is to be used for n > 3, and then the post loop
|
|
Packit |
5c3484 |
C offsets need fixing.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(bc): sub $2, n
|
|
Packit |
5c3484 |
jc L(end)
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
1: mov (up), %rax
|
|
Packit |
5c3484 |
mov 8(up), %rcx
|
|
Packit |
5c3484 |
dnl lea 16(up), up
|
|
Packit |
5c3484 |
mov %rax, (rp)
|
|
Packit |
5c3484 |
mov %rcx, 8(rp)
|
|
Packit |
5c3484 |
dnl lea 16(rp), rp
|
|
Packit |
5c3484 |
dnl sub $2, n
|
|
Packit |
5c3484 |
dnl jnc 1b
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
test $1, R8(n)
|
|
Packit |
5c3484 |
jz L(ret)
|
|
Packit |
5c3484 |
mov 16(up), %rax
|
|
Packit |
5c3484 |
mov %rax, 16(rp)
|
|
Packit |
5c3484 |
L(ret): FUNC_EXIT()
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
EPILOGUE()
|