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
, "5.8598744820488384738");
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
, "0.9092974268256816953");
47 // Check 2 (cos 1)^2 - 1 = cos 2 =
48 s1
= cf_new_cos1(); // TODO: Implement tee, use that instead.
50 mpz_set_si(addarray
[0], 2);
51 mpz_set_si(addarray
[3], -1);
52 mpz_set_si(addarray
[7], 1);
53 b
= cf_new_bihom(s1
, c1
, addarray
);
55 CF_EXPECT_DEC(b
, "-0.41614683654714238699");
61 for (int i
= 0; i
< 8; i
++) {
62 mpz_clear(addarray
[i
]);