version 1.1.1.1, 2000/09/09 14:12:42 |
version 1.1.1.2, 2003/08/25 16:06:28 |
|
|
dnl AMD K7 mpn_divexact_by3 -- mpn division by 3, expecting no remainder. |
dnl AMD K7 mpn_divexact_by3 -- mpn division by 3, expecting no remainder. |
dnl |
|
dnl K7: 8.0 cycles/limb |
|
|
|
|
dnl Copyright 2000, 2002 Free Software Foundation, Inc. |
dnl Copyright (C) 2000 Free Software Foundation, Inc. |
|
dnl |
dnl |
dnl This file is part of the GNU MP Library. |
dnl This file is part of the GNU MP Library. |
dnl |
dnl |
Line 22 dnl License along with the GNU MP Library; see the fi |
|
Line 19 dnl License along with the GNU MP Library; see the fi |
|
dnl not, write to the Free Software Foundation, Inc., 59 Temple Place - |
dnl not, write to the Free Software Foundation, Inc., 59 Temple Place - |
dnl Suite 330, Boston, MA 02111-1307, USA. |
dnl Suite 330, Boston, MA 02111-1307, USA. |
|
|
|
|
include(`../config.m4') |
include(`../config.m4') |
|
|
|
|
|
C K7: 8.0 cycles/limb |
|
|
|
|
C mp_limb_t mpn_divexact_by3c (mp_ptr dst, mp_srcptr src, mp_size_t size, |
C mp_limb_t mpn_divexact_by3c (mp_ptr dst, mp_srcptr src, mp_size_t size, |
C mp_limb_t carry); |
C mp_limb_t carry); |
|
|
Line 41 dnl ceil(b/3) and floor(b*2/3) where b=2^32 |
|
Line 40 dnl ceil(b/3) and floor(b*2/3) where b=2^32 |
|
deflit(ONE_THIRD_CEIL, 0x55555556) |
deflit(ONE_THIRD_CEIL, 0x55555556) |
deflit(TWO_THIRDS_FLOOR, 0xAAAAAAAA) |
deflit(TWO_THIRDS_FLOOR, 0xAAAAAAAA) |
|
|
.text |
TEXT |
ALIGN(32) |
ALIGN(32) |
|
|
PROLOGUE(mpn_divexact_by3c) |
PROLOGUE(mpn_divexact_by3c) |
|
|
C edi &dst[size-1] |
C edi &dst[size-1] |
C ebp |
C ebp |
|
|
imull $INVERSE_3, %eax |
imull $INVERSE_3, %eax, %eax |
|
|
cmpl $ONE_THIRD_CEIL, %eax |
cmpl $ONE_THIRD_CEIL, %eax |
movl %eax, (%edi) |
movl %eax, (%edi) |