mpir/mpn/x86w/p4/sse2/divexact_1.asm

108 lines
2.4 KiB
NASM
Raw Normal View History

2012-11-25 17:13:44 -05:00
; Copyright 2001, 2002 Free Software Foundation, Inc.
;
; This file is part of the GNU MP Library.
;
; The GNU MP Library is free software; you can redistribute it and/or
; modify it under the terms of the GNU Lesser General Public License as
; published by the Free Software Foundation; either version 2.1 of the
; License, or (at your option) any later version.
;
; The GNU MP Library is distributed in the hope that it will be useful,
; but WITHOUT ANY WARRANTY; without even the implied warranty of
; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
; Lesser General Public License for more details.
;
; You should have received a copy of the GNU Lesser General Public
; License along with the GNU MP Library; see the file COPYING.LIB. If
; not, write to the Free Software Foundation, Inc., 59 Temple Place -
; Suite 330, Boston, MA 02111-1307, USA.
;
; Translation of AT&T syntax code by Brian Gladman
%include "..\..\x86i.inc"
%define p_space 16
extern ___gmp_modlimb_invert_table
section .text
global ___gmpn_divexact_1
%ifdef DLL
export ___gmpn_divexact_1
%endif
align 16
___gmpn_divexact_1:
mov edx,[12+esp]
mov eax,[8+esp]
mov ecx,[16+esp]
sub edx,1
jnz two_or_more
mov eax,[eax]
xor edx,edx
div ecx
mov ecx,[4+esp]
mov [ecx],eax
ret
two_or_more:
mov eax,ecx
bsf ecx,ecx
shr eax,cl
movd mm6,eax
movd mm7,ecx
shr eax,1
and eax,127
movzx eax,byte [___gmp_modlimb_invert_table+eax]
movd mm5,eax
movd mm0,eax
pmuludq mm5,mm5
pmuludq mm5,mm6
paddd mm0,mm0
psubd mm0,mm5
pxor mm5,mm5
paddd mm5,mm0
pmuludq mm0,mm0
pcmpeqd mm4,mm4
psrlq mm4,32
pmuludq mm0,mm6
paddd mm5,mm5
mov eax,[8+esp]
mov ecx,[4+esp]
pxor mm1,mm1
psubd mm5,mm0
pxor mm0,mm0
top:
movd mm2,[eax]
movd mm3,[4+eax]
add eax,4
punpckldq mm2,mm3
psrlq mm2,mm7
pand mm2,mm4
psubq mm2,mm0
psubq mm2,mm1
movq mm0,mm2
psrlq mm0,63
pmuludq mm2,mm5
movd [ecx],mm2
add ecx,4
movq mm1,mm6
pmuludq mm1,mm2
psrlq mm1,32
sub edx,1
jnz top
done:
movd mm2,[eax]
psrlq mm2,mm7
psubq mm2,mm0
psubq mm2,mm1
pmuludq mm2,mm5
movd [ecx],mm2
emms
ret
end