private static void ModDivCore(ref BigIntegerBuilder regNum, ref BigIntegerBuilder regDen, bool fQuo, ref BigIntegerBuilder regQuo)
{
regQuo.Set((uint)0);
if (regNum._iuLast >= regDen._iuLast)
{
int num = regDen._iuLast + 1;
int num2 = regNum._iuLast - regDen._iuLast;
int cu = num2;
int index = regNum._iuLast;
while (true)
{
if (index < num2)
{
cu++;
break;
}
if (regDen._rgu[index - num2] != regNum._rgu[index])
{
if (regDen._rgu[index - num2] < regNum._rgu[index])
{
cu++;
}
break;
}
index--;
}
if (cu != 0)
{
if (fQuo)
{
regQuo.SetSizeLazy(cu);
}
uint u = regDen._rgu[num - 1];
uint num6 = regDen._rgu[num - 2];
int num7 = NumericsHelpers.CbitHighZero(u);
int num8 = 0x20 - num7;
if (num7 > 0)
{
u = (u << num7) | (num6 >> num8);
num6 = num6 << num7;
if (num > 2)
{
num6 |= regDen._rgu[num - 3] >> num8;
}
}
regNum.EnsureWritable();
int num9 = cu;
while (--num9 >= 0)
{
uint uHi = ((num9 + num) <= regNum._iuLast) ? regNum._rgu[num9 + num] : 0;
ulong num11 = NumericsHelpers.MakeUlong(uHi, regNum._rgu[(num9 + num) - 1]);
uint uLo = regNum._rgu[(num9 + num) - 2];
if (num7 > 0)
{
num11 = (num11 << num7) | (uLo >> num8);
uLo = uLo << num7;
if ((num9 + num) >= 3)
{
uLo |= regNum._rgu[(num9 + num) - 3] >> num8;
}
}
ulong num13 = num11 / ((ulong)u);
ulong num14 = (uint)(num11 % ((ulong)u));
if (num13 > 0xffffffffL)
{
num14 += (ulong)u * (num13 - 0xffffffffL);
num13 = 0xffffffffL;
}
while ((num14 <= 0xffffffffL) && ((num13 * num6) > NumericsHelpers.MakeUlong((uint)num14, uLo)))
{
num13 -= (ulong)1L;
num14 += (ulong)u;
}
if (num13 > 0L)
{
ulong num15 = 0L;
for (int i = 0; i < num; i++)
{
num15 += regDen._rgu[i] * num13;
uint num17 = (uint)num15;
num15 = num15 >> 0x20;
if (regNum._rgu[num9 + i] < num17)
{
num15 += (ulong)1L;
}
regNum._rgu[num9 + i] -= num17;
}
if (uHi < num15)
{
uint uCarry = 0;
for (int j = 0; j < num; j++)
{
uCarry = AddCarry(ref regNum._rgu[num9 + j], regDen._rgu[j], uCarry);
}
num13 -= (ulong)1L;
}
regNum._iuLast = (num9 + num) - 1;
}
if (fQuo)
{
if (cu == 1)
{
regQuo._uSmall = (uint)num13;
}
else
{
regQuo._rgu[num9] = (uint)num13;
}
}
}
regNum._iuLast = num - 1;
regNum.Trim();
}
}
}
private static void ModDivCore(ref BigIntegerBuilder regNum, ref BigIntegerBuilder regDen, bool fQuo, ref BigIntegerBuilder regQuo)
{
regQuo.Set(0);
if (regNum._iuLast < regDen._iuLast)
{
return;
}
int cuDen = regDen._iuLast + 1;
int cuDiff = regNum._iuLast - regDen._iuLast;
// Determine whether the result will have cuDiff "digits" or cuDiff+1 "digits".
int cuQuo = cuDiff;
for (int iu = regNum._iuLast; ; iu--)
{
if (iu < cuDiff)
{
cuQuo++;
break;
}
if (regDen._rgu[iu - cuDiff] != regNum._rgu[iu])
{
if (regDen._rgu[iu - cuDiff] < regNum._rgu[iu])
{
cuQuo++;
}
break;
}
}
if (cuQuo == 0)
{
return;
}
if (fQuo)
{
regQuo.SetSizeLazy(cuQuo);
}
// Get the uint to use for the trial divisions. We normalize so the high bit is set.
uint uDen = regDen._rgu[cuDen - 1];
uint uDenNext = regDen._rgu[cuDen - 2];
int cbitShiftLeft = NumericsHelpers.CbitHighZero(uDen);
int cbitShiftRight = kcbitUint - cbitShiftLeft;
if (cbitShiftLeft > 0)
{
uDen = (uDen << cbitShiftLeft) | (uDenNext >> cbitShiftRight);
uDenNext <<= cbitShiftLeft;
if (cuDen > 2)
{
uDenNext |= regDen._rgu[cuDen - 3] >> cbitShiftRight;
}
}
// Allocate and initialize working space.
regNum.EnsureWritable();
for (int iu = cuQuo; --iu >= 0;)
{
// Get the high (normalized) bits of regNum.
uint uNumHi = (iu + cuDen <= regNum._iuLast) ? regNum._rgu[iu + cuDen] : 0;
ulong uuNum = NumericsHelpers.MakeUlong(uNumHi, regNum._rgu[iu + cuDen - 1]);
uint uNumNext = regNum._rgu[iu + cuDen - 2];
if (cbitShiftLeft > 0)
{
uuNum = (uuNum << cbitShiftLeft) | (uNumNext >> cbitShiftRight);
uNumNext <<= cbitShiftLeft;
if (iu + cuDen >= 3)
{
uNumNext |= regNum._rgu[iu + cuDen - 3] >> cbitShiftRight;
}
}
// Divide to get the quotient digit.
ulong uuQuo = uuNum / uDen;
ulong uuRem = (uint)(uuNum % uDen);
if (uuQuo > uint.MaxValue)
{
uuRem += uDen * (uuQuo - uint.MaxValue);
uuQuo = uint.MaxValue;
}
while (uuRem <= uint.MaxValue && uuQuo * uDenNext > NumericsHelpers.MakeUlong((uint)uuRem, uNumNext))
{
uuQuo--;
uuRem += uDen;
}
// Multiply and subtract. Note that uuQuo may be 1 too large. If we have a borrow
// at the end, we'll add the denominator back on and decrement uuQuo.
if (uuQuo > 0)
{
ulong uuBorrow = 0;
for (int iu2 = 0; iu2 < cuDen; iu2++)
{
uuBorrow += regDen._rgu[iu2] * uuQuo;
uint uSub = (uint)uuBorrow;
uuBorrow >>= kcbitUint;
if (regNum._rgu[iu + iu2] < uSub)
{
uuBorrow++;
}
regNum._rgu[iu + iu2] -= uSub;
}
if (uNumHi < uuBorrow)
{
// Add, tracking carry.
uint uCarry = 0;
for (int iu2 = 0; iu2 < cuDen; iu2++)
{
uCarry = AddCarry(ref regNum._rgu[iu + iu2], regDen._rgu[iu2], uCarry);
}
uuQuo--;
}
regNum._iuLast = iu + cuDen - 1;
}
if (fQuo)
{
if (cuQuo == 1)
{
regQuo._uSmall = (uint)uuQuo;
}
else
{
regQuo._rgu[iu] = (uint)uuQuo;
}
}
}
regNum._iuLast = cuDen - 1;
regNum.Trim();
}
private static void ModDivCore(ref BigIntegerBuilder regNum, ref BigIntegerBuilder regDen, bool fQuo, ref BigIntegerBuilder regQuo) {
regQuo.Set((uint)0);
if (regNum._iuLast >= regDen._iuLast) {
int num = regDen._iuLast + 1;
int num2 = regNum._iuLast - regDen._iuLast;
int cu = num2;
int index = regNum._iuLast;
while (true) {
if (index < num2) {
cu++;
break;
}
if (regDen._rgu[index - num2] != regNum._rgu[index]) {
if (regDen._rgu[index - num2] < regNum._rgu[index]) {
cu++;
}
break;
}
index--;
}
if (cu != 0) {
if (fQuo) {
regQuo.SetSizeLazy(cu);
}
uint u = regDen._rgu[num - 1];
uint num6 = regDen._rgu[num - 2];
int num7 = NumericsHelpers.CbitHighZero(u);
int num8 = 0x20 - num7;
if (num7 > 0) {
u = (u << num7) | (num6 >> num8);
num6 = num6 << num7;
if (num > 2) {
num6 |= regDen._rgu[num - 3] >> num8;
}
}
regNum.EnsureWritable();
int num9 = cu;
while (--num9 >= 0) {
uint uHi = ((num9 + num) <= regNum._iuLast) ? regNum._rgu[num9 + num] : 0;
ulong num11 = NumericsHelpers.MakeUlong(uHi, regNum._rgu[(num9 + num) - 1]);
uint uLo = regNum._rgu[(num9 + num) - 2];
if (num7 > 0) {
num11 = (num11 << num7) | (uLo >> num8);
uLo = uLo << num7;
if ((num9 + num) >= 3) {
uLo |= regNum._rgu[(num9 + num) - 3] >> num8;
}
}
ulong num13 = num11 / ((ulong)u);
ulong num14 = (uint)(num11 % ((ulong)u));
if (num13 > 0xffffffffL) {
num14 += (ulong)u * (num13 - 0xffffffffL);
num13 = 0xffffffffL;
}
while ((num14 <= 0xffffffffL) && ((num13 * num6) > NumericsHelpers.MakeUlong((uint)num14, uLo))) {
num13 -= (ulong)1L;
num14 += (ulong)u;
}
if (num13 > 0L) {
ulong num15 = 0L;
for (int i = 0; i < num; i++) {
num15 += regDen._rgu[i] * num13;
uint num17 = (uint)num15;
num15 = num15 >> 0x20;
if (regNum._rgu[num9 + i] < num17) {
num15 += (ulong)1L;
}
regNum._rgu[num9 + i] -= num17;
}
if (uHi < num15) {
uint uCarry = 0;
for (int j = 0; j < num; j++) {
uCarry = AddCarry(ref regNum._rgu[num9 + j], regDen._rgu[j], uCarry);
}
num13 -= (ulong)1L;
}
regNum._iuLast = (num9 + num) - 1;
}
if (fQuo) {
if (cu == 1) {
regQuo._uSmall = (uint)num13;
}
else {
regQuo._rgu[num9] = (uint)num13;
}
}
}
regNum._iuLast = num - 1;
regNum.Trim();
}
}
}
private static void ModDivCore(ref BigIntegerBuilder regNum, ref BigIntegerBuilder regDen, bool fQuo, ref BigIntegerBuilder regQuo)
{
regQuo.Set(0);
if (regNum._iuLast < regDen._iuLast)
{
return;
}
var num1 = regDen._iuLast + 1;
var num2 = regNum._iuLast - regDen._iuLast;
var num3 = num2;
var num4 = regNum._iuLast;
while (true)
{
if (num4 < num2)
{
num3++;
break;
}
if (regDen._bits[num4 - num2] == regNum._bits[num4])
{
num4--;
}
else
{
if (regDen._bits[num4 - num2] < regNum._bits[num4])
{
num3++;
}
break;
}
}
if (num3 == 0)
{
return;
}
if (fQuo)
{
regQuo.SetSizeLazy(num3);
}
var num5 = regDen._bits[num1 - 1];
var num6 = regDen._bits[num1 - 2];
var num7 = NumericHelper.CbitHighZero(num5);
var num8 = 32 - num7;
if (num7 > 0)
{
num5 = (num5 << (num7 & 31)) | (num6 >> (num8 & 31));
num6 <<= num7 & 31;
if (num1 > 2)
{
num6 |= regDen._bits[num1 - 3] >> (num8 & 31);
}
}
regNum.EnsureWritable();
var num9 = num3;
while (true)
{
var num10 = num9 - 1;
num9 = num10;
if (num10 < 0)
{
break;
}
var num = num9 + num1 > regNum._iuLast ? 0 : regNum._bits[num9 + num1];
var num11 = num;
var num12 = NumericHelper.BuildUInt64(num11, regNum._bits[num9 + num1 - 1]);
var num13 = regNum._bits[num9 + num1 - 2];
if (num7 > 0)
{
num12 = (num12 << (num7 & 63)) | (num13 >> (num8 & 31));
num13 <<= num7 & 31;
if (num9 + num1 >= 3)
{
num13 |= regNum._bits[num9 + num1 - 3] >> (num8 & 31);
}
}
var num14 = num12 / num5;
var num15 = (ulong)(uint)(num12 % num5);
if (num14 > uint.MaxValue)
{
num15 += num5 * (num14 - uint.MaxValue);
num14 = uint.MaxValue;
}
while (num15 <= uint.MaxValue && num14 * num6 > NumericHelper.BuildUInt64((uint)num15, num13))
{
num14--;
num15 += num5;
}
if (num14 > 0)
{
var num16 = (ulong)0;
for (var i = 0; i < num1; i++)
{
num16 += regDen._bits[i] * num14;
var num17 = (uint)num16;
num16 >>= 32;
if (regNum._bits[num9 + i] < num17)
{
num16++;
}
regNum._bits[num9 + i] -= num17;
}
if (num11 < num16)
{
uint num18 = 0;
for (var j = 0; j < num1; j++)
{
num18 = AddCarry(ref regNum._bits[num9 + j], regDen._bits[j], num18);
}
num14--;
}
regNum._iuLast = num9 + num1 - 1;
}
if (!fQuo)
{
continue;
}
if (num3 != 1)
{
regQuo._bits[num9] = (uint)num14;
}
else
{
regQuo._uSmall = (uint)num14;
}
}
regNum._iuLast = num1 - 1;
regNum.Trim();
}
private static void ModDivCore(ref BigIntegerBuilder regNum, ref BigIntegerBuilder regDen, bool fQuo, ref BigIntegerBuilder regQuo) {
Contract.Assert(regNum._iuLast > 0 && regDen._iuLast > 0);
regQuo.Set(0);
if (regNum._iuLast < regDen._iuLast)
return;
Contract.Assert(0 < regDen._iuLast && regDen._iuLast <= regNum._iuLast);
int cuDen = regDen._iuLast + 1;
int cuDiff = regNum._iuLast - regDen._iuLast;
// Determine whether the result will have cuDiff "digits" or cuDiff+1 "digits".
int cuQuo = cuDiff;
for (int iu = regNum._iuLast; ; iu--) {
if (iu < cuDiff) {
cuQuo++;
break;
}
if (regDen._rgu[iu - cuDiff] != regNum._rgu[iu]) {
if (regDen._rgu[iu - cuDiff] < regNum._rgu[iu])
cuQuo++;
break;
}
}
if (cuQuo == 0)
return;
if (fQuo)
regQuo.SetSizeLazy(cuQuo);
// Get the uint to use for the trial divisions. We normalize so the high bit is set.
uint uDen = regDen._rgu[cuDen - 1];
uint uDenNext = regDen._rgu[cuDen - 2];
int cbitShiftLeft = NumericsHelpers.CbitHighZero(uDen);
int cbitShiftRight = kcbitUint - cbitShiftLeft;
if (cbitShiftLeft > 0) {
uDen = (uDen << cbitShiftLeft) | (uDenNext >> cbitShiftRight);
uDenNext <<= cbitShiftLeft;
if (cuDen > 2)
uDenNext |= regDen._rgu[cuDen - 3] >> cbitShiftRight;
}
Contract.Assert((uDen & 0x80000000) != 0);
// Allocate and initialize working space.
Contract.Assert(cuQuo + cuDen == regNum._iuLast + 1 || cuQuo + cuDen == regNum._iuLast + 2);
regNum.EnsureWritable();
for (int iu = cuQuo; --iu >= 0; ) {
// Get the high (normalized) bits of regNum.
uint uNumHi = (iu + cuDen <= regNum._iuLast) ? regNum._rgu[iu + cuDen] : 0;
Contract.Assert(uNumHi <= regDen._rgu[cuDen - 1]);
ulong uuNum = NumericsHelpers.MakeUlong(uNumHi, regNum._rgu[iu + cuDen - 1]);
uint uNumNext = regNum._rgu[iu + cuDen - 2];
if (cbitShiftLeft > 0) {
uuNum = (uuNum << cbitShiftLeft) | (uNumNext >> cbitShiftRight);
uNumNext <<= cbitShiftLeft;
if (iu + cuDen >= 3)
uNumNext |= regNum._rgu[iu + cuDen - 3] >> cbitShiftRight;
}
// Divide to get the quotient digit.
ulong uuQuo = uuNum / uDen;
ulong uuRem = (uint)(uuNum % uDen);
Contract.Assert(uuQuo <= (ulong)uint.MaxValue + 2);
if (uuQuo > uint.MaxValue) {
uuRem += uDen * (uuQuo - uint.MaxValue);
uuQuo = uint.MaxValue;
}
while (uuRem <= uint.MaxValue && uuQuo * uDenNext > NumericsHelpers.MakeUlong((uint)uuRem, uNumNext)) {
uuQuo--;
uuRem += uDen;
}
// Multiply and subtract. Note that uuQuo may be 1 too large. If we have a borrow
// at the end, we'll add the denominator back on and decrement uuQuo.
if (uuQuo > 0) {
ulong uuBorrow = 0;
for (int iu2 = 0; iu2 < cuDen; iu2++) {
uuBorrow += regDen._rgu[iu2] * uuQuo;
uint uSub = (uint)uuBorrow;
uuBorrow >>= kcbitUint;
if (regNum._rgu[iu + iu2] < uSub)
uuBorrow++;
regNum._rgu[iu + iu2] -= uSub;
}
Contract.Assert(uNumHi == uuBorrow || uNumHi == uuBorrow - 1);
if (uNumHi < uuBorrow) {
// Add, tracking carry.
uint uCarry = 0;
for (int iu2 = 0; iu2 < cuDen; iu2++) {
uCarry = AddCarry(ref regNum._rgu[iu + iu2], regDen._rgu[iu2], uCarry);
Contract.Assert(uCarry <= 1);
}
Contract.Assert(uCarry == 1);
uuQuo--;
}
regNum._iuLast = iu + cuDen - 1;
}
if (fQuo) {
if (cuQuo == 1)
regQuo._uSmall = (uint)uuQuo;
else
regQuo._rgu[iu] = (uint)uuQuo;
}
}
Contract.Assert(cuDen > 1 && regNum._iuLast > 0);
regNum._iuLast = cuDen - 1;
regNum.Trim();
}
