| Index: gcc/gmp/mpn/x86_64/popham.asm
|
| diff --git a/gcc/gmp/mpn/x86_64/popham.asm b/gcc/gmp/mpn/x86_64/popham.asm
|
| deleted file mode 100644
|
| index e2bdb1a0b882b28b68b2535789ffd5d705b14c2f..0000000000000000000000000000000000000000
|
| --- a/gcc/gmp/mpn/x86_64/popham.asm
|
| +++ /dev/null
|
| @@ -1,157 +0,0 @@
|
| -dnl AMD64 mpn_popcount, mpn_hamdist -- population count and hamming distance.
|
| -
|
| -dnl Copyright 2004, 2005, 2007 Free Software Foundation, Inc.
|
| -
|
| -dnl This file is part of the GNU MP Library.
|
| -
|
| -dnl The GNU MP Library is free software; you can redistribute it and/or modify
|
| -dnl it under the terms of the GNU Lesser General Public License as published
|
| -dnl by the Free Software Foundation; either version 3 of the License, or (at
|
| -dnl your option) any later version.
|
| -
|
| -dnl The GNU MP Library is distributed in the hope that it will be useful, but
|
| -dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
|
| -dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
|
| -dnl License for more details.
|
| -
|
| -dnl You should have received a copy of the GNU Lesser General Public License
|
| -dnl along with the GNU MP Library. If not, see http://www.gnu.org/licenses/.
|
| -
|
| -
|
| -include(`../config.m4')
|
| -
|
| -
|
| -C popcount hamdist
|
| -C cycles/limb cycles/limb
|
| -C K8,K9: 6 7
|
| -C K10: 6 7
|
| -C P4: 12 14.3
|
| -C P6-15: 7 8
|
| -
|
| -C TODO
|
| -C * Tune. It should be possible to reach 5 c/l for popcount and 6 c/l for
|
| -C hamdist for K8/K9.
|
| -
|
| -
|
| -ifdef(`OPERATION_popcount',`
|
| - define(`func',`mpn_popcount')
|
| - define(`up', `%rdi')
|
| - define(`n', `%rsi')
|
| - define(`h55555555', `%r10')
|
| - define(`h33333333', `%r11')
|
| - define(`h0f0f0f0f', `%rcx')
|
| - define(`h01010101', `%rdx')
|
| - define(`HAM', `dnl')
|
| -')
|
| -ifdef(`OPERATION_hamdist',`
|
| - define(`func',`mpn_hamdist')
|
| - define(`up', `%rdi')
|
| - define(`vp', `%rsi')
|
| - define(`n', `%rdx')
|
| - define(`h55555555', `%r10')
|
| - define(`h33333333', `%r11')
|
| - define(`h0f0f0f0f', `%rcx')
|
| - define(`h01010101', `%r14')
|
| - define(`HAM', `$1')
|
| -')
|
| -
|
| -
|
| -MULFUNC_PROLOGUE(mpn_popcount mpn_hamdist)
|
| -
|
| -ASM_START()
|
| - TEXT
|
| - ALIGN(32)
|
| -PROLOGUE(func)
|
| -
|
| - pushq %r12
|
| - pushq %r13
|
| - HAM(` pushq %r14 ')
|
| -
|
| - movq $0x5555555555555555, h55555555
|
| - movq $0x3333333333333333, h33333333
|
| - movq $0x0f0f0f0f0f0f0f0f, h0f0f0f0f
|
| - movq $0x0101010101010101, h01010101
|
| -
|
| - leaq (up,n,8), up
|
| - HAM(` leaq (vp,n,8), vp ')
|
| - negq n
|
| -
|
| - xorl %eax, %eax
|
| -
|
| - btq $0, n
|
| - jnc L(oop)
|
| -
|
| - movq (up,n,8), %r8
|
| - HAM(` xorq (vp,n,8), %r8 ')
|
| -
|
| - movq %r8, %r9
|
| - shrq %r8
|
| - andq h55555555, %r8
|
| - subq %r8, %r9
|
| -
|
| - movq %r9, %r8
|
| - shrq $2, %r9
|
| - andq h33333333, %r8
|
| - andq h33333333, %r9
|
| - addq %r8, %r9 C 16 4-bit fields (0..4)
|
| -
|
| - movq %r9, %r8
|
| - shrq $4, %r9
|
| - andq h0f0f0f0f, %r8
|
| - andq h0f0f0f0f, %r9
|
| - addq %r8, %r9 C 8 8-bit fields (0..16)
|
| -
|
| - imulq h01010101, %r9 C sum the 8 fields in high 8 bits
|
| - shrq $56, %r9
|
| -
|
| - addq %r9, %rax C add to total
|
| - addq $1, n
|
| - jz L(done)
|
| -
|
| - ALIGN(16)
|
| -L(oop): movq (up,n,8), %r8
|
| - movq 8(up,n,8), %r12
|
| - HAM(` xorq (vp,n,8), %r8 ')
|
| - HAM(` xorq 8(vp,n,8), %r12 ')
|
| -
|
| - movq %r8, %r9
|
| - movq %r12, %r13
|
| - shrq %r8
|
| - shrq %r12
|
| - andq h55555555, %r8
|
| - andq h55555555, %r12
|
| - subq %r8, %r9
|
| - subq %r12, %r13
|
| -
|
| - movq %r9, %r8
|
| - movq %r13, %r12
|
| - shrq $2, %r9
|
| - shrq $2, %r13
|
| - andq h33333333, %r8
|
| - andq h33333333, %r9
|
| - andq h33333333, %r12
|
| - andq h33333333, %r13
|
| - addq %r8, %r9 C 16 4-bit fields (0..4)
|
| - addq %r12, %r13 C 16 4-bit fields (0..4)
|
| -
|
| - addq %r13, %r9 C 16 4-bit fields (0..8)
|
| - movq %r9, %r8
|
| - shrq $4, %r9
|
| - andq h0f0f0f0f, %r8
|
| - andq h0f0f0f0f, %r9
|
| - addq %r8, %r9 C 8 8-bit fields (0..16)
|
| -
|
| - imulq h01010101, %r9 C sum the 8 fields in high 8 bits
|
| - shrq $56, %r9
|
| -
|
| - addq %r9, %rax C add to total
|
| - addq $2, n
|
| - jnc L(oop)
|
| -
|
| -L(done):
|
| - HAM(` popq %r14 ')
|
| - popq %r13
|
| - popq %r12
|
| - ret
|
| -
|
| -EPILOGUE()
|
|
|