sysdeps/hppa/sub_n.s

   1 ;! HP-PA  __mpn_sub_n -- Subtract two limb vectors of the same length > 0 and
   2 ;! store difference in a third limb vector.
   3
   4 ;! Copyright (C) 1992, 1994 Free Software Foundation, Inc.
   5
   6 ;! This file is part of the GNU MP Library.
   7
   8 ;! The GNU MP Library is free software; you can redistribute it and/or modify
   9 ;! it under the terms of the GNU Lesser General Public License as published by
  10 ;! the Free Software Foundation; either version 2.1 of the License, or (at your
  11 ;! option) any later version.
  12
  13 ;! The GNU MP Library is distributed in the hope that it will be useful, but
  14 ;! WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  15 ;! or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  16 ;! License for more details.
  17
  18 ;! You should have received a copy of the GNU Lesser General Public License
  19 ;! along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
  20 ;! the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
  21 ;! MA 02111-1307, USA.
  22
  23
  24 ;! INPUT PARAMETERS
  25 ;! res_ptr      gr26
  26 ;! s1_ptr       gr25
  27 ;! s2_ptr       gr24
  28 ;! size         gr23
  29
  30 ;! One might want to unroll this as for other processors, but it turns
  31 ;! out that the data cache contention after a store makes such
  32 ;! unrolling useless.  We can't come under 5 cycles/limb anyway.
  33
  34         .text
  35         .export         __mpn_sub_n
  36 __mpn_sub_n:
  37         .proc
  38         .callinfo       frame=0,no_calls
  39         .entry
  40
  41         ldws,ma         4(%r25),%r21
  42         ldws,ma         4(%r24),%r20
  43
  44         addib,=         -1,%r23,L$end   ;! check for (SIZE == 1)
  45          sub            %r21,%r20,%r28  ;! subtract first limbs ignoring cy
  46
  47 L$loop: ldws,ma         4(%r25),%r21
  48         ldws,ma         4(%r24),%r20
  49         stws,ma         %r28,4(%r26)
  50         addib,<>        -1,%r23,L$loop
  51          subb           %r21,%r20,%r28
  52
  53 L$end:  stws            %r28,0(%r26)
  54         addc            %r0,%r0,%r28
  55         bv              0(%r2)
  56          subi           1,%r28,%r28
  57
  58         .exit
  59         .procend