diff options
Diffstat (limited to 'gmp/mpn/alpha/lshift.asm')
| -rw-r--r-- | gmp/mpn/alpha/lshift.asm | 193 |
1 files changed, 54 insertions, 139 deletions
diff --git a/gmp/mpn/alpha/lshift.asm b/gmp/mpn/alpha/lshift.asm index c62a856aea..eb5b2a0b68 100644 --- a/gmp/mpn/alpha/lshift.asm +++ b/gmp/mpn/alpha/lshift.asm @@ -1,39 +1,28 @@ dnl Alpha mpn_lshift -- Shift a number left. -dnl Copyright 1994, 1995, 2000, 2003, 2009 Free Software Foundation, Inc. +dnl Copyright 1994, 1995, 2000, 2002, 2003 Free Software Foundation, Inc. dnl This file is part of the GNU MP Library. -dnl + dnl The GNU MP Library is free software; you can redistribute it and/or modify -dnl it under the terms of either: -dnl -dnl * the GNU Lesser General Public License as published by the Free -dnl Software Foundation; either version 3 of the License, or (at your -dnl option) any later version. -dnl -dnl or -dnl -dnl * the GNU General Public License as published by the Free Software -dnl Foundation; either version 2 of the License, or (at your option) any -dnl later version. -dnl -dnl or both in parallel, as here. -dnl +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 General Public License -dnl for more details. -dnl -dnl You should have received copies of the GNU General Public License and the -dnl GNU Lesser General Public License along with the GNU MP Library. If not, -dnl see https://www.gnu.org/licenses/. +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 cycles/limb -C EV4: ? -C EV5: 3.25 -C EV6: 1.75 +C EV4: 4.75 +C EV5: 4 +C EV6: 2 C INPUT PARAMETERS C rp r16 @@ -46,137 +35,63 @@ ASM_START() PROLOGUE(mpn_lshift) s8addq r18,r17,r17 C make r17 point at end of s1 ldq r4,-8(r17) C load first limb - subq r31,r19,r20 + subq r17,8,r17 + subq r31,r19,r7 s8addq r18,r16,r16 C make r16 point at end of RES subq r18,1,r18 - and r18,4-1,r28 C number of limbs in first loop - srl r4,r20,r0 C compute function result + and r18,4-1,r20 C number of limbs in first loop + srl r4,r7,r0 C compute function result - beq r28,L(L0) - subq r18,r28,r18 + beq r20,$L0 + subq r18,r20,r18 ALIGN(8) -L(top0): - ldq r3,-16(r17) +$Loop0: ldq r3,-8(r17) subq r16,8,r16 - sll r4,r19,r5 subq r17,8,r17 - subq r28,1,r28 - srl r3,r20,r6 + subq r20,1,r20 + sll r4,r19,r5 + srl r3,r7,r6 bis r3,r3,r4 bis r5,r6,r8 stq r8,0(r16) - bne r28,L(top0) + bne r20,$Loop0 -L(L0): sll r4,r19,r24 - beq r18,L(end) -C warm up phase 1 - ldq r1,-16(r17) - subq r18,4,r18 - ldq r2,-24(r17) - ldq r3,-32(r17) - ldq r4,-40(r17) -C warm up phase 2 - srl r1,r20,r7 - sll r1,r19,r21 - srl r2,r20,r8 - beq r18,L(end1) - ldq r1,-48(r17) - sll r2,r19,r22 - ldq r2,-56(r17) - srl r3,r20,r5 - bis r7,r24,r7 - sll r3,r19,r23 - bis r8,r21,r8 - srl r4,r20,r6 - ldq r3,-64(r17) - sll r4,r19,r24 - ldq r4,-72(r17) - subq r18,4,r18 - beq r18,L(end2) - ALIGN(16) -C main loop -L(top): stq r7,-8(r16) - bis r5,r22,r5 - stq r8,-16(r16) - bis r6,r23,r6 - - srl r1,r20,r7 - subq r18,4,r18 - sll r1,r19,r21 - unop C ldq r31,-96(r17) - - srl r2,r20,r8 - ldq r1,-80(r17) - sll r2,r19,r22 - ldq r2,-88(r17) - - stq r5,-24(r16) - bis r7,r24,r7 - stq r6,-32(r16) - bis r8,r21,r8 - - srl r3,r20,r5 - unop C ldq r31,-96(r17) - sll r3,r19,r23 +$L0: beq r18,$Lend + + ALIGN(8) +$Loop: ldq r3,-8(r17) subq r16,32,r16 + subq r18,4,r18 + sll r4,r19,r5 + srl r3,r7,r6 + + ldq r4,-16(r17) + sll r3,r19,r1 + bis r5,r6,r8 + stq r8,24(r16) + srl r4,r7,r2 + + ldq r3,-24(r17) + sll r4,r19,r5 + bis r1,r2,r8 + stq r8,16(r16) + srl r3,r7,r6 - srl r4,r20,r6 - ldq r3,-96(r17) - sll r4,r19,r24 - ldq r4,-104(r17) + ldq r4,-32(r17) + sll r3,r19,r1 + bis r5,r6,r8 + stq r8,8(r16) + srl r4,r7,r2 subq r17,32,r17 - bne r18,L(top) -C cool down phase 2/1 -L(end2): - stq r7,-8(r16) - bis r5,r22,r5 - stq r8,-16(r16) - bis r6,r23,r6 - srl r1,r20,r7 - sll r1,r19,r21 - srl r2,r20,r8 - sll r2,r19,r22 - stq r5,-24(r16) - bis r7,r24,r7 - stq r6,-32(r16) - bis r8,r21,r8 - srl r3,r20,r5 - sll r3,r19,r23 - srl r4,r20,r6 - sll r4,r19,r24 -C cool down phase 2/2 - stq r7,-40(r16) - bis r5,r22,r5 - stq r8,-48(r16) - bis r6,r23,r6 - stq r5,-56(r16) - stq r6,-64(r16) -C cool down phase 2/3 - stq r24,-72(r16) - ret r31,(r26),1 + bis r1,r2,r8 + stq r8,0(r16) -C cool down phase 1/1 -L(end1): - sll r2,r19,r22 - srl r3,r20,r5 - bis r7,r24,r7 - sll r3,r19,r23 - bis r8,r21,r8 - srl r4,r20,r6 - sll r4,r19,r24 -C cool down phase 1/2 - stq r7,-8(r16) - bis r5,r22,r5 - stq r8,-16(r16) - bis r6,r23,r6 - stq r5,-24(r16) - stq r6,-32(r16) - stq r24,-40(r16) - ret r31,(r26),1 + bgt r18,$Loop -L(end): stq r24,-8(r16) +$Lend: sll r4,r19,r8 + stq r8,-8(r16) ret r31,(r26),1 EPILOGUE(mpn_lshift) ASM_END() |