46 lines
1.3 KiB
ArmAsm
46 lines
1.3 KiB
ArmAsm
|
# Pyramid __gmpn_addmul_1 -- Multiply a limb vector with a limb and add
|
||
|
# the result to a second limb vector.
|
||
|
|
||
|
# Copyright 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., 51 Franklin Street, Fifth Floor, Boston,
|
||
|
# MA 02110-1301, USA.
|
||
|
|
||
|
.text
|
||
|
.align 2
|
||
|
.globl ___gmpn_addmul_1
|
||
|
___gmpn_addmul_1:
|
||
|
mova (pr0)[pr2*4],pr0
|
||
|
mova (pr1)[pr2*4],pr1
|
||
|
mnegw pr2,pr2
|
||
|
movw $0,tr3
|
||
|
|
||
|
Loop: movw (pr1)[pr2*4],tr1
|
||
|
uemul pr3,tr0
|
||
|
addw tr3,tr1
|
||
|
movw $0,tr3
|
||
|
addwc tr0,tr3
|
||
|
movw (pr0)[pr2*0x4],tr0
|
||
|
addw tr0,tr1
|
||
|
addwc $0,tr3
|
||
|
movw tr1,(pr0)[pr2*4]
|
||
|
addw $1,pr2
|
||
|
bne Loop
|
||
|
|
||
|
movw tr3,pr0
|
||
|
ret
|