Subversion Repositories planix.SVN

Rev

Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
26 7u83 1
/*
2
 *	mpvecadd(mpdigit *a, int alen, mpdigit *b, int blen, mpdigit *sum)
3
 *
4
 *		sum[0:alen] = a[0:alen-1] + b[0:blen-1]
5
 *
6
 *	prereq: alen >= blen, sum has room for alen+1 digits
7
 */
8
TEXT	mpvecadd(SB),$0
9
 
10
	MOVL	alen+8(FP),DX
11
	MOVL	blen+24(FP),CX
12
	MOVQ	RARG, SI
13
	MOVQ	b+16(FP),BX
14
	SUBL	CX,DX
15
	MOVQ	sum+32(FP),DI
16
	XORL	BP,BP			/* this also sets carry to 0 */
17
 
18
	/* skip addition if b is zero */
19
	TESTL	CX,CX
20
	JZ	_add1
21
 
22
	/* sum[0:blen-1],carry = a[0:blen-1] + b[0:blen-1] */
23
_addloop1:
24
	MOVL	(SI)(BP*4), AX
25
	ADCL	(BX)(BP*4), AX
26
	MOVL	AX,(DI)(BP*4)
27
	INCL	BP
28
	LOOP	_addloop1
29
 
30
_add1:
31
	/* jump if alen > blen */
32
	INCL	DX
33
	MOVL	DX,CX
34
	LOOP	_addloop2
35
 
36
	/* sum[alen] = carry */
37
 
38
_addend:
39
	ADCL	$0, CX
40
	MOVL	CX,(DI)(BP*4)
41
	RET
42
 
43
	/* sum[blen:alen-1],carry = a[blen:alen-1] + 0 */
44
_addloop2:
45
	MOVL	(SI)(BP*4),AX
46
	ADCL	$0,AX
47
	MOVL	AX,(DI)(BP*4)
48
	INCL	BP
49
	LOOP	_addloop2
50
	JMP	_addend
51