Function PositiveDivide

src/cryptlib/integer.cpp:1876–1909

*/

void PositiveDivide(Integer &remainder, Integer &quotient,
                   const Integer &a, const Integer &b)
{
    unsigned aSize = a.WordCount();
    unsigned bSize = b.WordCount();

    if (!bSize)
    {
#ifdef THROW_EXCEPTIONS
        throw Integer::DivideByZero();
#else
        return;
#endif
    }

    if (a.PositiveCompare(b) == -1)
    {
        remainder = a;
        remainder.sign = Integer::POSITIVE;
        quotient = Integer::Zero();
        return;
    }

    aSize += aSize%2;   // round up to next even number
    bSize += bSize%2;

    remainder.reg.CleanNew(RoundupSize(bSize));
    remainder.sign = Integer::POSITIVE;
    quotient.reg.CleanNew(RoundupSize(aSize-bSize+2));
    quotient.sign = Integer::POSITIVE;

    SecWordBlock T(aSize+2*bSize+4);
    Divide(remainder.reg, quotient.reg, T, a.reg, aSize, b.reg, bSize);
}

void Integer::Divide(Integer &remainder, Integer &quotient, const Integer &dividend, const Integer &divisor)
{