1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
; Clipper __gmpn_mul_1 -- Multiply a limb vector with a limb and store
; the result in a second limb vector.
; Copyright (C) 1995, 2000 Free Software Foundation, Inc.
; This file is part of the GNU MP Library.
; The GNU MP Library is free software; you can redistribute it and/or modify
; 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
; 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
; or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
; License for more details.
; You should have received a copy of the GNU Lesser General Public License
; along with the GNU MP Library; see the file COPYING.LIB. If not, write to
; the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
; MA 02111-1307, USA.
.text
.align 16
.globl ___gmpn_mul_1
___gmpn_mul_1:
subq $8,sp
storw r6,(sp)
loadw 12(sp),r2
loadw 16(sp),r3
loadq $0,r6 ; clear carry limb
.Loop: loadw (r1),r4
mulwux r3,r4
addw r6,r4 ; add old carry limb into low product limb
loadq $0,r6
addwc r5,r6 ; propagate cy into high product limb
storw r4,(r0)
addq $4,r0
addq $4,r1
subq $1,r2
brne .Loop
movw r6,r0
loadw 0(sp),r6
addq $8,sp
ret sp
|