Index: gcc/gmp/mpn/ia64/divrem_1.asm |
diff --git a/gcc/gmp/mpn/ia64/divrem_1.asm b/gcc/gmp/mpn/ia64/divrem_1.asm |
deleted file mode 100644 |
index aa50ac902bceb43dbcce898209b4d1b4abe43db6..0000000000000000000000000000000000000000 |
--- a/gcc/gmp/mpn/ia64/divrem_1.asm |
+++ /dev/null |
@@ -1,464 +0,0 @@ |
-dnl IA-64 mpn_divrem_1 and mpn_preinv_divrem_1 -- Divide an mpn number by an |
-dnl unnormalized limb. |
- |
-dnl Copyright 2002, 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 Itanium: 40-42 |
-C Itanium 2: 29-30 |
- |
-C This was generated by gcc, then the loops were optimized. The preinv entry |
-C point was shoehorned into the file. Lots of things outside the loops could |
-C be streamlined. It would probably be a good idea to merge the loops for |
-C normalized and unnormalized divisor, since the shifting stuff is done for |
-C free in parallel with other operations. It would even be possible to merge |
-C all loops, if the ld8 were made conditional. |
- |
-C TODO |
-C * Consider delaying inversion for normalized mpn_divrem_1 entry till after |
-C computing leading limb. |
-C * Inline and interleave limb inversion code with loop setup code. |
- |
-ASM_START() |
- |
-C HP's assembler requires these declarations for importing mpn_invert_limb |
- .global mpn_invert_limb |
- .type mpn_invert_limb,@function |
- |
-C INPUT PARAMETERS |
-C rp = r32 |
-C qxn = r33 |
-C up = r34 |
-C n = r35 |
-C vl = r36 |
-C vlinv = r37 (preinv only) |
-C cnt = r38 (preinv only) |
- |
-PROLOGUE(mpn_preinv_divrem_1) |
- .prologue |
- .save ar.pfs, r42 |
- alloc r42 = ar.pfs, 7, 8, 1, 0 |
- .save ar.lc, r44 |
- mov r44 = ar.lc |
- .save rp, r41 |
- mov r41 = b0 |
- .body |
-ifdef(`HAVE_ABI_32', |
-` addp4 r32 = 0, r32 |
- sxt4 r33 = r33 |
- addp4 r34 = 0, r34 |
- sxt4 r35 = r35 |
- ;; |
-') |
- mov r40 = r38 |
- shladd r34 = r35, 3, r34 |
- ;; |
- adds r34 = -8, r34 |
- ;; |
- ld8 r39 = [r34], -8 |
- ;; |
- |
- add r15 = r35, r33 |
- ;; |
- mov r8 = r37 |
- shladd r32 = r15, 3, r32 C r32 = rp + n + qxn |
- cmp.le p8, p0 = 0, r36 |
- ;; |
- adds r32 = -8, r32 C r32 = rp + n + qxn - 1 |
- cmp.leu p6, p7 = r36, r39 |
- (p8) br.cond.dpnt .Lpunnorm |
- ;; |
- |
- (p6) addl r15 = 1, r0 |
- (p7) mov r15 = r0 |
- ;; |
- (p6) sub r38 = r39, r36 |
- (p7) mov r38 = r39 |
- st8 [r32] = r15, -8 |
- adds r35 = -2, r35 C un -= 2 |
- br .Lpn |
- |
-.Lpunnorm: |
- (p6) add r34 = 8, r34 |
- mov r38 = 0 C r = 0 |
- shl r36 = r36, r40 |
- (p6) br.cond.dptk .Lpu |
- ;; |
- shl r38 = r39, r40 C r = ahigh << cnt |
- cmp.ne p8, p0 = 1, r35 |
- st8 [r32] = r0, -8 |
- adds r35 = -1, r35 C un-- |
- (p8) br.cond.dpnt .Lpu |
- |
- mov r23 = 1 |
- ;; |
- setf.sig f6 = r8 |
- setf.sig f12 = r23 |
- br .L435 |
-EPILOGUE() |
- |
- |
-PROLOGUE(mpn_divrem_1) |
- .prologue |
- .save ar.pfs, r42 |
- alloc r42 = ar.pfs, 5, 8, 1, 0 |
- .save ar.lc, r44 |
- mov r44 = ar.lc |
- .save rp, r41 |
- mov r41 = b0 |
- .body |
-ifdef(`HAVE_ABI_32', |
-` addp4 r32 = 0, r32 |
- sxt4 r33 = r33 |
- addp4 r34 = 0, r34 |
- sxt4 r35 = r35 |
- ;; |
-') |
- mov r38 = r0 |
- add r15 = r35, r33 |
- ;; |
- cmp.ne p6, p7 = 0, r15 |
- ;; |
- (p7) mov r8 = r0 |
- (p7) br.cond.dpnt .Lret |
- shladd r14 = r15, 3, r32 C r14 = rp + n + qxn |
- cmp.le p6, p7 = 0, r36 |
- ;; |
- adds r32 = -8, r14 C r32 = rp + n + qxn - 1 |
- (p6) br.cond.dpnt .Lunnorm |
- cmp.eq p6, p7 = 0, r35 |
- (p6) br.cond.dpnt .L179 |
- shladd r14 = r35, 3, r34 |
- ;; |
- adds r14 = -8, r14 |
- adds r35 = -1, r35 |
- ;; |
- ld8 r38 = [r14] |
- ;; |
- cmp.leu p6, p7 = r36, r38 |
- ;; |
- (p6) addl r15 = 1, r0 |
- (p7) mov r15 = r0 |
- ;; |
- st8 [r32] = r15, -8 |
- (p6) sub r38 = r38, r36 |
- |
-.L179: |
- mov r45 = r36 |
- adds r35 = -1, r35 |
- br.call.sptk.many b0 = mpn_invert_limb |
- ;; |
- shladd r34 = r35, 3, r34 |
-.Lpn: |
- mov r23 = 1 |
- ;; |
- setf.sig f6 = r8 |
- setf.sig f12 = r23 |
- cmp.le p6, p7 = 0, r35 |
- mov r40 = 0 |
- (p7) br.cond.dpnt .L435 |
- setf.sig f10 = r36 |
- mov ar.lc = r35 |
- setf.sig f7 = r38 |
- ;; |
- sub r28 = -1, r36 |
-C Develop quotient limbs for normalized divisor |
-.Loop1: C 00 C q=r18 nh=r38/f7 |
- ld8 r20 = [r34], -8 |
- xma.hu f11 = f7, f6, f0 |
- ;; C 04 |
- xma.l f8 = f11, f12, f7 C q = q + nh |
- ;; C 08 |
- getf.sig r18 = f8 |
- xma.hu f9 = f8, f10, f0 |
- xma.l f8 = f8, f10, f0 |
- ;; C 12 |
- getf.sig r16 = f9 |
- C 13 |
- getf.sig r15 = f8 |
- ;; C 18 |
- cmp.ltu p6, p7 = r20, r15 |
- sub r15 = r20, r15 |
- sub r16 = r38, r16 |
- ;; C 19 |
- (p6) cmp.ne p8, p9 = 1, r16 C is rH != 0? |
- (p7) cmp.ne p8, p9 = 0, r16 C is rH != 0? |
- (p6) add r16 = -1, r16 |
- (p0) cmp.ne.unc p6, p7 = r0, r0 |
- ;; C 20 |
- (p8) cmp.ltu p6, p7 = r15, r36 |
- (p8) sub r15 = r15, r36 |
- (p8) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; C 21 |
- .pred.rel "mutex",p6,p7 |
- (p6) cmp.ne p8, p9 = 1, r16 C is rH != 0 still? |
- (p7) cmp.ne p8, p9 = 0, r16 C is rH != 0 still? |
- cmp.ltu p6, p7 = r15, r36 C speculative |
- sub r28 = r15, r36 C speculative, just for cmp |
- ;; C 22 |
- (p8) cmp.ltu p6, p7 = r28, r36 C redo last cmp if needed |
- (p8) mov r15 = r28 |
- (p8) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; C 23 |
- (p6) setf.sig f7 = r15 |
- (p7) sub r15 = r15, r36 |
- (p7) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; C 24 |
- (p7) setf.sig f7 = r15 |
- st8 [r32] = r18, -8 |
- mov r38 = r15 |
- br.cloop.dptk .Loop1 |
- C 29/30 |
- br.sptk .L435 |
- ;; |
-.Lunnorm: |
- mux1 r16 = r36, @rev |
- cmp.eq p6, p7 = 0, r35 |
- (p6) br.cond.dpnt .L322 |
- shladd r34 = r35, 3, r34 |
- ;; |
- adds r34 = -8, r34 |
- ;; |
- ld8 r39 = [r34] |
- ;; |
- cmp.leu p6, p7 = r36, r39 |
- (p6) br.cond.dptk .L322 |
- adds r34 = -8, r34 |
- ;; |
- mov r38 = r39 |
- ;; |
- cmp.ne p6, p7 = 1, r15 |
- st8 [r32] = r0, -8 |
- ;; |
- (p7) mov r8 = r38 |
- (p7) br.cond.dpnt .Lret |
- adds r35 = -1, r35 |
-.L322: |
- sub r14 = r0, r16 |
- ;; |
- or r14 = r16, r14 |
- ;; |
- mov r16 = -8 |
- czx1.l r14 = r14 |
- ;; |
- shladd r16 = r14, 3, r16 |
- ;; |
- shr.u r14 = r36, r16 |
- ;; |
- cmp.geu p6, p7 = 15, r14 |
- ;; |
- (p7) shr.u r14 = r14, 4 |
- (p7) adds r16 = 4, r16 |
- ;; |
- cmp.geu p6, p7 = 3, r14 |
- ;; |
- (p7) shr.u r14 = r14, 2 |
- (p7) adds r16 = 2, r16 |
- ;; |
- tbit.nz p6, p7 = r14, 1 |
- ;; |
- .pred.rel "mutex",p6,p7 |
- (p6) sub r40 = 62, r16 |
- (p7) sub r40 = 63, r16 |
- ;; |
- shl r45 = r36, r40 |
- shl r36 = r36, r40 |
- shl r38 = r38, r40 |
- br.call.sptk.many b0 = mpn_invert_limb |
- ;; |
-.Lpu: |
- mov r23 = 1 |
- ;; |
- setf.sig f6 = r8 |
- setf.sig f12 = r23 |
- cmp.eq p6, p7 = 0, r35 |
- (p6) br.cond.dpnt .L435 |
- sub r16 = 64, r40 |
- adds r35 = -2, r35 |
- ;; |
- ld8 r39 = [r34], -8 |
- cmp.le p6, p7 = 0, r35 |
- ;; |
- shr.u r14 = r39, r16 |
- ;; |
- or r38 = r14, r38 |
- (p7) br.cond.dpnt .Lend3 |
- ;; |
- mov r22 = r16 |
- setf.sig f10 = r36 |
- setf.sig f7 = r38 |
- mov ar.lc = r35 |
- ;; |
-C Develop quotient limbs for unnormalized divisor |
-.Loop3: |
- ld8 r14 = [r34], -8 |
- xma.hu f11 = f7, f6, f0 |
- ;; |
- xma.l f8 = f11, f12, f7 C q = q + nh |
- ;; |
- getf.sig r18 = f8 |
- xma.hu f9 = f8, f10, f0 |
- shl r20 = r39, r40 |
- xma.l f8 = f8, f10, f0 |
- shr.u r24 = r14, r22 |
- ;; |
- getf.sig r16 = f9 |
- getf.sig r15 = f8 |
- or r20 = r24, r20 |
- ;; |
- cmp.ltu p6, p7 = r20, r15 |
- sub r15 = r20, r15 |
- sub r16 = r38, r16 |
- ;; |
- (p6) cmp.ne p8, p9 = 1, r16 C is rH != 0? |
- (p7) cmp.ne p8, p9 = 0, r16 C is rH != 0? |
- (p6) add r16 = -1, r16 |
- (p0) cmp.ne.unc p6, p7 = r0, r0 |
- ;; |
- (p8) cmp.ltu p6, p7 = r15, r36 |
- (p8) sub r15 = r15, r36 |
- (p8) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; |
- .pred.rel "mutex",p6,p7 |
- (p6) cmp.ne p8, p9 = 1, r16 C is rH != 0 still? |
- (p7) cmp.ne p8, p9 = 0, r16 C is rH != 0 still? |
- cmp.ltu p6, p7 = r15, r36 C speculative |
- sub r28 = r15, r36 C speculative, just for cmp |
- ;; |
- (p8) cmp.ltu p6, p7 = r28, r36 C redo last cmp if needed |
- (p8) mov r15 = r28 |
- (p8) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; |
- (p6) setf.sig f7 = r15 |
- (p7) sub r15 = r15, r36 |
- (p7) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; |
- (p7) setf.sig f7 = r15 |
- st8 [r32] = r18, -8 |
- mov r39 = r14 |
- mov r38 = r15 |
- br.cloop.dptk .Loop3 |
- ;; |
-.Lend3: |
- setf.sig f10 = r36 |
- setf.sig f7 = r38 |
- ;; |
- xma.hu f11 = f7, f6, f0 |
- ;; |
- xma.l f8 = f11, f12, f7 C q = q + nh |
- ;; |
- getf.sig r18 = f8 |
- xma.hu f9 = f8, f10, f0 |
- shl r20 = r39, r40 |
- xma.l f8 = f8, f10, f0 |
- ;; |
- getf.sig r16 = f9 |
- getf.sig r15 = f8 |
- ;; |
- cmp.ltu p6, p7 = r20, r15 |
- sub r15 = r20, r15 |
- sub r16 = r38, r16 |
- ;; |
- (p6) cmp.ne p8, p9 = 1, r16 C is rH != 0? |
- (p7) cmp.ne p8, p9 = 0, r16 C is rH != 0? |
- (p6) add r16 = -1, r16 |
- (p0) cmp.ne.unc p6, p7 = r0, r0 |
- ;; |
- (p8) cmp.ltu p6, p7 = r15, r36 |
- (p8) sub r15 = r15, r36 |
- (p8) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; |
- .pred.rel "mutex",p6,p7 |
- (p6) cmp.ne p8, p9 = 1, r16 C is rH != 0 still? |
- (p7) cmp.ne p8, p9 = 0, r16 C is rH != 0 still? |
- ;; |
- (p8) sub r15 = r15, r36 |
- (p8) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; |
- cmp.ltu p6, p7 = r15, r36 |
- ;; |
- (p7) sub r15 = r15, r36 |
- (p7) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; |
- st8 [r32] = r18, -8 |
- mov r38 = r15 |
-.L435: |
- adds r35 = -1, r33 |
- cmp.le p6, p7 = 1, r33 |
- (p7) br.cond.dpnt .Lend4 |
- ;; |
- setf.sig f7 = r38 |
- setf.sig f10 = r36 |
- mov ar.lc = r35 |
- ;; |
-.Loop4: |
- xma.hu f11 = f7, f6, f0 |
- ;; |
- xma.l f8 = f11, f12, f7 C q = q + nh |
- ;; |
- getf.sig r18 = f8 |
- xma.hu f9 = f8, f10, f0 |
- xma.l f8 = f8, f10, f0 |
- ;; |
- getf.sig r16 = f9 |
- getf.sig r15 = f8 |
- ;; |
- cmp.ltu p6, p7 = 0, r15 |
- sub r15 = 0, r15 |
- sub r16 = r38, r16 |
- ;; |
- (p6) cmp.ne p8, p9 = 1, r16 C is rH != 0? |
- (p7) cmp.ne p8, p9 = 0, r16 C is rH != 0? |
- (p6) add r16 = -1, r16 |
- (p0) cmp.ne.unc p6, p7 = r0, r0 |
- ;; |
- (p8) cmp.ltu p6, p7 = r15, r36 |
- (p8) sub r15 = r15, r36 |
- (p8) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; |
- .pred.rel "mutex",p6,p7 |
- (p6) cmp.ne p8, p9 = 1, r16 C is rH != 0 still? |
- (p7) cmp.ne p8, p9 = 0, r16 C is rH != 0 still? |
- cmp.ltu p6, p7 = r15, r36 C speculative |
- sub r28 = r15, r36 C speculative, just for cmp |
- ;; |
- (p8) cmp.ltu p6, p7 = r28, r36 C redo last cmp if needed |
- (p8) mov r15 = r28 |
- (p8) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; |
- (p6) setf.sig f7 = r15 |
- (p7) sub r15 = r15, r36 |
- (p7) add r18 = 1, r18 C q = q + 1; done if: rH > 0 |
- ;; |
- (p7) setf.sig f7 = r15 |
- st8 [r32] = r18, -8 |
- mov r38 = r15 |
- br.cloop.dptk .Loop4 |
- ;; |
-.Lend4: |
- shr.u r8 = r38, r40 |
-.Lret: |
- mov ar.pfs = r42 |
- mov ar.lc = r44 |
- mov b0 = r41 |
- br.ret.sptk.many b0 |
-EPILOGUE() |
-ASM_END() |