| version 1.1.1.1, 2000/09/09 14:12:41 |
version 1.1.1.2, 2003/08/25 16:06:26 |
|
|
| ! SPARC v9 __gmpn_rshift -- |
dnl SPARC v9 mpn_rshift |
| |
|
| ! Copyright (C) 1996, 2000 Free Software Foundation, Inc. |
dnl Copyright 1996, 2000, 2001, 2002 Free Software Foundation, Inc. |
| |
|
| ! This file is part of the GNU MP Library. |
dnl This file is part of the GNU MP Library. |
| |
|
| ! The GNU MP Library is free software; you can redistribute it and/or modify |
dnl 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 |
dnl 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 |
dnl the Free Software Foundation; either version 2.1 of the License, or (at your |
| ! option) any later version. |
dnl option) any later version. |
| |
|
| ! The GNU MP Library is distributed in the hope that it will be useful, but |
dnl 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 |
dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY |
| ! or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public |
dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public |
| ! License for more details. |
dnl License for more details. |
| |
|
| ! You should have received a copy of the GNU Lesser General Public License |
dnl 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 |
dnl 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, |
dnl the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, |
| ! MA 02111-1307, USA. |
dnl MA 02111-1307, USA. |
| |
|
| |
|
| ! INPUT PARAMETERS |
|
| ! res_ptr %o0 |
|
| ! src_ptr %o1 |
|
| ! size %o2 |
|
| ! cnt %o3 |
|
| |
|
| include(`../config.m4') |
include(`../config.m4') |
| |
|
| ASM_START() |
C INPUT PARAMETERS |
| .register %g2,#scratch |
define(`rp',`%i0') |
| .register %g3,#scratch |
define(`up',`%i1') |
| PROLOGUE(mpn_rshift) |
define(`n',`%i2') |
| ldx [%o1],%g2 ! load first limb |
define(`cnt',`%i3') |
| sub %g0,%o3,%o5 ! negate shift count |
|
| add %o2,-1,%o2 |
|
| and %o2,4-1,%g4 ! number of limbs in first loop |
|
| sllx %g2,%o5,%g1 ! compute function result |
|
| brz,pn %g4,L(0) ! if multiple of 4 limbs, skip first loop |
|
| mov %g1,%g5 |
|
| |
|
| sub %o2,%g4,%o2 ! adjust count for main loop |
define(`u0',`%l0') |
| |
define(`u1',`%l2') |
| |
define(`u2',`%l4') |
| |
define(`u3',`%l6') |
| |
|
| L(loop0): |
define(`tnc',`%i4') |
| ldx [%o1+8],%g3 |
|
| add %o0,8,%o0 |
|
| add %o1,8,%o1 |
|
| add %g4,-1,%g4 |
|
| srlx %g2,%o3,%o4 |
|
| sllx %g3,%o5,%g1 |
|
| mov %g3,%g2 |
|
| or %o4,%g1,%o4 |
|
| brnz,pt %g4,L(loop0) |
|
| stx %o4,[%o0-8] |
|
| |
|
| L(0): brz,pn %o2,L(end) |
define(`fanop',`fitod %f0,%f2') dnl A quasi nop running in the FA pipe |
| nop |
define(`fmnop',`fmuld %f0,%f0,%f4') dnl A quasi nop running in the FM pipe |
| |
|
| L(loop1): |
ASM_START() |
| ldx [%o1+8],%g3 |
REGISTER(%g2,#scratch) |
| add %o0,32,%o0 |
REGISTER(%g3,#scratch) |
| add %o2,-4,%o2 |
PROLOGUE(mpn_rshift) |
| srlx %g2,%o3,%o4 |
save %sp,-160,%sp |
| sllx %g3,%o5,%g1 |
|
| |
|
| ldx [%o1+16],%g2 |
sub %g0,cnt,tnc C negate shift count |
| srlx %g3,%o3,%g4 |
ldx [up],u3 C load first limb |
| or %o4,%g1,%o4 |
subcc n,5,n |
| stx %o4,[%o0-32] |
sllx u3,tnc,%i5 C compute function result |
| sllx %g2,%o5,%g1 |
srlx u3,cnt,%g3 |
| |
bl,pn %icc,.Lend1234 |
| |
fanop |
| |
|
| ldx [%o1+24],%g3 |
subcc n,4,n |
| srlx %g2,%o3,%o4 |
ldx [up+8],u0 |
| or %g4,%g1,%g4 |
ldx [up+16],u1 |
| stx %g4,[%o0-24] |
add up,32,up |
| sllx %g3,%o5,%g1 |
ldx [up-8],u2 |
| |
ldx [up+0],u3 |
| |
sllx u0,tnc,%g2 |
| |
|
| ldx [%o1+32],%g2 |
bl,pn %icc,.Lend5678 |
| srlx %g3,%o3,%g4 |
fanop |
| or %o4,%g1,%o4 |
|
| stx %o4,[%o0-16] |
|
| sllx %g2,%o5,%g1 |
|
| |
|
| add %o1,32,%o1 |
b,a .Loop |
| or %g4,%g1,%g4 |
.align 16 |
| brnz %o2,L(loop1) |
.Loop: |
| stx %g4,[%o0-8] |
srlx u0,cnt,%g1 |
| |
or %g3,%g2,%g3 |
| |
ldx [up+8],u0 |
| |
fanop |
| |
C -- |
| |
sllx u1,tnc,%g2 |
| |
subcc n,4,n |
| |
stx %g3,[rp+0] |
| |
fanop |
| |
C -- |
| |
srlx u1,cnt,%g3 |
| |
or %g1,%g2,%g1 |
| |
ldx [up+16],u1 |
| |
fanop |
| |
C -- |
| |
sllx u2,tnc,%g2 |
| |
stx %g1,[rp+8] |
| |
add up,32,up |
| |
fanop |
| |
C -- |
| |
srlx u2,cnt,%g1 |
| |
or %g3,%g2,%g3 |
| |
ldx [up-8],u2 |
| |
fanop |
| |
C -- |
| |
sllx u3,tnc,%g2 |
| |
stx %g3,[rp+16] |
| |
add rp,32,rp |
| |
fanop |
| |
C -- |
| |
srlx u3,cnt,%g3 |
| |
or %g1,%g2,%g1 |
| |
ldx [up+0],u3 |
| |
fanop |
| |
C -- |
| |
sllx u0,tnc,%g2 |
| |
stx %g1,[rp-8] |
| |
bge,pt %icc,.Loop |
| |
fanop |
| |
C -- |
| |
.Lend5678: |
| |
srlx u0,cnt,%g1 |
| |
or %g3,%g2,%g3 |
| |
sllx u1,tnc,%g2 |
| |
stx %g3,[rp+0] |
| |
srlx u1,cnt,%g3 |
| |
or %g1,%g2,%g1 |
| |
sllx u2,tnc,%g2 |
| |
stx %g1,[rp+8] |
| |
srlx u2,cnt,%g1 |
| |
or %g3,%g2,%g3 |
| |
sllx u3,tnc,%g2 |
| |
stx %g3,[rp+16] |
| |
add rp,32,rp |
| |
srlx u3,cnt,%g3 C carry... |
| |
or %g1,%g2,%g1 |
| |
stx %g1,[rp-8] |
| |
|
| L(end): srlx %g2,%o3,%g2 |
.Lend1234: |
| stx %g2,[%o0-0] |
addcc n,4,n |
| retl |
bz,pn %icc,.Lret |
| mov %g5,%o0 |
fanop |
| |
.Loop0: |
| |
add rp,8,rp |
| |
subcc n,1,n |
| |
ldx [up+8],u3 |
| |
add up,8,up |
| |
sllx u3,tnc,%g2 |
| |
or %g3,%g2,%g3 |
| |
stx %g3,[rp-8] |
| |
srlx u3,cnt,%g3 |
| |
bnz,pt %icc,.Loop0 |
| |
fanop |
| |
.Lret: |
| |
stx %g3,[rp+0] |
| |
mov %i5,%i0 |
| |
ret |
| |
restore |
| EPILOGUE(mpn_rshift) |
EPILOGUE(mpn_rshift) |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to rshift.asm CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / mpn / sparc64