ref: 96c2da5fef670527ba6a7655917803bc22e83e79
dir: /3rd/mp/mpvectscmp.c/
#include "platform.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; }