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

Diff for /OpenXM_contrib/gmp/mpn/powerpc32/Attic/sub_n.asm between version 1.1.1.1 and 1.1.1.2

version 1.1.1.1, 2000/09/09 14:12:38 version 1.1.1.2, 2003/08/25 16:06:24
Line 1 
Line 1 
 dnl PowerPC-32 mpn_sub_n -- Subtract two limb vectors of the same length > 0  dnl  PowerPC 750 mpn_sub_n -- subtract limb vectors.
 dnl and store difference in a third limb vector.  
   
 dnl Copyright (C) 1995, 1997, 2000 Free Software Foundation, Inc.  dnl  Copyright 2002 Free Software Foundation, Inc.
   dnl
   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
   dnl  modify it under the terms of the GNU Lesser General Public License as
   dnl  published by the Free Software Foundation; either version 2.1 of the
   dnl  License, or (at your option) any later version.
   dnl
   dnl  The GNU MP Library is distributed in the hope that it will be useful,
   dnl  but WITHOUT ANY WARRANTY; without even the implied warranty of
   dnl  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
   dnl  Lesser General Public License for more details.
   dnl
   dnl  You should have received a copy of the GNU Lesser General Public
   dnl  License along with the GNU MP Library; see the file COPYING.LIB.  If
   dnl  not, write to the Free Software Foundation, Inc., 59 Temple Place -
   dnl  Suite 330, Boston, MA 02111-1307, USA.
   
 dnl This file is part of the GNU MP Library.  include(`../config.m4')
   
 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 the Free Software Foundation; either version 2.1 of the License, or (at your  
 dnl option) any later version.  
   
 dnl The GNU MP Library is distributed in the hope that it will be useful, but  C       cycles/limb
 dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY  C 604e:     4.0
 dnl or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public  C 750:      4.0
 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; see the file COPYING.LIB.  If not, write to  
 dnl the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,  
 dnl MA 02111-1307, USA.  
   
   C mp_limb_t mpn_sub_n (mp_ptr dst, mp_srcptr src1, mp_srcptr src2,
   C                      mp_size_t size);
   C
   C Same style as mpn_add_n.
   
 dnl INPUT PARAMETERS  
 dnl res_ptr     r3  
 dnl s1_ptr      r4  
 dnl s2_ptr      r5  
 dnl size        r6  
   
 include(`../config.m4')  
   
 ASM_START()  ASM_START()
 PROLOGUE(mpn_sub_n)  PROLOGUE(mpn_sub_n)
         mtctr   r6              C copy size into CTR  
         addic   r0,r6,-1        C set cy  
         lwz     r8,0(r4)        C load least significant s1 limb  
         lwz     r0,0(r5)        C load least significant s2 limb  
         addi    r3,r3,-4        C offset res_ptr, it's updated before it's used  
         bdz     .Lend           C If done, skip loop  
 .Loop:  lwz     r9,4(r4)        C load s1 limb  
         lwz     r10,4(r5)       C load s2 limb  
         subfe   r7,r0,r8        C subtract limbs with cy, set cy  
         stw     r7,4(r3)        C store result limb  
         bdz     .Lexit          C decrement CTR and exit if done  
         lwzu    r8,8(r4)        C load s1 limb and update s1_ptr  
         lwzu    r0,8(r5)        C load s2 limb and update s2_ptr  
         subfe   r7,r10,r9       C subtract limbs with cy, set cy  
         stwu    r7,8(r3)        C store result limb and update res_ptr  
         bdnz    .Loop           C decrement CTR and loop back  
   
 .Lend:  subfe   r7,r0,r8          C r3    wp
         stw     r7,4(r3)        C store ultimate result limb          C r4    xp
         subfe   r3,r0,r0        C load !cy into ...          C r5    yp
         subfic  r3,r3,0         C ... return value register          C r6    size
   
           mtctr   r6              C size
           lwz     r6, 0(r4)       C xp[0]
   
           lwz     r7, 0(r5)       C yp[0]
           sub     r4, r4, r3      C xp-wp
   
           sub     r5, r5, r3      C yp-wp
           addi    r4, r4, 4       C xp-wp-4
   
           addi    r5, r5, 4       C yp-wp-4
           subfc   r6, r7, r6
   
           stw     r6, 0(r3)       C wp[0]
           bdz     L(done)
   
   
   L(top):
           C r3    wp, incrementing
           C r4    xp-wp+4
           C r5    yp-wp+4
           C r6    scratch
           C r7    scratch
   
           lwzx    r6, r4,r3       C xp[i]
   
           lwzx    r7, r5,r3       C yp[i]
   
           C serialize
           subfe   r6, r7, r6
   
           stwu    r6, 4(r3)       C wp[i]
           bdnz    L(top)
   
   
   L(done):
           subfe   r3, r0, r0      C 0 or -1
   
           subfic  r3, r3, 0       C 0 or 1
         blr          blr
 .Lexit: subfe   r7,r10,r9  
         stw     r7,8(r3)  EPILOGUE()
         subfe   r3,r0,r0        C load !cy into ...  
         subfic  r3,r3,0         C ... return value register  
         blr  
 EPILOGUE(mpn_sub_n)  
Legend:
Removed from v.1.1.1.1  
changed lines
  Added in v.1.1.1.2
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>