Annotation of OpenXM_contrib/gmp/mpn/sh/sh2/submul_1.s, Revision 1.1.1.3
1.1.1.2 maekawa 1: ! SH2 __gmpn_submul_1 -- Multiply a limb vector with a limb and subtract
1.1 maekawa 2: ! the result from a second limb vector.
3:
1.1.1.3 ! ohara 4: ! Copyright 1995, 2000 Free Software Foundation, Inc.
1.1 maekawa 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
1.1.1.2 maekawa 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
1.1 maekawa 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
1.1.1.2 maekawa 15: ! or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
1.1 maekawa 16: ! License for more details.
17:
1.1.1.2 maekawa 18: ! You should have received a copy of the GNU Lesser General Public License
1.1 maekawa 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.
22:
23:
24: ! INPUT PARAMETERS
25: ! res_ptr r4
26: ! s1_ptr r5
27: ! size r6
28: ! s2_limb r7
29:
30: .text
31: .align 1
1.1.1.2 maekawa 32: .global ___gmpn_submul_1
33: ___gmpn_submul_1:
1.1 maekawa 34: mov #0,r2 ! cy_limb = 0
35: mov #0,r0 ! Keep r0 = 0 for entire loop
36: clrt
37:
38: Loop: mov.l @r5+,r3
39: dmulu.l r3,r7
40: sts macl,r1
41: addc r2,r1 ! lo_prod += old cy_limb
42: sts mach,r2 ! new cy_limb = hi_prod
43: mov.l @r4,r3
44: addc r0,r2 ! cy_limb += T, T = 0
45: subc r3,r1
46: addc r0,r2 ! cy_limb += T, T = 0
47: dt r6
48: mov.l r1,@r4
49: bf.s Loop
50: add #4,r4
51:
52: rts
53: mov r2,r0
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>