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
|
#include "os.h"
#include <mp.h>
#include "dat.h"
// return neg, 0, pos as abs(b1)-abs(b2) is neg, 0, pos
int
mpmagcmp(mpint *b1, mpint *b2)
{
int i;
i = b1->flags | b2->flags;
if(i & MPtimesafe)
return mpvectscmp(b1->p, b1->top, b2->p, b2->top);
if(i & MPnorm){
i = b1->top - b2->top;
if(i)
return i;
}
return mpveccmp(b1->p, b1->top, b2->p, b2->top);
}
// return neg, 0, pos as b1-b2 is neg, 0, pos
int
mpcmp(mpint *b1, mpint *b2)
{
int sign;
sign = (b1->sign - b2->sign) >> 1; // -1, 0, 1
return sign | (((sign&1)-1) & mpmagcmp(b1, b2)*b1->sign);
}
|
