1/* Pentium __mpn_submul_1 -- Multiply a limb vector with a limb and subtract
2   the result from a second limb vector.
3   Copyright (C) 1992, 94, 96, 97, 98, 00 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, write to
18   the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
19   MA 02111-1307, USA. */
20
21#include "sysdep.h"
22#include "asm-syntax.h"
23#include "bp-sym.h"
24#include "bp-asm.h"
25
26#define PARMS	LINKAGE+16	/* space for 4 saved regs */
27#define RES	PARMS
28#define S1	RES+PTR_SIZE
29#define SIZE	S1+PTR_SIZE
30#define S2LIMB	SIZE+4
31
32#define res_ptr edi
33#define s1_ptr esi
34#define size ecx
35#define s2_limb ebx
36
37	.text
38ENTRY (BP_SYM (__mpn_submul_1))
39	ENTER
40
41	pushl	%edi
42	pushl	%esi
43	pushl	%ebp
44	pushl	%ebx
45
46	movl	RES(%esp), %res_ptr
47	movl	S1(%esp), %s1_ptr
48	movl	SIZE(%esp), %size
49	movl	S2LIMB(%esp), %s2_limb
50#if __BOUNDED_POINTERS__
51	shll	$2, %sizeP	/* convert limbs to bytes */
52	CHECK_BOUNDS_BOTH_WIDE (%res_ptr, RES(%esp), %sizeP)
53	CHECK_BOUNDS_BOTH_WIDE (%s1_ptr, S1(%esp), %sizeP)
54	shrl	$2, %sizeP
55#endif
56	leal	(%res_ptr,%size,4), %res_ptr
57	leal	(%s1_ptr,%size,4), %s1_ptr
58	negl	%size
59	xorl	%ebp, %ebp
60	ALIGN (3)
61
62L(oop):	adcl	$0, %ebp
63	movl	(%s1_ptr,%size,4), %eax
64
65	mull	%s2_limb
66
67	addl	%ebp, %eax
68	movl	(%res_ptr,%size,4), %ebp
69
70	adcl	$0, %edx
71	subl	%eax, %ebp
72
73	movl	%ebp, (%res_ptr,%size,4)
74	incl	%size
75
76	movl	%edx, %ebp
77	jnz	L(oop)
78
79	adcl	$0, %ebp
80	movl	%ebp, %eax
81	popl	%ebx
82	popl	%ebp
83	popl	%esi
84	popl	%edi
85
86	LEAVE
87	ret
88#undef size
89END (BP_SYM (__mpn_submul_1))
90