35   T q2 = min / ad;               // Init. q2 = 2**p/|d|.  in SignedDivisionByConstant()  local
36   T r2 = min - q2 * ad;          // Init. r2 = rem(2**p, |d|).  in SignedDivisionByConstant()
46     q2 = 2 * q2;     // Update q2 = 2**p/|d|.  in SignedDivisionByConstant()
49       q2 = q2 + 1;  in SignedDivisionByConstant()
54   T mul = q2 + 1;  in SignedDivisionByConstant()
73   T q2 = max / d;         // Init. q2 = (2**p - 1)/d.  in UnsignedDivisionByConstant()  local
74   T r2 = max - q2 * d;    // Init. r2 = rem(2**p - 1, d).  in UnsignedDivisionByConstant()
86       if (q2 >= max) a = true;  in UnsignedDivisionByConstant()
87       q2 = 2 * q2 + 1;  in UnsignedDivisionByConstant()
90       if (q2 >= min) a = true;  in UnsignedDivisionByConstant()
91       q2 = 2 * q2;  in UnsignedDivisionByConstant()
96   return {q2 + 1, p - bits, a};  in UnsignedDivisionByConstant()