xref: /netbsd-src/external/lgpl3/gmp/dist/mpn/sparc64/ultrasparc1234/submul_1.asm (revision c34236556bea94afcaca1782d7d228301edc3ea0)
1dnl  SPARC v9 64-bit mpn_submul_1 -- Multiply a limb vector with a limb and
2dnl  subtract the result from a second limb vector.
3
4dnl  Copyright 2001, 2002, 2003 Free Software Foundation, Inc.
5
6dnl  This file is part of the GNU MP Library.
7
8dnl  The GNU MP Library is free software; you can redistribute it and/or modify
9dnl  it under the terms of the GNU Lesser General Public License as published
10dnl  by the Free Software Foundation; either version 3 of the License, or (at
11dnl  your option) any later version.
12
13dnl  The GNU MP Library is distributed in the hope that it will be useful, but
14dnl  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
15dnl  or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
16dnl  License for more details.
17
18dnl  You should have received a copy of the GNU Lesser General Public License
19dnl  along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.
20
21include(`../config.m4')
22
23C		   cycles/limb
24C UltraSPARC 1&2:     18
25C UltraSPARC 3:	      23
26
27C INPUT PARAMETERS
28C rp	i0
29C up	i1
30C n	i2
31C v	i3
32
33ASM_START()
34	REGISTER(%g2,#scratch)
35
36PROLOGUE(mpn_submul_1)
37	save	%sp,-176,%sp
38
39	sllx	%i2, 3, %g2
40	or	%g0, %i1, %o1
41	add	%g2, 15, %o0
42	or	%g0, %i2, %o2
43	and	%o0, -16, %o0
44	sub	%sp, %o0, %sp
45	add	%sp, 2223, %o0
46	or	%g0, %o0, %l0
47	call	mpn_mul_1
48	or	%g0, %i3, %o3
49	or	%g0, %o0, %l1		C preserve carry value from mpn_mul_1
50	or	%g0, %i0, %o0
51	or	%g0, %i0, %o1
52	or	%g0, %l0, %o2
53	call	mpn_sub_n
54	or	%g0, %i2, %o3
55	ret
56	restore	%l1, %o0, %o0		C sum carry values
57EPILOGUE(mpn_submul_1)
58