constant_time_msb_bn
return constant_time_msb_bn(a ^ ((a ^ b) | ((a - b) ^ b)));
return constant_time_msb_bn(~a & (a - 1));