pub(crate) const fn mul_rem(a: Limb, b: Limb, d: NonZero<Limb>) -> Limb
Computes (a * b) % d.
(a * b) % d