aors_n.asm
上传用户:qaz666999
上传日期:2022-08-06
资源大小:2570k
文件大小:15k
- dnl IA-64 mpn_add_n/mpn_sub_n -- mpn addition and subtraction.
- dnl Copyright 2003, 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: 2.67
- C Itanium 2: 1.25
- C TODO
- C * Consider using special code for small n, using something like
- C "switch (8 * (n >= 8) + (n mod 8))" to enter it and feed-in code.
- C INPUT PARAMETERS
- define(`rp',`r32')
- define(`up',`r33')
- define(`vp',`r34')
- define(`n',`r35')
- ifdef(`OPERATION_add_n',`
- define(ADDSUB, add)
- define(PRED, ltu)
- define(INCR, 1)
- define(LIM, -1)
- define(func, mpn_add_n)
- ')
- ifdef(`OPERATION_sub_n',`
- define(ADDSUB, sub)
- define(PRED, gtu)
- define(INCR, -1)
- define(LIM, 0)
- define(func, mpn_sub_n)
- ')
- C Some useful aliases for registers we use
- define(`u0',`r14') define(`u1',`r15') define(`u2',`r16') define(`u3',`r17')
- define(`u4',`r18') define(`u5',`r19') define(`u6',`r20') define(`u7',`r21')
- define(`v0',`r24') define(`v1',`r25') define(`v2',`r26') define(`v3',`r27')
- define(`v4',`r28') define(`v5',`r29') define(`v6',`r30') define(`v7',`r31')
- define(`w0',`r22') define(`w1',`r9') define(`w2',`r8') define(`w3',`r23')
- define(`w4',`r22') define(`w5',`r9') define(`w6',`r8') define(`w7',`r23')
- define(`rpx',`r3')
- MULFUNC_PROLOGUE(mpn_add_n mpn_sub_n)
- ASM_START()
- PROLOGUE(func)
- .prologue
- .save ar.lc, r2
- .body
- ifdef(`HAVE_ABI_32',`
- addp4 rp = 0, rp C M I
- addp4 up = 0, up C M I
- addp4 vp = 0, vp C M I
- zxt4 n = n C I
- ;;
- ')
- {.mmi C 00
- ld8 r11 = [vp], 8 C M01
- ld8 r10 = [up], 8 C M01
- mov.i r2 = ar.lc C I0
- }
- {.mmi
- and r14 = 7, n C M I
- cmp.lt p15, p14 = 8, n C M I
- add n = -8, n C M I
- ;;
- }
- {.mmi C 01
- cmp.eq p6, p0 = 1, r14 C M I
- cmp.eq p7, p0 = 2, r14 C M I
- cmp.eq p8, p0 = 3, r14 C M I
- }
- {.bbb
- (p6) br.dptk .Lb001 C B
- (p7) br.dptk .Lb010 C B
- (p8) br.dptk .Lb011 C B
- ;;
- }
- {.mmi C 02
- cmp.eq p9, p0 = 4, r14 C M I
- cmp.eq p10, p0 = 5, r14 C M I
- cmp.eq p11, p0 = 6, r14 C M I
- }
- {.bbb
- (p9) br.dptk .Lb100 C B
- (p10) br.dptk .Lb101 C B
- (p11) br.dptk .Lb110 C B
- ;;
- } C 03
- {.mmb
- cmp.eq p12, p0 = 7, r14 C M I
- add n = -1, n C loop count M I
- (p12) br.dptk .Lb111 C B
- }
- .Lb000: ld8 v2 = [vp], 8 C M01
- ld8 u2 = [up], 8 C M01
- add rpx = 8, rp C M I
- ;;
- ld8 v3 = [vp], 8 C M01
- ld8 u3 = [up], 8 C M01
- ADDSUB w1 = r10, r11 C M I
- ;;
- ld8 v4 = [vp], 8 C M01
- ld8 u4 = [up], 8 C M01
- cmp.PRED p7, p0 = w1, r10 C M I
- ;;
- ld8 v5 = [vp], 8 C M01
- ld8 u5 = [up], 8 C M01
- ADDSUB w2 = u2, v2 C M I
- ;;
- ld8 v6 = [vp], 8 C M01
- ld8 u6 = [up], 8 C M01
- cmp.PRED p8, p0 = w2, u2 C M I
- ;;
- ld8 v7 = [vp], 8 C M01
- ld8 u7 = [up], 8 C M01
- ADDSUB w3 = u3, v3 C M I
- ;;
- ld8 v0 = [vp], 8 C M01
- ld8 u0 = [up], 8 C M01
- cmp.PRED p9, p0 = w3, u3 C M I
- (p7) cmp.eq.or p8, p0 = LIM, w2 C M I
- (p7) add w2 = INCR, w2 C M I
- (p14) br.cond.dptk .Lcj8 C B
- ;;
- .grt8: ld8 v1 = [vp], 8 C M01
- ld8 u1 = [up], 8 C M01
- shr.u n = n, 3 C I0
- ;;
- add r11 = 512, vp
- ld8 v2 = [vp], 8 C M01
- add r10 = 512, up
- ld8 u2 = [up], 8 C M01
- nop.i 0
- nop.b 0
- ;;
- ld8 v3 = [vp], 8 C M01
- ld8 u3 = [up], 8 C M01
- mov.i ar.lc = n C I0
- br .LL000 C B
- .Lb001: add rpx = 16, rp C M I
- ADDSUB w0 = r10, r11 C M I
- (p15) br.cond.dpnt .grt1 C B
- ;;
- cmp.PRED p6, p0 = w0, r10 C M I
- mov r8 = 0 C M I
- br .Lcj1 C B
- .grt1: ld8 v1 = [vp], 8 C M01
- ld8 u1 = [up], 8 C M01
- shr.u n = n, 3 C I0
- ;;
- ld8 v2 = [vp], 8 C M01
- ld8 u2 = [up], 8 C M01
- cmp.ne p9, p0 = r0, r0 C read near Loop
- ;;
- ld8 v3 = [vp], 8 C M01
- ld8 u3 = [up], 8 C M01
- mov.i ar.lc = n C I0
- ;;
- ld8 v4 = [vp], 8 C M01
- ld8 u4 = [up], 8 C M01
- cmp.PRED p6, p0 = w0, r10 C M I
- ;;
- ld8 v5 = [vp], 8 C M01
- ld8 u5 = [up], 8 C M01
- ADDSUB w1 = u1, v1 C M I
- ;;
- ld8 v6 = [vp], 8 C M01
- ld8 u6 = [up], 8 C M01
- cmp.PRED p7, p0 = w1, u1 C M I
- ;;
- ld8 v7 = [vp], 8 C M01
- ld8 u7 = [up], 8 C M01
- ADDSUB w2 = u2, v2 C M I
- ;;
- add r11 = 512, vp
- ld8 v0 = [vp], 8 C M01
- add r10 = 512, up
- ld8 u0 = [up], 8 C M01
- br.cloop.dptk .Loop C B
- br .Lcj9 C B
- .Lb010: ld8 v0 = [vp], 8 C M01
- ld8 u0 = [up], 8 C M01
- add rpx = 24, rp C M I
- ADDSUB w7 = r10, r11 C M I
- (p15) br.cond.dpnt .grt2 C B
- ;;
- cmp.PRED p9, p0 = w7, r10 C M I
- ADDSUB w0 = u0, v0 C M I
- br .Lcj2 C B
- .grt2: ld8 v1 = [vp], 8 C M01
- ld8 u1 = [up], 8 C M01
- shr.u n = n, 3 C I0
- ;;
- ld8 v2 = [vp], 8 C M01
- ld8 u2 = [up], 8 C M01
- ;;
- ld8 v3 = [vp], 8 C M01
- ld8 u3 = [up], 8 C M01
- mov.i ar.lc = n C I0
- ;;
- ld8 v4 = [vp], 8 C M01
- ld8 u4 = [up], 8 C M01
- ;;
- ld8 v5 = [vp], 8 C M01
- ld8 u5 = [up], 8 C M01
- cmp.PRED p9, p0 = w7, r10 C M I
- ;;
- ld8 v6 = [vp], 8 C M01
- ld8 u6 = [up], 8 C M01
- ADDSUB w0 = u0, v0 C M I
- ;;
- add r11 = 512, vp
- ld8 v7 = [vp], 8 C M01
- add r10 = 512, up
- ld8 u7 = [up], 8 C M01
- br .LL01x C B
- .Lb011: ld8 v7 = [vp], 8 C M01
- ld8 u7 = [up], 8 C M01
- ADDSUB w6 = r10, r11 C M I
- ;;
- ld8 v0 = [vp], 8 C M01
- ld8 u0 = [up], 8 C M01
- (p15) br.cond.dpnt .grt3 C B
- ;;
- cmp.PRED p8, p0 = w6, r10 C M I
- ADDSUB w7 = u7, v7 C M I
- ;;
- st8 [rp] = w6, 8 C M23
- cmp.PRED p9, p0 = w7, u7 C M I
- br .Lcj3 C B
- .grt3: ld8 v1 = [vp], 8 C M01
- ld8 u1 = [up], 8 C M01
- add rpx = 32, rp C M I
- ;;
- ld8 v2 = [vp], 8 C M01
- ld8 u2 = [up], 8 C M01
- shr.u n = n, 3 C I0
- ;;
- ld8 v3 = [vp], 8 C M01
- ld8 u3 = [up], 8 C M01
- cmp.PRED p8, p0 = w6, r10 C M I
- ;;
- ld8 v4 = [vp], 8 C M01
- ld8 u4 = [up], 8 C M01
- mov.i ar.lc = n C I0
- ADDSUB w7 = u7, v7 C M I
- nop.i 0
- nop.b 0
- ;;
- ld8 v5 = [vp], 8 C M01
- ld8 u5 = [up], 8 C M01
- cmp.PRED p9, p0 = w7, u7 C M I
- ;;
- add r11 = 512, vp
- ld8 v6 = [vp], 8 C M01
- add r10 = 512, up
- ld8 u6 = [up], 8 C M01
- (p8) cmp.eq.or p9, p0 = LIM, w7 C M I
- ;;
- ld8 v7 = [vp], 8 C M01
- ld8 u7 = [up], 8 C M01
- (p8) add w7 = INCR, w7 C M I
- st8 [rp] = w6, 8 C M23
- ADDSUB w0 = u0, v0 C M I
- br .LL01x C B
- .Lb100: ld8 v6 = [vp], 8 C M01
- ld8 u6 = [up], 8 C M01
- add rpx = 8, rp C M I
- ;;
- ld8 v7 = [vp], 8 C M01
- ld8 u7 = [up], 8 C M01
- ADDSUB w5 = r10, r11 C M I
- ;;
- ld8 v0 = [vp], 8 C M01
- ld8 u0 = [up], 8 C M01
- (p15) br.cond.dpnt .grt4 C B
- ;;
- cmp.PRED p7, p0 = w5, r10 C M I
- ADDSUB w6 = u6, v6 C M I
- ;;
- cmp.PRED p8, p0 = w6, u6 C M I
- ADDSUB w7 = u7, v7 C M I
- br .Lcj4 C B
- .grt4: ld8 v1 = [vp], 8 C M01
- ld8 u1 = [up], 8 C M01
- shr.u n = n, 3 C I0
- cmp.PRED p7, p0 = w5, r10 C M I
- ;;
- ld8 v2 = [vp], 8 C M01
- ld8 u2 = [up], 8 C M01
- ADDSUB w6 = u6, v6 C M I
- ;;
- ld8 v3 = [vp], 8 C M01
- ld8 u3 = [up], 8 C M01
- cmp.PRED p8, p0 = w6, u6 C M I
- ;;
- ld8 v4 = [vp], 8 C M01
- ld8 u4 = [up], 8 C M01
- mov.i ar.lc = n C I0
- ;;
- ld8 v5 = [vp], 8 C M01
- ld8 u5 = [up], 8 C M01
- ADDSUB w7 = u7, v7 C M I
- ;;
- add r11 = 512, vp
- ld8 v6 = [vp], 8 C M01
- add r10 = 512, up
- ld8 u6 = [up], 8 C M01
- cmp.PRED p9, p0 = w7, u7 C M I
- ;;
- ld8 v7 = [vp], 8 C M01
- ld8 u7 = [up], 8 C M01
- (p7) cmp.eq.or p8, p0 = LIM, w6 C M I
- (p7) add w6 = INCR, w6 C M I
- br .LL100 C B
- .Lb101: ld8 v5 = [vp], 8 C M01
- ld8 u5 = [up], 8 C M01
- add rpx = 16, rp C M I
- ;;
- ld8 v6 = [vp], 8 C M01
- ld8 u6 = [up], 8 C M01
- ADDSUB w4 = r10, r11 C M I
- ;;
- ld8 v7 = [vp], 8 C M01
- ld8 u7 = [up], 8 C M01
- cmp.PRED p6, p0 = w4, r10 C M I
- ;;
- ld8 v0 = [vp], 8 C M01
- ld8 u0 = [up], 8 C M01
- ADDSUB w5 = u5, v5 C M I
- shr.u n = n, 3 C I0
- (p15) br.cond.dpnt .grt5 C B
- ;;
- cmp.PRED p7, p0 = w5, u5 C M I
- ADDSUB w6 = u6, v6 C M I
- br .Lcj5 C B
- .grt5: ld8 v1 = [vp], 8 C M01
- ld8 u1 = [up], 8 C M01
- ;;
- ld8 v2 = [vp], 8 C M01
- ld8 u2 = [up], 8 C M01
- mov.i ar.lc = n C I0
- ;;
- ld8 v3 = [vp], 8 C M01
- ld8 u3 = [up], 8 C M01
- cmp.PRED p7, p0 = w5, u5 C M I
- ;;
- ld8 v4 = [vp], 8 C M01
- ld8 u4 = [up], 8 C M01
- ADDSUB w6 = u6, v6 C M I
- ;;
- add r11 = 512, vp
- ld8 v5 = [vp], 8 C M01
- add r10 = 512, up
- ld8 u5 = [up], 8 C M01
- br .LL101 C B
- .Lb110: ld8 v4 = [vp], 8 C M01
- ld8 u4 = [up], 8 C M01
- add rpx = 24, rp C M I
- ;;
- ld8 v5 = [vp], 8 C M01
- ld8 u5 = [up], 8 C M01
- ADDSUB w3 = r10, r11 C M I
- ;;
- ld8 v6 = [vp], 8 C M01
- ld8 u6 = [up], 8 C M01
- shr.u n = n, 3 C I0
- ;;
- ld8 v7 = [vp], 8 C M01
- ld8 u7 = [up], 8 C M01
- cmp.PRED p9, p0 = w3, r10 C M I
- ;;
- ld8 v0 = [vp], 8 C M01
- ld8 u0 = [up], 8 C M01
- ADDSUB w4 = u4, v4 C M I
- (p14) br.cond.dptk .Lcj67 C B
- ;;
- .grt6: ld8 v1 = [vp], 8 C M01
- ld8 u1 = [up], 8 C M01
- mov.i ar.lc = n C I0
- cmp.PRED p9, p0 = w3, r10 C M I
- nop.i 0
- nop.b 0
- ;;
- ld8 v2 = [vp], 8 C M01
- ld8 u2 = [up], 8 C M01
- ADDSUB w4 = u4, v4 C M I
- ;;
- add r11 = 512, vp
- ld8 v3 = [vp], 8 C M01
- add r10 = 512, up
- ld8 u3 = [up], 8 C M01
- br .LL11x C B
- .Lb111: ld8 v3 = [vp], 8 C M01
- ld8 u3 = [up], 8 C M01
- add rpx = 32, rp C M I
- ;;
- ld8 v4 = [vp], 8 C M01
- ld8 u4 = [up], 8 C M01
- ADDSUB w2 = r10, r11 C M I
- ;;
- ld8 v5 = [vp], 8 C M01
- ld8 u5 = [up], 8 C M01
- cmp.PRED p8, p0 = w2, r10 C M I
- ;;
- ld8 v6 = [vp], 8 C M01
- ld8 u6 = [up], 8 C M01
- ADDSUB w3 = u3, v3 C M I
- ;;
- ld8 v7 = [vp], 8 C M01
- ld8 u7 = [up], 8 C M01
- cmp.PRED p9, p0 = w3, u3 C M I
- ;;
- ld8 v0 = [vp], 8 C M01
- ld8 u0 = [up], 8 C M01
- (p15) br.cond.dpnt .grt7 C B
- ;;
- st8 [rp] = w2, 8 C M23
- (p8) cmp.eq.or p9, p0 = LIM, w3 C M I
- (p8) add w3 = INCR, w3 C M I
- ADDSUB w4 = u4, v4 C M I
- br .Lcj67 C B
- .grt7: ld8 v1 = [vp], 8 C M01
- ld8 u1 = [up], 8 C M01
- shr.u n = n, 3 C I0
- (p8) cmp.eq.or p9, p0 = LIM, w3 C M I
- nop.i 0
- nop.b 0
- ;;
- add r11 = 512, vp
- ld8 v2 = [vp], 8 C M01
- add r10 = 512, up
- ld8 u2 = [up], 8 C M01
- (p8) add w3 = INCR, w3 C M I
- nop.b 0
- ;;
- ld8 v3 = [vp], 8 C M01
- ld8 u3 = [up], 8 C M01
- mov.i ar.lc = n C I0
- st8 [rp] = w2, 8 C M23
- ADDSUB w4 = u4, v4 C M I
- br .LL11x C B
- C *** MAIN LOOP START ***
- ALIGN(32)
- .Loop: ld8 v1 = [vp], 8 C M01
- cmp.PRED p7, p0 = w1, u1 C M I
- (p9) cmp.eq.or p6, p0 = LIM, w0 C M I
- ld8 u1 = [up], 8 C M01
- (p9) add w0 = INCR, w0 C M I
- ADDSUB w2 = u2, v2 C M I
- ;;
- ld8 v2 = [vp], 8 C M01
- cmp.PRED p8, p0 = w2, u2 C M I
- (p6) cmp.eq.or p7, p0 = LIM, w1 C M I
- ld8 u2 = [up], 8 C M01
- (p6) add w1 = INCR, w1 C M I
- ADDSUB w3 = u3, v3 C M I
- ;;
- st8 [rp] = w0, 8 C M23
- ld8 v3 = [vp], 8 C M01
- cmp.PRED p9, p0 = w3, u3 C M I
- (p7) cmp.eq.or p8, p0 = LIM, w2 C M I
- ld8 u3 = [up], 8 C M01
- (p7) add w2 = INCR, w2 C M I
- ;;
- .LL000: st8 [rp] = w1, 16 C M23
- st8 [rpx] = w2, 32 C M23
- (p8) cmp.eq.or p9, p0 = LIM, w3 C M I
- lfetch [r10], 64
- (p8) add w3 = INCR, w3 C M I
- ADDSUB w4 = u4, v4 C M I
- ;;
- .LL11x: st8 [rp] = w3, 8 C M23
- ld8 v4 = [vp], 8 C M01
- cmp.PRED p6, p0 = w4, u4 C M I
- ld8 u4 = [up], 8 C M01
- ADDSUB w5 = u5, v5 C M I
- ;;
- ld8 v5 = [vp], 8 C M01
- cmp.PRED p7, p0 = w5, u5 C M I
- (p9) cmp.eq.or p6, p0 = LIM, w4 C M I
- ld8 u5 = [up], 8 C M01
- (p9) add w4 = INCR, w4 C M I
- ADDSUB w6 = u6, v6 C M I
- ;;
- .LL101: ld8 v6 = [vp], 8 C M01
- cmp.PRED p8, p0 = w6, u6 C M I
- (p6) cmp.eq.or p7, p0 = LIM, w5 C M I
- ld8 u6 = [up], 8 C M01
- (p6) add w5 = INCR, w5 C M I
- ADDSUB w7 = u7, v7 C M I
- ;;
- st8 [rp] = w4, 8 C M23
- ld8 v7 = [vp], 8 C M01
- cmp.PRED p9, p0 = w7, u7 C M I
- (p7) cmp.eq.or p8, p0 = LIM, w6 C M I
- ld8 u7 = [up], 8 C M01
- (p7) add w6 = INCR, w6 C M I
- ;;
- .LL100: st8 [rp] = w5, 16 C M23
- st8 [rpx] = w6, 32 C M23
- (p8) cmp.eq.or p9, p0 = LIM, w7 C M I
- lfetch [r11], 64
- (p8) add w7 = INCR, w7 C M I
- ADDSUB w0 = u0, v0 C M I
- ;;
- .LL01x: st8 [rp] = w7, 8 C M23
- ld8 v0 = [vp], 8 C M01
- cmp.PRED p6, p0 = w0, u0 C M I
- ld8 u0 = [up], 8 C M01
- ADDSUB w1 = u1, v1 C M I
- br.cloop.dptk .Loop C B
- ;;
- C *** MAIN LOOP END ***
- cmp.PRED p7, p0 = w1, u1 C M I
- (p9) cmp.eq.or p6, p0 = LIM, w0 C M I
- (p9) add w0 = INCR, w0 C M I
- ADDSUB w2 = u2, v2 C M I
- ;;
- .Lcj9: cmp.PRED p8, p0 = w2, u2 C M I
- (p6) cmp.eq.or p7, p0 = LIM, w1 C M I
- st8 [rp] = w0, 8 C M23
- (p6) add w1 = INCR, w1 C M I
- ADDSUB w3 = u3, v3 C M I
- ;;
- cmp.PRED p9, p0 = w3, u3 C M I
- (p7) cmp.eq.or p8, p0 = LIM, w2 C M I
- (p7) add w2 = INCR, w2 C M I
- ;;
- .Lcj8: st8 [rp] = w1, 16 C M23
- st8 [rpx] = w2, 32 C M23
- (p8) cmp.eq.or p9, p0 = LIM, w3 C M I
- (p8) add w3 = INCR, w3 C M I
- ADDSUB w4 = u4, v4 C M I
- ;;
- .Lcj67: st8 [rp] = w3, 8 C M23
- cmp.PRED p6, p0 = w4, u4 C M I
- ADDSUB w5 = u5, v5 C M I
- ;;
- cmp.PRED p7, p0 = w5, u5 C M I
- (p9) cmp.eq.or p6, p0 = LIM, w4 C M I
- (p9) add w4 = INCR, w4 C M I
- ADDSUB w6 = u6, v6 C M I
- ;;
- .Lcj5: cmp.PRED p8, p0 = w6, u6 C M I
- (p6) cmp.eq.or p7, p0 = LIM, w5 C M I
- st8 [rp] = w4, 8 C M23
- (p6) add w5 = INCR, w5 C M I
- ADDSUB w7 = u7, v7 C M I
- ;;
- .Lcj4: cmp.PRED p9, p0 = w7, u7 C M I
- (p7) cmp.eq.or p8, p0 = LIM, w6 C M I
- (p7) add w6 = INCR, w6 C M I
- ;;
- st8 [rp] = w5, 16 C M23
- st8 [rpx] = w6, 32 C M23
- .Lcj3:
- (p8) cmp.eq.or p9, p0 = LIM, w7 C M I
- (p8) add w7 = INCR, w7 C M I
- ADDSUB w0 = u0, v0 C M I
- ;;
- .Lcj2: st8 [rp] = w7, 8 C M23
- cmp.PRED p6, p0 = w0, u0 C M I
- ;;
- (p9) cmp.eq.or p6, p0 = LIM, w0 C M I
- (p9) add w0 = INCR, w0 C M I
- mov r8 = 0 C M I
- ;;
- .Lcj1: st8 [rp] = w0, 8 C M23
- mov.i ar.lc = r2 C I0
- (p6) mov r8 = 1 C M I
- br.ret.sptk.many b0 C B
- EPILOGUE()
- ASM_END()