1# Pyramid __gmpn_addmul_1 -- Multiply a limb vector with a limb and add 2# the result to a second limb vector. 3 4# Copyright 1995, 2000 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 3 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. If not, see http://www.gnu.org/licenses/. 20 21.text 22 .align 2 23.globl ___gmpn_addmul_1 24___gmpn_addmul_1: 25 mova (pr0)[pr2*4],pr0 26 mova (pr1)[pr2*4],pr1 27 mnegw pr2,pr2 28 movw $0,tr3 29 30Loop: movw (pr1)[pr2*4],tr1 31 uemul pr3,tr0 32 addw tr3,tr1 33 movw $0,tr3 34 addwc tr0,tr3 35 movw (pr0)[pr2*0x4],tr0 36 addw tr0,tr1 37 addwc $0,tr3 38 movw tr1,(pr0)[pr2*4] 39 addw $1,pr2 40 bne Loop 41 42 movw tr3,pr0 43 ret 44