ML_DSA_D_BITS
*r1 = r >> ML_DSA_D_BITS; /* top 13 bits */
*r0 = r - (*r1 << ML_DSA_D_BITS); /* The remainder mod q */
r0_adjusted = mod_sub(*r0, 1 << ML_DSA_D_BITS);
mask = constant_time_lt((uint32_t)(1 << (ML_DSA_D_BITS - 1)), *r0);
out->coeff[i] = (in->coeff[i] << ML_DSA_D_BITS);