glibc/sysdeps/m88k/m88110/mul_1.s

60 lines
1.6 KiB
ArmAsm
Raw Normal View History

1995-10-16 01:18:40 +00:00
; mc88110 __mpn_mul_1 -- Multiply a limb vector with a single limb and
; store the product in a second limb vector.
1996-03-01 18:43:45 +00:00
; Copyright (C) 1992, 1994, 1995 Free Software Foundation, Inc.
1995-10-16 01:18:40 +00:00
; This file is part of the GNU MP Library.
; The GNU MP Library is free software; you can redistribute it and/or modify
2001-07-07 10:13:33 +00:00
; it under the terms of the GNU Lesser General Public License as published by
; the Free Software Foundation; either version 2.1 of the License, or (at your
1995-10-16 01:18:40 +00:00
; option) any later version.
; The GNU MP Library is distributed in the hope that it will be useful, but
; WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
2001-07-07 10:13:33 +00:00
; or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
1995-10-16 01:18:40 +00:00
; License for more details.
2001-07-07 10:13:33 +00:00
; You should have received a copy of the GNU Lesser General Public License
1995-10-16 01:18:40 +00:00
; along with the GNU MP Library; see the file COPYING.LIB. If not, write to
1996-05-07 21:04:45 +00:00
; the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
; MA 02111-1307, USA.
1995-10-16 01:18:40 +00:00
; INPUT PARAMETERS
; res_ptr r2
; s1_ptr r3
; size r4
; s2_limb r5
text
align 16
global ___mpn_mul_1
___mpn_mul_1:
; Make S1_PTR and RES_PTR point at the end of their blocks
; and negate SIZE.
lda r3,r3[r4]
lda r8,r2[r4] ; RES_PTR in r8 since r2 is retval
subu r4,r0,r4
addu.co r2,r0,r0 ; r2 = cy = 0
ld r6,r3[r4]
addu r4,r4,1
mulu.d r10,r6,r5
bcnd.n eq0,r4,Lend
subu r8,r8,8
Loop: ld r6,r3[r4]
addu.cio r9,r11,r2
or r2,r10,r0 ; could be avoided if unrolled
addu r4,r4,1
mulu.d r10,r6,r5
bcnd.n ne0,r4,Loop
st r9,r8[r4]
Lend: addu.cio r9,r11,r2
st r9,r8,4
jmp.n r1
addu.ci r2,r10,r0