1/* i80386 __mpn_submul_1 -- Multiply a limb vector with a limb and subtract
2   the result from a second limb vector.
3   Copyright (C) 1992-2021 Free Software Foundation, Inc.
4   This file is part of the GNU MP Library.
5
6   The GNU MP Library is free software; you can redistribute it and/or modify
7   it under the terms of the GNU Lesser General Public License as published by
8   the Free Software Foundation; either version 2.1 of the License, or (at your
9   option) any later version.
10
11   The GNU MP Library is distributed in the hope that it will be useful, but
12   WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
13   or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
14   License for more details.
15
16   You should have received a copy of the GNU Lesser General Public License
17   along with the GNU MP Library; see the file COPYING.LIB.  If not,
18   see <https://www.gnu.org/licenses/>.  */
19
20#include "sysdep.h"
21#include "asm-syntax.h"
22
23#define PARMS	4+16	/* space for 4 saved regs */
24#define RES	PARMS
25#define S1	RES+4
26#define SIZE	S1+4
27#define S2LIMB	SIZE+4
28
29#define res_ptr edi
30#define s1_ptr esi
31#define sizeP ecx
32#define s2_limb ebx
33
34	.text
35ENTRY (__mpn_submul_1)
36
37	pushl	%edi
38	cfi_adjust_cfa_offset (4)
39	pushl	%esi
40	cfi_adjust_cfa_offset (4)
41	pushl	%ebp
42	cfi_adjust_cfa_offset (4)
43	pushl	%ebx
44	cfi_adjust_cfa_offset (4)
45	cfi_rel_offset (edi, 12)
46	cfi_rel_offset (esi, 8)
47	cfi_rel_offset (ebp, 4)
48	cfi_rel_offset (ebx, 0)
49
50	movl	RES(%esp), %res_ptr
51	movl	S1(%esp), %s1_ptr
52	movl	SIZE(%esp), %sizeP
53	movl	S2LIMB(%esp), %s2_limb
54	leal	(%res_ptr,%sizeP,4), %res_ptr
55	leal	(%s1_ptr,%sizeP,4), %s1_ptr
56	negl	%sizeP
57	xorl	%ebp, %ebp
58	ALIGN (3)
59L(oop):
60	movl	(%s1_ptr,%sizeP,4), %eax
61	mull	%s2_limb
62	addl	%ebp, %eax
63	adcl	$0, %edx
64	subl	%eax, (%res_ptr,%sizeP,4)
65	adcl	$0, %edx
66	movl	%edx, %ebp
67
68	incl	%sizeP
69	jnz	L(oop)
70	movl	%ebp, %eax
71
72	popl	%ebx
73	cfi_adjust_cfa_offset (-4)
74	cfi_restore (ebx)
75	popl	%ebp
76	cfi_adjust_cfa_offset (-4)
77	cfi_restore (ebp)
78	popl	%esi
79	cfi_adjust_cfa_offset (-4)
80	cfi_restore (esi)
81	popl	%edi
82	cfi_adjust_cfa_offset (-4)
83	cfi_restore (edi)
84
85	ret
86END (__mpn_submul_1)
87