Date: Sat, 17 Feb 2024 15:00:39 +0000

Double-wide remainders are important for safe modular arithmetic. For

example, std::linear_congruential_generator is commonly implemented

using double-wide arithmetic through __int128. However, it would be

desirable to have a general solution for any integer type.

Your div_wide function does not cover this use case because it has the

precondition is_div_wide_defined. A rem_wide function would have a

wider contract. Only division by zero and the edge case of HUGE / -1

are problematic, but the remainder always fits into a 64-bit number.

> FWIW, for std::linear_congruential_generator the integers involved are unsigned and the remainder is a compile-time constant that can't be zero, so neither of those cases is a problem. That doesn't mean it couldn't be useful in other situations, but isn't needed for std::linear_congruential_generator.

That was my thought, being a modular math, I had assumed it was unsigned, and for unsigned you can extend the domain of usability for the remainder using the following algorithm:

//assume dividend != 0

div_result<uint64_t> temporary = dive_wide(0, divisor_hi, dividend);

temporary = dive_wide(temporary.remainder, divisor_low, dividend);

return temporary.remainder;

Thatâ€™s why these are so useful.

Received on 2024-02-17 15:00:43