pub fn div_nx1_normalized(u: &mut [u64], d: u64) -> u64
⚠️ Compute single limb normalized division.
The divisor must be normalized. See algorithm 7 from MG10.