Subversion Repositories planix.SVN

Rev

Rev 33 | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
2 - 1 
#include "os.h"
2 
#include <mp.h>
3 
#include "dat.h"
4 
 
5 
// sum = abs(b1) + abs(b2), i.e., add the magnitudes
6 
void
7 
mpmagadd(mpint *b1, mpint *b2, mpint *sum)
8 
{
9 
	int m, n;
10 
	mpint *t;
11 
 
33 7u83 12 
	sum->flags |= (b1->flags | b2->flags) & MPtimesafe;
13 
 
2 - 14 
	// get the sizes right
15 
	if(b2->top > b1->top){
16 
		t = b1;
17 
		b1 = b2;
18 
		b2 = t;
19 
	}
20 
	n = b1->top;
21 
	m = b2->top;
22 
	if(n == 0){
23 
		mpassign(mpzero, sum);
24 
		return;
25 
	}
26 
	if(m == 0){
27 
		mpassign(b1, sum);
33 7u83 28 
		sum->sign = 1;
2 - 29 
		return;
30 
	}
31 
	mpbits(sum, (n+1)*Dbits);
32 
	sum->top = n+1;
33 
 
34 
	mpvecadd(b1->p, n, b2->p, m, sum->p);
35 
	sum->sign = 1;
36 
 
37 
	mpnorm(sum);
38 
}
39 
 
40 
// sum = b1 + b2
41 
void
42 
mpadd(mpint *b1, mpint *b2, mpint *sum)
43 
{
44 
	int sign;
45 
 
46 
	if(b1->sign != b2->sign){
33 7u83 47 
		assert(((b1->flags | b2->flags | sum->flags) & MPtimesafe) == 0);
2 - 48 
		if(b1->sign < 0)
49 
			mpmagsub(b2, b1, sum);
50 
		else
51 
			mpmagsub(b1, b2, sum);
52 
	} else {
53 
		sign = b1->sign;
54 
		mpmagadd(b1, b2, sum);
55 
		if(sum->top != 0)
56 
			sum->sign = sign;
57 
	}
58 
}