Index: gcc/gmp/mpn/alpha/diveby3.asm |
diff --git a/gcc/gmp/mpn/alpha/diveby3.asm b/gcc/gmp/mpn/alpha/diveby3.asm |
deleted file mode 100644 |
index e2d1c6beee9e9c6f2af8fb891133616f99342f89..0000000000000000000000000000000000000000 |
--- a/gcc/gmp/mpn/alpha/diveby3.asm |
+++ /dev/null |
@@ -1,322 +0,0 @@ |
-dnl Alpha mpn_divexact_by3c -- mpn division by 3, expecting no remainder. |
- |
-dnl Copyright 2004, 2005 Free Software Foundation, Inc. |
- |
-dnl This file is part of the GNU MP Library. |
- |
-dnl The GNU MP Library is free software; you can redistribute it and/or modify |
-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 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: 22 |
-C EV5: 11.5 |
-C EV6: 6.3 |
- |
-C TODO |
-C * Trim this to 6.0 c/l for ev6. |
-C * Write special ev5 version, should reach 9 c/l, and could be smaller. |
-C * Try prefetch for destination, using lds. |
-C * Improve feed-in code, by moving initial mulq earlier; make initial load |
-C to u0/u0 to save some copying. |
-C * Combine u0 and u2, u1 and u3. |
- |
-C INPUT PARAMETERS |
-define(`rp', `r16') |
-define(`up', `r17') |
-define(`n', `r18') |
-define(`cy', `r19') |
- |
-ASM_START() |
- |
-DATASTART(L(LC)) |
- .quad 0xAAAAAAAAAAAAAAAB |
- .quad 0x5555555555555555 |
- .quad 0xAAAAAAAAAAAAAAAA |
-DATAEND() |
- |
-define(`xAAAAAAAAAAAAAAAB', `r20') |
-define(`x5555555555555555', `r21') |
-define(`xAAAAAAAAAAAAAAAA', `r22') |
-define(`u0', `r0') define(`u1', `r1') |
-define(`u2', `r2') define(`u3', `r3') |
-define(`l0', `r25') define(`x', `r8') |
-define(`q0', `r4') define(`q1', `r5') |
-define(`p6', `r6') define(`p7', `r7') |
-define(`t0', `r23') define(`t1', `r24') |
-define(`cymask',`r28') |
- |
- |
-PROLOGUE(mpn_divexact_by3c,gp) |
- |
- ldq r28, 0(up) C load first limb early |
- |
-C Put magic constants in registers |
- lda r0, L(LC) |
- ldq xAAAAAAAAAAAAAAAB, 0(r0) |
- ldq x5555555555555555, 8(r0) |
- ldq xAAAAAAAAAAAAAAAA, 16(r0) |
- |
-C Compute initial l0 value |
- cmpeq cy, 1, p6 |
- cmpeq cy, 2, p7 |
- negq p6, p6 |
- and p6, x5555555555555555, l0 |
- cmovne p7, xAAAAAAAAAAAAAAAA, l0 |
- |
-C Feed-in depending on (n mod 4) |
- and n, 3, r8 |
- lda n, -3(n) |
- cmpeq r8, 1, r4 |
- cmpeq r8, 2, r5 |
- bne r4, $Lb01 |
- bne r5, $Lb10 |
- beq r8, $Lb00 |
- |
-$Lb11: ldq u3, 8(up) |
- lda up, -24(up) |
- lda rp, -24(rp) |
- mulq r28, xAAAAAAAAAAAAAAAB, q0 |
- mov r28, u2 |
- br r31, $L11 |
- |
-$Lb00: ldq u2, 8(up) |
- lda up, -16(up) |
- lda rp, -16(rp) |
- mulq r28, xAAAAAAAAAAAAAAAB, q1 |
- mov r28, u1 |
- br r31, $L00 |
- |
-$Lb01: lda rp, -8(rp) |
- mulq r28, xAAAAAAAAAAAAAAAB, q0 |
- mov r28, u0 |
- blt n, $Lcj1 |
- ldq u1, 8(up) |
- lda up, -8(up) |
- br r31, $L01 |
- |
-$Lb10: ldq u0, 8(up) |
- mulq r28, xAAAAAAAAAAAAAAAB, q1 |
- mov r28, u3 |
- blt n, $Lend |
- |
- ALIGN(16) |
-$Ltop: |
-C 0 |
- cmpult u3, cy, cy C L0 |
- mulq u0, xAAAAAAAAAAAAAAAB, q0 C U1 |
- ldq u1, 16(up) C L1 |
- addq q1, l0, x C U0 |
-C 1 |
- negq cy, cymask C L0 |
- unop C U1 |
- unop C L1 |
- cmpult x5555555555555555, x, p6 C U0 |
-C 2 |
- cmpult xAAAAAAAAAAAAAAAA, x, p7 C U1 |
- unop |
- unop |
- negq p6, t0 C L0 |
-C 3 |
- negq p7, t1 C L0 |
- and cymask, x5555555555555555, l0 C U1 |
- addq p6, cy, cy |
- and t0, x5555555555555555, t0 |
-C 4 |
- and t1, x5555555555555555, t1 |
- addq p7, cy, cy |
- unop |
- addq t0, l0, l0 |
-C 5 |
- addq t1, l0, l0 |
- unop |
- stq x, 0(rp) C L1 |
- unop |
-$L01: |
-C 0 |
- cmpult u0, cy, cy C L0 |
- mulq u1, xAAAAAAAAAAAAAAAB, q1 C U1 |
- ldq u2, 24(up) C L1 |
- addq q0, l0, x C U0 |
-C 1 |
- negq cy, cymask C L0 |
- unop C U1 |
- unop C L1 |
- cmpult x5555555555555555, x, p6 C U0 |
-C 2 |
- cmpult xAAAAAAAAAAAAAAAA, x, p7 C U1 |
- unop |
- unop |
- negq p6, t0 C L0 |
-C 3 |
- negq p7, t1 C L0 |
- and cymask, x5555555555555555, l0 C U1 |
- addq p6, cy, cy |
- and t0, x5555555555555555, t0 |
-C 4 |
- and t1, x5555555555555555, t1 |
- addq p7, cy, cy |
- unop |
- addq t0, l0, l0 |
-C 5 |
- addq t1, l0, l0 |
- unop |
- stq x, 8(rp) C L1 |
- unop |
-$L00: |
-C 0 |
- cmpult u1, cy, cy C L0 |
- mulq u2, xAAAAAAAAAAAAAAAB, q0 C U1 |
- ldq u3, 32(up) C L1 |
- addq q1, l0, x C U0 |
-C 1 |
- negq cy, cymask C L0 |
- unop C U1 |
- unop C L1 |
- cmpult x5555555555555555, x, p6 C U0 |
-C 2 |
- cmpult xAAAAAAAAAAAAAAAA, x, p7 C U1 |
- unop |
- unop |
- negq p6, t0 C L0 |
-C 3 |
- negq p7, t1 C L0 |
- and cymask, x5555555555555555, l0 C U1 |
- addq p6, cy, cy |
- and t0, x5555555555555555, t0 |
-C 4 |
- and t1, x5555555555555555, t1 |
- addq p7, cy, cy |
- unop |
- addq t0, l0, l0 |
-C 5 |
- addq t1, l0, l0 |
- unop |
- stq x, 16(rp) C L1 |
- unop |
-$L11: |
-C 0 |
- cmpult u2, cy, cy C L0 |
- mulq u3, xAAAAAAAAAAAAAAAB, q1 C U1 |
- ldq u0, 40(up) C L1 |
- addq q0, l0, x C U0 |
-C 1 |
- negq cy, cymask C L0 |
- unop C U1 |
- unop C L1 |
- cmpult x5555555555555555, x, p6 C U0 |
-C 2 |
- cmpult xAAAAAAAAAAAAAAAA, x, p7 C U1 |
- lda n, -4(n) C L1 bookkeeping |
- unop |
- negq p6, t0 C L0 |
-C 3 |
- negq p7, t1 C L0 |
- and cymask, x5555555555555555, l0 C U1 |
- addq p6, cy, cy |
- and t0, x5555555555555555, t0 |
-C 4 |
- and t1, x5555555555555555, t1 |
- addq p7, cy, cy |
- unop |
- addq t0, l0, l0 |
-C 5 |
- addq t1, l0, l0 |
- unop |
- stq x, 24(rp) C L1 |
- lda up, 32(up) |
-C |
- ldl r31, 256(up) C prefetch |
- unop |
- lda rp, 32(rp) |
- bge n, $Ltop C U1 |
-C *** MAIN LOOP END *** |
-$Lend: |
- |
- cmpult u3, cy, cy C L0 |
- mulq u0, xAAAAAAAAAAAAAAAB, q0 C U1 |
- unop |
- addq q1, l0, x C U0 |
-C 1 |
- negq cy, cymask C L0 |
- unop C U1 |
- unop C L1 |
- cmpult x5555555555555555, x, p6 C U0 |
-C 2 |
- cmpult xAAAAAAAAAAAAAAAA, x, p7 C U1 |
- unop |
- unop |
- negq p6, t0 C L0 |
-C 3 |
- negq p7, t1 C L0 |
- and cymask, x5555555555555555, l0 C U1 |
- addq p6, cy, cy |
- and t0, x5555555555555555, t0 |
-C 4 |
- and t1, x5555555555555555, t1 |
- addq p7, cy, cy |
- unop |
- addq t0, l0, l0 |
-C 5 |
- addq t1, l0, l0 |
- unop |
- stq x, 0(rp) C L1 |
- unop |
-$Lcj1: |
- cmpult u0, cy, cy C L0 |
- addq q0, l0, x C U0 |
- cmpult x5555555555555555, x, p6 C U0 |
- cmpult xAAAAAAAAAAAAAAAA, x, p7 C U1 |
- addq p6, cy, cy |
- addq p7, cy, r0 |
- stq x, 8(rp) C L1 |
- |
- ret r31,(r26),1 |
-EPILOGUE() |
-ASM_END() |
- |
-C This is useful for playing with various schedules. |
-C Expand as: one(0)one(1)one(2)one(3) |
-define(`one',` |
-C 0 |
- cmpult `$'eval(($1+3)%4), cy, cy C L0 |
- mulq `$'$1, xAAAAAAAAAAAAAAAB, `$'eval(4+$1%2) C U1 |
- ldq `$'eval(($1+1)%4), eval($1*8+16)(up) C L1 |
- addq `$'eval(4+($1+1)%2), l0, x C U0 |
-C 1 |
- negq cy, cymask C L0 |
- unop C U1 |
- unop C L1 |
- cmpult x5555555555555555, x, p6 C U0 |
-C 2 |
- cmpult xAAAAAAAAAAAAAAAA, x, p7 C U1 |
- unop |
- unop |
- negq p6, t0 C L0 |
-C 3 |
- negq p7, t1 C L0 |
- and cymask, x5555555555555555, l0 C U1 |
- addq p6, cy, cy |
- and t0, x5555555555555555, t0 |
-C 4 |
- and t1, x5555555555555555, t1 |
- addq p7, cy, cy |
- unop |
- addq t0, l0, l0 |
-C 5 |
- addq t1, l0, l0 |
- unop |
- stq x, eval($1*8)(rp) C L1 |
- unop |
-') |