constant_time_msb
return constant_time_msb(a ^ ((a ^ b) | ((a - b) ^ b)));
return constant_time_msb(~a & (a - 1));
static inline unsigned int constant_time_msb(unsigned int a);