void testDivRand512Bit() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
//512 bits
BigInt num("22479199231365811817124860008034"
"86229160630675450396401816569207"
"37104752501976414843519588079979"
"01710354041585586289392397420347"
"216513925045972181584884921");
//500 bits
BigInt denom("351508594074011603394279375147"
"626833956407011572250308540816"
"854541429976970595757637034702"
"593071101429245293222692961211"
"2896690197169368617519259819528");
start = std::chrono::system_clock::now();
end = std::chrono::system_clock::now();
BigInt result = num / denom;
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "DivRand512Bit took: " << elapsed_time.count() << " computing " << result << std::endl;
#endif
std::vector<limb_t> actual{6395};
std::cout << "512b num / 500b num Correct? " << testEquals(result, actual) << std::endl;
}
void test(){
testAdd();
testSub();
testMult();
testScale();
testEquals();
testDivision();
testAbsolute();
testArgument();
printf("All passed!\n");
}
void test60Fact() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
start = std::chrono::system_clock::now();
BigInt fact60 = Factorial(60);
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "fact60 took: " << elapsed_time.count() << " computing " << fact60 << std::endl;
#endif
std::vector<limb_t> actual{18396530,1065373609,393476149,835766904,1949166040,1155626993,1901353954,234881024,0};
std::cout << "60! Correct? " << testEquals(fact60, actual) << std::endl;
}
void testLongStringConstructor() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
start = std::chrono::system_clock::now();
BigInt longString("12345678998765432101234567899876543210");
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "longString took: " << elapsed_time.count() << " computing " << longString << std::endl;
#endif
std::vector<limb_t> actual{1246595009,748116907,1035237483,1632488170};
std::cout << "String Constructor Correct? " << testEquals(longString, actual) << std::endl;
}
void test500Fib() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
start = std::chrono::system_clock::now();
BigInt fib500 = Fibonacci(500);
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "fib500 took: " << elapsed_time.count() << " computing " << fib500 << std::endl;
#endif
std::vector<limb_t> actual{31, 268010877,1041980198,1505931748,939335682,475289477,
384828411,411953712,1000216971,1090861368,1011151204,315178285};
std::cout << "Fib(500) Correct? " << testEquals(fib500, actual) << std::endl;
}
void test2Pow1024(){
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
start = std::chrono::system_clock::now();
BigInt pow(2);
pow = pow.pow(1024);
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "2Pow1024 took: " << elapsed_time.count() << " computing " << pow << std::endl;
#endif
std::vector<limb_t> actual{2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
std::cout << "2^1024 Correct? " << testEquals(pow, actual) << std::endl;
}
void test20FactQuad() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
start = std::chrono::system_clock::now();
BigInt fact = Factorial(20);
BigInt factSquared = fact * fact;
factSquared *= factSquared;
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "fact20quad took: " << elapsed_time.count() << " computing " << factSquared << std::endl;
#endif
std::vector<limb_t> actual{166337208,764515097,443847783,29956932,1554513582,703611904,0,0};
std::cout << "(20!)^4 Correct? " << testEquals(factSquared, actual) << std::endl;
}
void testRemainderlessMod() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
BigInt num("10000000000000000000000000000000000");
BigInt modulus("10000000000000000");
start = std::chrono::system_clock::now();
BigInt q = num % modulus;
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "testRemainderlessMod took: " << elapsed_time.count() << " computing " << q << std::endl;
#endif
std::vector<limb_t> actual{0};
std::cout << "Remainderless mod Correct? " << testEquals(q, actual) << std::endl;
}
void testDivAddBack() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
BigInt num("12345678998765432112345678901111111111222222222");
BigInt denom("12345678998765431212340871");
start = std::chrono::system_clock::now();
BigInt q = num / denom;
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "testDivAddBack took: " << elapsed_time.count() << " computing " << q << std::endl;
#endif
std::vector<limb_t> actual{216,1804819339,1587616964};
std::cout << "Division\'s 'add back' case Correct? " << testEquals(q, actual) << std::endl;
}
void test2Pow31Pow2Pow10() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
BigInt twoPow31(2);
twoPow31 = twoPow31.pow(31);
start = std::chrono::system_clock::now();
twoPow31 = twoPow31.pow(1024);
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "2Pow31Pow2Pow10 took: " << elapsed_time.count() << " computing " << twoPow31 << std::endl;
#endif
std::vector<limb_t> actual{ 1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
std::cout << "(2^31)^(2^10) Correct? " << testEquals(twoPow31, actual) << std::endl;
}
void testSmallModExp() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
BigInt base(3);
BigInt exp(6);
BigInt mod(7);
start = std::chrono::system_clock::now();
BigInt result = base.pow(exp, mod);
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "testSmallExp took: " << elapsed_time.count() << " computing " << result << std::endl;
#endif
std::vector<limb_t> actual{1};
std::cout << "3^6 mod 7 Correct? " << testEquals(result, actual) << std::endl;
}
void test10Pow50() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
start = std::chrono::system_clock::now();
BigInt pow(10);
pow = pow.pow(50);
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "10Pow50 took: " << elapsed_time.count() << " computing " << pow << std::endl;
#endif
std::vector<limb_t> actual{2189,1135697817,322190998,269800438,1791492096,0};
std::cout << "10^50 Correct? " << testEquals(pow, actual) << std::endl;
}
void testVeryLongStringConstructor() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
start = std::chrono::system_clock::now();
BigInt longString(
#include "bignum.txt"
);
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "VeryLongString took: " << elapsed_time.count() << " computing " << longString << std::endl;
#endif
std::vector<limb_t> actual{ 1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
std::cout << "Very-long String Constructor Correct? " << testEquals(longString, actual) << std::endl;
}
void test3Pow512Div2Pow512() {
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
BigInt pow2(2);
pow2 = pow2.pow(512);
BigInt pow3(3);
pow3 = pow3.pow(512);
start = std::chrono::system_clock::now();
BigInt pow = pow3 / pow2;
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "3Pow512Div2Pow512 took: " << elapsed_time.count() << " computing " << pow << std::endl;
#endif
std::vector<limb_t> actual{1483733,662018629,1236000823,1498554168,837056997,892058291,616689941,1251421047,1182766622,442154954};
std::cout << "3^512 / 2^512 Correct? " << testEquals(pow, actual) << std::endl;
}
void testMulRandom1024Bit() {
BigInt a("276601373654928975871754674997064595262751078392398403"
"140078759572403299733275944955495590985079635677358326"
"437768288362804324337725761906816178599761686458110239"
"720130379865178259907677503176456269594835847095657438"
"537293270309725637889873193915190724261885472342077067"
"844166147815436649001795340162949998568");
BigInt b("908270564241715128928475652095664817489270259211612103"
"998361943067766201208031918379136862415371024909615030"
"639734030078204866227763421426612690401640345569552459"
"305270873069670449962612360682955706072635452841457707"
"721202565127178715538792506240083189193164965778468706"
"43406500278810140967997440157799795150");
std::chrono::time_point<std::chrono::system_clock> start, end;
std::chrono::duration<double> elapsed_time;
start = std::chrono::system_clock::now();
a *= b;
end = std::chrono::system_clock::now();
elapsed_time = end - start;
#ifdef _PRINT_VALS
std::cout<< "MulRandom1024Bit took: " << elapsed_time.count() << " computing " << a.ToDecimal() << std::endl;
#endif
// std::cout << "Avg time: " << elapsed_time.count() / 100000 << std::endl;
std::vector<limb_t> actual{3,235270659,568911404,922708634,1696526154,1010925488,110054269,
1158468645,1442079807,1577453968,299033823,519200492,1286835655,
1419899891,1067667472,1499562046,108550203,1558147072,1198215629,
78649644,329564816,225153618,1824298395,790487858,1994024926,1000868029,
1812712781,1156872405,1574232354,425287274,962686480,152563639,234031561,
964340897,896011126,672551734,473018745,1981102438,768108938,571079643,
276511614,1067268780,904046332,1878006326,943061062,675882021,1428995088,
41648230,344962290,897227716,95647102,355939536,379735193,346761280,
548236789,445511246,409511175,1229211920,1463486602,1827214088,1092666381,
1236026364,628671149,562427339,201968067,1990132260,823182512,};
std::cout << "Mul random 1024bit numbers Correct? " << testEquals(a, actual) << std::endl;
}
