pub fn limb_mul(a: u64, b: u64) -> u128
The limb multiply mag_mul iterates: one 64×64→128 partial product.
mag_mul