code: drawterm

ref: 0f3766f44e554d73c7a3d8bab55cc32dc13f1c26
dir: /libmp/mpvectscmp.c/

View raw version
#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;
}