si_mult_sat
int64_t ll = si_mult_sat(tp, a.smin, b.smin);
int64_t lu = si_mult_sat(tp, a.smin, b.smax);
int64_t ul = si_mult_sat(tp, a.smax, b.smin);
int64_t uu = si_mult_sat(tp, a.smax, b.smax);