|
Packit |
5c3484 |
dnl AMD K7 mpn_mul_1.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dnl Copyright 1999-2002, 2005, 2008 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 |
|
|
Packit |
5c3484 |
C cycles/limb
|
|
Packit |
5c3484 |
C P5
|
|
Packit |
5c3484 |
C P6 model 0-8,10-12)
|
|
Packit |
5c3484 |
C P6 model 9 (Banias)
|
|
Packit |
5c3484 |
C P6 model 13 (Dothan)
|
|
Packit |
5c3484 |
C P4 model 0 (Willamette)
|
|
Packit |
5c3484 |
C P4 model 1 (?)
|
|
Packit |
5c3484 |
C P4 model 2 (Northwood)
|
|
Packit |
5c3484 |
C P4 model 3 (Prescott)
|
|
Packit |
5c3484 |
C P4 model 4 (Nocona)
|
|
Packit |
5c3484 |
C AMD K6
|
|
Packit |
5c3484 |
C AMD K7 3.25
|
|
Packit |
5c3484 |
C AMD K8
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
C TODO
|
|
Packit |
5c3484 |
C * Improve feed-in and wind-down code. We beat the old code for all n != 1,
|
|
Packit |
5c3484 |
C but we might be able to do even better.
|
|
Packit |
5c3484 |
C * The feed-in code for mul_1c is crude.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ASM_START()
|
|
Packit |
5c3484 |
TEXT
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
PROLOGUE(mpn_mul_1c)
|
|
Packit |
5c3484 |
add $-16, %esp
|
|
Packit |
5c3484 |
mov %ebp, (%esp)
|
|
Packit |
5c3484 |
mov %ebx, 4(%esp)
|
|
Packit |
5c3484 |
mov %esi, 8(%esp)
|
|
Packit |
5c3484 |
mov %edi, 12(%esp)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov 20(%esp), %edi
|
|
Packit |
5c3484 |
mov 24(%esp), %esi
|
|
Packit |
5c3484 |
mov 28(%esp), %ebp
|
|
Packit |
5c3484 |
mov 32(%esp), %ecx
|
|
Packit |
5c3484 |
mov %ebp, %ebx
|
|
Packit |
5c3484 |
shr $2, %ebp
|
|
Packit |
5c3484 |
mov %ebp, 28(%esp)
|
|
Packit |
5c3484 |
mov (%esi), %eax
|
|
Packit |
5c3484 |
and $3, %ebx
|
|
Packit |
5c3484 |
jz L(c0)
|
|
Packit |
5c3484 |
cmp $2, %ebx
|
|
Packit |
5c3484 |
mov 36(%esp), %ebx
|
|
Packit |
5c3484 |
jz L(c2)
|
|
Packit |
5c3484 |
jg L(c3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(c1): lea -4(%edi), %edi
|
|
Packit |
5c3484 |
mul %ecx
|
|
Packit |
5c3484 |
test %ebp, %ebp
|
|
Packit |
5c3484 |
jnz 1f
|
|
Packit |
5c3484 |
add %ebx, %eax
|
|
Packit |
5c3484 |
mov %eax, 4(%edi)
|
|
Packit |
5c3484 |
mov %edx, %eax
|
|
Packit |
5c3484 |
adc %ebp, %eax
|
|
Packit |
5c3484 |
jmp L(rt)
|
|
Packit |
5c3484 |
1: add %eax, %ebx
|
|
Packit |
5c3484 |
mov $0, %ebp
|
|
Packit |
5c3484 |
adc %edx, %ebp
|
|
Packit |
5c3484 |
mov 4(%esi), %eax
|
|
Packit |
5c3484 |
jmp L(1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(c2): lea 4(%esi), %esi
|
|
Packit |
5c3484 |
mul %ecx
|
|
Packit |
5c3484 |
test %ebp, %ebp
|
|
Packit |
5c3484 |
mov %ebx, %ebp
|
|
Packit |
5c3484 |
jnz 2f
|
|
Packit |
5c3484 |
add %eax, %ebp
|
|
Packit |
5c3484 |
mov $0, %ebx
|
|
Packit |
5c3484 |
adc %edx, %ebx
|
|
Packit |
5c3484 |
mov (%esi), %eax
|
|
Packit |
5c3484 |
jmp L(cj2)
|
|
Packit |
5c3484 |
2: add %eax, %ebp
|
|
Packit |
5c3484 |
mov $0, %ebx
|
|
Packit |
5c3484 |
adc %edx, %ebx
|
|
Packit |
5c3484 |
mov (%esi), %eax
|
|
Packit |
5c3484 |
jmp L(2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(c3): lea 8(%esi), %esi
|
|
Packit |
5c3484 |
lea -12(%edi), %edi
|
|
Packit |
5c3484 |
mul %ecx
|
|
Packit |
5c3484 |
add %eax, %ebx
|
|
Packit |
5c3484 |
mov $0, %ebp
|
|
Packit |
5c3484 |
adc %edx, %ebp
|
|
Packit |
5c3484 |
mov -4(%esi), %eax
|
|
Packit |
5c3484 |
incl 28(%esp)
|
|
Packit |
5c3484 |
jmp L(3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(c0): mov 36(%esp), %ebx
|
|
Packit |
5c3484 |
lea -4(%esi), %esi
|
|
Packit |
5c3484 |
lea -8(%edi), %edi
|
|
Packit |
5c3484 |
mul %ecx
|
|
Packit |
5c3484 |
mov %ebx, %ebp
|
|
Packit |
5c3484 |
add %eax, %ebp
|
|
Packit |
5c3484 |
mov $0, %ebx
|
|
Packit |
5c3484 |
adc %edx, %ebx
|
|
Packit |
5c3484 |
mov 8(%esi), %eax
|
|
Packit |
5c3484 |
jmp L(0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
EPILOGUE()
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
PROLOGUE(mpn_mul_1)
|
|
Packit |
5c3484 |
add $-16, %esp
|
|
Packit |
5c3484 |
mov %ebp, (%esp)
|
|
Packit |
5c3484 |
mov %ebx, 4(%esp)
|
|
Packit |
5c3484 |
mov %esi, 8(%esp)
|
|
Packit |
5c3484 |
mov %edi, 12(%esp)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mov 20(%esp), %edi
|
|
Packit |
5c3484 |
mov 24(%esp), %esi
|
|
Packit |
5c3484 |
mov 28(%esp), %ebp
|
|
Packit |
5c3484 |
mov 32(%esp), %ecx
|
|
Packit |
5c3484 |
mov %ebp, %ebx
|
|
Packit |
5c3484 |
shr $2, %ebp
|
|
Packit |
5c3484 |
mov %ebp, 28(%esp)
|
|
Packit |
5c3484 |
mov (%esi), %eax
|
|
Packit |
5c3484 |
and $3, %ebx
|
|
Packit |
5c3484 |
jz L(b0)
|
|
Packit |
5c3484 |
cmp $2, %ebx
|
|
Packit |
5c3484 |
jz L(b2)
|
|
Packit |
5c3484 |
jg L(b3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b1): lea -4(%edi), %edi
|
|
Packit |
5c3484 |
mul %ecx
|
|
Packit |
5c3484 |
test %ebp, %ebp
|
|
Packit |
5c3484 |
jnz L(gt1)
|
|
Packit |
5c3484 |
mov %eax, 4(%edi)
|
|
Packit |
5c3484 |
mov %edx, %eax
|
|
Packit |
5c3484 |
jmp L(rt)
|
|
Packit |
5c3484 |
L(gt1): mov %eax, %ebx
|
|
Packit |
5c3484 |
mov %edx, %ebp
|
|
Packit |
5c3484 |
mov 4(%esi), %eax
|
|
Packit |
5c3484 |
jmp L(1)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b2): lea 4(%esi), %esi
|
|
Packit |
5c3484 |
mul %ecx
|
|
Packit |
5c3484 |
test %ebp, %ebp
|
|
Packit |
5c3484 |
mov %eax, %ebp
|
|
Packit |
5c3484 |
mov %edx, %ebx
|
|
Packit |
5c3484 |
mov (%esi), %eax
|
|
Packit |
5c3484 |
jnz L(2)
|
|
Packit |
5c3484 |
jmp L(cj2)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b3): lea 8(%esi), %esi
|
|
Packit |
5c3484 |
lea -12(%edi), %edi
|
|
Packit |
5c3484 |
mul %ecx
|
|
Packit |
5c3484 |
mov %eax, %ebx
|
|
Packit |
5c3484 |
mov %edx, %ebp
|
|
Packit |
5c3484 |
mov -4(%esi), %eax
|
|
Packit |
5c3484 |
incl 28(%esp)
|
|
Packit |
5c3484 |
jmp L(3)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(b0): lea -4(%esi), %esi
|
|
Packit |
5c3484 |
lea -8(%edi), %edi
|
|
Packit |
5c3484 |
mul %ecx
|
|
Packit |
5c3484 |
mov %eax, %ebp
|
|
Packit |
5c3484 |
mov %edx, %ebx
|
|
Packit |
5c3484 |
mov 8(%esi), %eax
|
|
Packit |
5c3484 |
jmp L(0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ALIGN(16)
|
|
Packit |
5c3484 |
L(top): mov $0, %ebx
|
|
Packit |
5c3484 |
adc %edx, %ebx
|
|
Packit |
5c3484 |
L(2): mul %ecx
|
|
Packit |
5c3484 |
add %eax, %ebx
|
|
Packit |
5c3484 |
mov %ebp, 0(%edi)
|
|
Packit |
5c3484 |
mov 4(%esi), %eax
|
|
Packit |
5c3484 |
mov $0, %ebp
|
|
Packit |
5c3484 |
adc %edx, %ebp
|
|
Packit |
5c3484 |
L(1): mul %ecx
|
|
Packit |
5c3484 |
add %eax, %ebp
|
|
Packit |
5c3484 |
mov 8(%esi), %eax
|
|
Packit |
5c3484 |
mov %ebx, 4(%edi)
|
|
Packit |
5c3484 |
mov $0, %ebx
|
|
Packit |
5c3484 |
adc %edx, %ebx
|
|
Packit |
5c3484 |
L(0): mov %ebp, 8(%edi)
|
|
Packit |
5c3484 |
mul %ecx
|
|
Packit |
5c3484 |
add %eax, %ebx
|
|
Packit |
5c3484 |
mov 12(%esi), %eax
|
|
Packit |
5c3484 |
lea 16(%esi), %esi
|
|
Packit |
5c3484 |
mov $0, %ebp
|
|
Packit |
5c3484 |
adc %edx, %ebp
|
|
Packit |
5c3484 |
L(3): mov %ebx, 12(%edi)
|
|
Packit |
5c3484 |
mul %ecx
|
|
Packit |
5c3484 |
lea 16(%edi), %edi
|
|
Packit |
5c3484 |
add %eax, %ebp
|
|
Packit |
5c3484 |
decl 28(%esp)
|
|
Packit |
5c3484 |
mov 0(%esi), %eax
|
|
Packit |
5c3484 |
jnz L(top)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(end): mov $0, %ebx
|
|
Packit |
5c3484 |
adc %edx, %ebx
|
|
Packit |
5c3484 |
L(cj2): mul %ecx
|
|
Packit |
5c3484 |
add %eax, %ebx
|
|
Packit |
5c3484 |
mov %ebp, (%edi)
|
|
Packit |
5c3484 |
L(cj1): mov %ebx, 4(%edi)
|
|
Packit |
5c3484 |
adc $0, %edx
|
|
Packit |
5c3484 |
mov %edx, %eax
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
L(rt): mov (%esp), %ebp
|
|
Packit |
5c3484 |
mov 4(%esp), %ebx
|
|
Packit |
5c3484 |
mov 8(%esp), %esi
|
|
Packit |
5c3484 |
mov 12(%esp), %edi
|
|
Packit |
5c3484 |
add $16, %esp
|
|
Packit |
5c3484 |
ret
|
|
Packit |
5c3484 |
EPILOGUE()
|
|
Packit |
5c3484 |
ASM_END()
|