blob: 894d7f80e4d5f296ef903d9e0921b27c3a78c70d (
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
|
#include "os.h"
#include <mp.h>
#include "dat.h"
int
mpvectscmp(mpdigit *a, int alen, mpdigit *b, int blen)
{
mpdigit x, y, z, v;
int m, p;
if(alen > blen){
v = 0;
while(alen > blen)
v |= a[--alen];
m = p = ((-v^v)|v)>>(Dbits-1);
} else if(blen > alen){
v = 0;
while(blen > alen)
v |= b[--blen];
m = ((-v^v)|v)>>(Dbits-1);
p = m^1;
} else
m = p = 0;
while(alen-- > 0){
x = a[alen];
y = b[alen];
z = x - y;
x = ~x;
v = (((-z^z)|z)>>(Dbits-1)) & ~m;
p = ((~((x&y)|(x&z)|(y&z))>>(Dbits-1)) & v) | (p & ~v);
m |= v;
}
return (p-m) | m;
}
|
