aboutsummaryrefslogtreecommitdiff
path: root/src/libmp/386/mpvecadd-Darwin.s
blob: 2f68dbda19ff61d2f59027fe00af5618e0308e64 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
/*	mpvecadd(mpdigit *a, int alen, mpdigit *b, int blen, mpdigit *sum) */
/*		sum[0:alen] = a[0:alen-1] + b[0:blen-1] */
/*	prereq: alen >= blen, sum has room for alen+1 digits */
/* (very old gnu assembler doesn't allow multiline comments) */

.text

.p2align 2,0x90
.globl _mpvecadd
_mpvecadd:
	/* Prelude */
	pushl %ebp		/* save on stack */
	pushl %ebx
	pushl %esi
	pushl %edi

	leal 20(%esp), %ebp		/* %ebp = FP for now */

	movl	4(%ebp), %edx		/* alen */
	movl	12(%ebp), %ecx		/* blen */
	movl	0(%ebp), %esi		/* a */
	movl	8(%ebp), %ebx		/* b */
	subl	%ecx, %edx
	movl	16(%ebp), %edi		/* sum */
	xorl	%ebp, %ebp		/* this also sets carry to 0 */

	/* skip addition if b is zero */
	testl	%ecx,%ecx
	je	2f

	/* sum[0:blen-1],carry = a[0:blen-1] + b[0:blen-1] */
1:
	movl	(%esi, %ebp, 4), %eax
	adcl	(%ebx, %ebp, 4), %eax
	movl	%eax, (%edi, %ebp, 4)
	incl	%ebp
	loop	1b

2:
	/* jump if alen > blen */
	incl	%edx
	movl	%edx, %ecx
	loop	5f

	/* sum[alen] = carry */
3:
	jb	4f
	movl	$0, (%edi, %ebp, 4)
	jmp 6f

4:
	movl	$1, (%edi, %ebp, 4)
	jmp 6f

	/* sum[blen:alen-1],carry = a[blen:alen-1] + 0 */
5:
	movl	(%esi, %ebp, 4),%eax
	adcl	$0, %eax
	movl	%eax, (%edi, %ebp, 4)
	incl	%ebp
	loop	5b
	jmp	3b

6:
	/* Postlude */
	popl %edi
	popl %esi
	popl %ebx
	popl %ebp
	ret