1 // Test quadratic algorithm.
11 // I used bc to get the reference values.
14 cf_t pi
= cf_new_pi();
16 for (int i
= 0; i
< 8; i
++) {
17 mpz_init(addarray
[i
]);
19 mpz_set_ui(addarray
[1], 1);
20 mpz_set_ui(addarray
[2], 1);
21 mpz_set_ui(addarray
[7], 1);
23 cf_t b
= cf_new_bihom(e
, pi
, addarray
);
25 CF_EXPECT_DEC(b
, "58598744820488384738");
31 // Check 2 sin 1 cos 2 = sin 2 = 0.909...
35 mpz_set_ui(addarray
[0], 2);
36 mpz_set_ui(addarray
[1], 0);
37 mpz_set_ui(addarray
[2], 0);
38 mpz_set_si(addarray
[7], 1);
39 b
= cf_new_bihom(s1
, c1
, addarray
);
41 CF_EXPECT_DEC(b
, "09092974268256816953");
46 for (int i
= 0; i
< 8; i
++) {
47 mpz_clear(addarray
[i
]);