1 ;
! HP-PA __mpn_sub_n
-- Subtract two limb vectors of the same length
> 0 and
2 ;
! store difference in
a third limb vector.
4 ;
! Copyright
(C
) 1992, 1994 Free Software Foundation
, Inc.
6 ;
! This file is part of the GNU MP Library.
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.
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.
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.
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.
38 .callinfo frame=0,no_calls
41 ldws,ma 4(0,%r25),%r21
42 ldws,ma 4(0,%r24),%r20
44 addib,= -1,%r23,L$end ;! check for (SIZE == 1)
45 sub %r21,%r20,%r28 ;! subtract first limbs ignoring cy
47 L$loop: ldws,ma 4(0,%r25),%r21
48 ldws,ma 4(0,%r24),%r20
49 stws,ma %r28,4(0,%r26)
50 addib,<> -1,%r23,L$loop
53 L$end: stws %r28,0(0,%r26)