1 #include <cln/integer.h>
2 #include <cln/integer_io.h>
10 std::cerr << "Assertion failed! File " << __FILE__ << ", line " << __LINE__ << endl; \
20 #define num_elements(array) (sizeof(array)/sizeof(array[0]))
22 // Note: This macro differs slightly from the one in "exam.h".
23 #define DO_BINOP_TEST(typename,type,rtype,opname) \
24 static int test_##typename##_##opname (void) \
27 for (unsigned int i = 0; i < num_elements(typename##_##opname##_tests); i++) { \
28 opname##_test& test = typename##_##opname##_tests[i]; \
29 type arg1 = type(test.arg1); \
30 type arg2 = type(test.arg2); \
31 rtype computed_result = opname(arg1,arg2); \
32 rtype result = rtype(test.result); \
33 if (computed_result != result) { \
34 std::cerr << "Error in " #typename "_" #opname "_tests[" << i << "] !" << endl; \
35 std::cerr << "Result should be: " << result << endl; \
36 std::cerr << "Result computed : " << computed_result << endl << endl; \
43 static gcd_test integer_gcd_tests[] = {
44 { "123456789", "345", "3" },
45 { "345", "123456789", "3" },
48 { "2523533737", "855322739", "1" },
49 { "855322739", "2523533737", "1" },
50 { "101611479673163974026724715741235467160607959655653420075620", "533177863832047932237026621580126811198495699416238676294977", "1" },
51 { "30729415811", "323233683197", "31071199" },
52 { "77874422", "32223899", "1" },
53 { "974507656412513757857315037382926980395082974811562770185617915360", "-1539496810360685510909469177732386446833404488164283", "1" },
54 { "2823618260546496405819033080103700734250203999069672146446", "18374686479688400895", "1" }
57 DO_BINOP_TEST(integer,cl_I,cl_I,gcd)
62 error |= test_integer_gcd();
74 cl_I g = xgcd(a,b, &u,&v);
77 ASSERT(u == -9206830);
78 ASSERT(v == 22249839);
85 cl_I g = xgcd(a,b, &u,&v);
88 ASSERT(u == 77165803);
89 ASSERT(v == -138134388);
92 cl_I a = "#x80000000";
93 cl_I b = "#x-C0000000";
96 cl_I g = xgcd(a,b, &u,&v);
97 ASSERT(g == (cl_I)"#x40000000");
103 cl_I a = "974507656412513757857315037382926980395082974811562770185617915360";
104 cl_I b = "-1539496810360685510909469177732386446833404488164283";
107 cl_I g = xgcd(a,b, &u,&v);
109 ASSERT(g == a*u+b*v);