[BACK]Return to mul_1.asm CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gmp / mpn / alpha

Diff for /OpenXM_contrib/gmp/mpn/alpha/Attic/mul_1.asm between version 1.1.1.1 and 1.1.1.2

version 1.1.1.1, 2000/09/09 14:12:22 version 1.1.1.2, 2003/08/25 16:06:18
Line 1 
Line 1 
 dnl  Alpha __gmpn_mul_1 -- Multiply a limb vector with a limb and store  dnl  Alpha __gmpn_mul_1 -- Multiply a limb vector with a limb and store
 dnl  the result in a second limb vector.  dnl  the result in a second limb vector.
   
 dnl  Copyright (C) 1992, 1994, 1995, 2000 Free Software Foundation, Inc.  dnl  Copyright 1992, 1994, 1995, 2000, 2002 Free Software Foundation, Inc.
   
 dnl  This file is part of the GNU MP Library.  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  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 by  dnl  it under the terms of the GNU Lesser General Public License as published
 dnl  the Free Software Foundation; either version 2.1 of the License, or (at your  dnl  by the Free Software Foundation; either version 2.1 of the License, or (at
 dnl  option) any later version.  dnl  your option) any later version.
   
 dnl  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
 dnl  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY  dnl  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
Line 22  dnl  MA 02111-1307, USA.
Line 22  dnl  MA 02111-1307, USA.
   
 include(`../config.m4')  include(`../config.m4')
   
 dnl  INPUT PARAMETERS  C  INPUT PARAMETERS
 dnl  res_ptr    r16  C  rp   r16
 dnl  s1_ptr     r17  C  up   r17
 dnl  size       r18  C  n    r18
 dnl  s2_limb    r19  C  vl   r19
   C  cl   r20
   
 dnl  This code runs at 42 cycles/limb on EV4, 18 cycles/limb on EV5, and 7  C  This code runs at 42 cycles/limb on EV4, 18 cycles/limb on EV5, and 7
 dnl  cycles/limb on EV6.  C  cycles/limb on EV6.
   
 ASM_START()  ASM_START()
   PROLOGUE(mpn_mul_1c)
           ldq     r2,0(r17)       C r2 = s1_limb
           lda     r18,-1(r18)     C size--
           mulq    r2,r19,r3       C r3 = prod_low
           umulh   r2,r19,r4       C r4 = prod_high
           beq     r18,$Le1c       C jump if size was == 1
           ldq     r2,8(r17)       C r2 = s1_limb
           lda     r18,-1(r18)     C size--
           addq    r3,r20,r3       C r3 = cy_limb + cl
           stq     r3,0(r16)
           cmpult  r3,r20,r0       C r0 = carry from (cy_limb + cl)
           bne     r18,$Loop       C jump if size was == 2
           br      r31,$Le2
   $Le1c:  addq    r3,r20,r3       C r3 = cy_limb + cl
           cmpult  r3,r20,r0       C r0 = carry from (cy_limb + cl)
   $Le1:   stq     r3,0(r16)
           addq    r4,r0,r0
           ret     r31,(r26),1
   EPILOGUE(mpn_mul_1c)
   
 PROLOGUE(mpn_mul_1)  PROLOGUE(mpn_mul_1)
         ldq     r2,0(r17)       C r2 = s1_limb          ldq     r2,0(r17)       C r2 = s1_limb
         subq    r18,1,r18       C size--          lda     r18,-1(r18)     C size--
         mulq    r2,r19,r3       C r3 = prod_low          mulq    r2,r19,r3       C r3 = prod_low
         bic     r31,r31,r4      C clear cy_limb          bic     r31,r31,r0      C clear cy_limb
         umulh   r2,r19,r0       C r0 = prod_high          umulh   r2,r19,r4       C r4 = prod_high
         beq     r18,$Lend1      C jump if size was == 1          beq     r18,$Le1        C jump if size was == 1
         ldq     r2,8(r17)       C r2 = s1_limb          ldq     r2,8(r17)       C r2 = s1_limb
         subq    r18,1,r18       C size--          lda     r18,-1(r18)     C size--
         stq     r3,0(r16)          stq     r3,0(r16)
         beq     r18,$Lend2      C jump if size was == 2          beq     r18,$Le2        C jump if size was == 2
   
         ALIGN(8)          ALIGN(8)
 $Loop:  mulq    r2,r19,r3       C r3 = prod_low  $Loop:  mulq    r2,r19,r3       C r3 = prod_low
         addq    r4,r0,r0        C cy_limb = cy_limb + 'cy'          addq    r4,r0,r0        C cy_limb = cy_limb + 'cy'
         subq    r18,1,r18       C size--          lda     r18,-1(r18)     C size--
         umulh   r2,r19,r4       C r4 = cy_limb          umulh   r2,r19,r4       C r4 = prod_high
         ldq     r2,16(r17)      C r2 = s1_limb          ldq     r2,16(r17)      C r2 = s1_limb
         addq    r17,8,r17       C s1_ptr++          lda     r17,8(r17)      C s1_ptr++
         addq    r3,r0,r3        C r3 = cy_limb + prod_low          addq    r3,r0,r3        C r3 = cy_limb + prod_low
         stq     r3,8(r16)          stq     r3,8(r16)
         cmpult  r3,r0,r0        C r0 = carry from (cy_limb + prod_low)          cmpult  r3,r0,r0        C r0 = carry from (cy_limb + prod_low)
         addq    r16,8,r16       C res_ptr++          lda     r16,8(r16)      C res_ptr++
         bne     r18,$Loop          bne     r18,$Loop
   
 $Lend2: mulq    r2,r19,r3       C r3 = prod_low  $Le2:   mulq    r2,r19,r3       C r3 = prod_low
         addq    r4,r0,r0        C cy_limb = cy_limb + 'cy'          addq    r4,r0,r0        C cy_limb = cy_limb + 'cy'
         umulh   r2,r19,r4       C r4 = cy_limb          umulh   r2,r19,r4       C r4 = prod_high
         addq    r3,r0,r3        C r3 = cy_limb + prod_low          addq    r3,r0,r3        C r3 = cy_limb + prod_low
         cmpult  r3,r0,r0        C r0 = carry from (cy_limb + prod_low)          cmpult  r3,r0,r0        C r0 = carry from (cy_limb + prod_low)
         stq     r3,8(r16)          stq     r3,8(r16)
         addq    r4,r0,r0        C cy_limb = prod_high + cy          addq    r4,r0,r0        C cy_limb = prod_high + cy
         ret     r31,(r26),1  
 $Lend1: stq     r3,0(r16)  
         ret     r31,(r26),1          ret     r31,(r26),1
 EPILOGUE(mpn_mul_1)  EPILOGUE(mpn_mul_1)
 ASM_END()  ASM_END()
Legend:
Removed from v.1.1.1.1  
changed lines
  Added in v.1.1.1.2
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>