int main() { long iter; flint_rand_t state; printf("root...."); fflush(stdout); flint_randinit(state); for (iter = 0; iter < 100000; iter++) { fmprb_t a, b, c; ulong k; long prec; prec = 2 + n_randint(state, 200); k = 1 + n_randint(state, 10); fmprb_init(a); fmprb_init(b); fmprb_init(c); fmprb_randtest(a, state, 1 + n_randint(state, 200), 10); fmprb_root(b, a, k, prec); fmprb_pow_ui(c, b, k, prec); if (!fmprb_contains(c, a)) { printf("FAIL: containment\n\n"); printf("k = %lu\n", k); printf("a = "); fmprb_print(a); printf("\n\n"); printf("b = "); fmprb_print(b); printf("\n\n"); printf("c = "); fmprb_print(c); printf("\n\n"); abort(); } fmprb_root(a, a, k, prec); if (!fmprb_equal(a, b)) { printf("FAIL: aliasing\n\n"); abort(); } fmprb_clear(a); fmprb_clear(b); fmprb_clear(c); } flint_randclear(state); flint_cleanup(); printf("PASS\n"); return EXIT_SUCCESS; }
int main() { slong iter; flint_rand_t state; flint_printf("rsqrt...."); fflush(stdout); flint_randinit(state); for (iter = 0; iter < 100000; iter++) { fmprb_t a, b, c; slong prec = 2 + n_randint(state, 200); fmprb_init(a); fmprb_init(b); fmprb_init(c); fmprb_randtest(a, state, 1 + n_randint(state, 200), 10); fmprb_rsqrt(b, a, prec); fmprb_inv(c, b, prec); fmprb_mul(c, c, c, prec); if (!fmprb_contains(c, a)) { flint_printf("FAIL: containment\n\n"); flint_printf("a = "); fmprb_print(a); flint_printf("\n\n"); flint_printf("b = "); fmprb_print(b); flint_printf("\n\n"); flint_printf("c = "); fmprb_print(c); flint_printf("\n\n"); abort(); } fmprb_rsqrt(a, a, prec); if (!fmprb_equal(a, b)) { flint_printf("FAIL: aliasing\n\n"); abort(); } fmprb_clear(a); fmprb_clear(b); fmprb_clear(c); } flint_randclear(state); flint_cleanup(); flint_printf("PASS\n"); return EXIT_SUCCESS; }
int fmprb_poly_contains(const fmprb_poly_t poly1, const fmprb_poly_t poly2) { long i; if (poly2->length > poly1->length) return 0; for (i = 0; i < poly2->length; i++) { if (!fmprb_contains(poly1->coeffs + i, poly2->coeffs + i)) return 0; } for (i = poly2->length; i < poly1->length; i++) if (!fmprb_contains_zero(poly1->coeffs + i)) return 0; return 1; }
int main() { slong iter; flint_rand_t state; flint_printf("union...."); fflush(stdout); flint_randinit(state); for (iter = 0; iter < 100000; iter++) { fmprb_t x, y, z; slong prec; int alias; fmprb_init(x); fmprb_init(y); fmprb_init(z); fmprb_randtest_special(x, state, 200, 10); fmprb_randtest_special(y, state, 200, 10); fmprb_randtest_special(z, state, 200, 10); prec = 2 + n_randint(state, 200); fmprb_union(z, x, y, prec); if (!fmprb_contains(z, x) || !fmprb_contains(z, y)) { flint_printf("FAIL:\n\n"); flint_printf("x = "); fmprb_print(x); flint_printf("\n\n"); flint_printf("y = "); fmprb_print(y); flint_printf("\n\n"); flint_printf("z = "); fmprb_print(z); flint_printf("\n\n"); abort(); } if (n_randint(state, 2)) { fmprb_union(x, x, y, prec); alias = fmprb_equal(x, z); } else { fmprb_union(y, x, y, prec); alias = fmprb_equal(y, z); } if (!alias) { flint_printf("FAIL (aliasing):\n\n"); flint_printf("x = "); fmprb_print(x); flint_printf("\n\n"); flint_printf("y = "); fmprb_print(y); flint_printf("\n\n"); flint_printf("z = "); fmprb_print(z); flint_printf("\n\n"); abort(); } fmprb_clear(x); fmprb_clear(y); fmprb_clear(z); } flint_randclear(state); flint_cleanup(); flint_printf("PASS\n"); return EXIT_SUCCESS; }
int main() { long iter; flint_rand_t state; printf("acos...."); fflush(stdout); flint_randinit(state); for (iter = 0; iter < 10000; iter++) { fmprb_t a, b; fmpq_t q; mpfr_t t; long prec = 2 + n_randint(state, 200); fmprb_init(a); fmprb_init(b); fmpq_init(q); mpfr_init2(t, prec + 100); fmprb_randtest(a, state, 1 + n_randint(state, 200), 3); fmprb_randtest(b, state, 1 + n_randint(state, 200), 3); fmprb_get_rand_fmpq(q, state, a, 1 + n_randint(state, 200)); fmpq_get_mpfr(t, q, MPFR_RNDN); mpfr_acos(t, t, MPFR_RNDN); fmprb_acos(b, a, prec); if (!fmprb_contains_mpfr(b, t)) { printf("FAIL: containment\n\n"); printf("a = "); fmprb_printd(a, 100); printf("\n\n"); printf("b = "); fmprb_printd(b, 100); printf("\n\n"); abort(); } fmprb_acos(a, a, prec); if (!fmprb_equal(a, b)) { printf("FAIL: aliasing\n\n"); abort(); } fmprb_clear(a); fmprb_clear(b); fmpq_clear(q); mpfr_clear(t); } /* check large arguments */ for (iter = 0; iter < 10000; iter++) { fmprb_t a, b, c; long prec1, prec2; prec1 = 2 + n_randint(state, 1000); prec2 = prec1 + 30; fmprb_init(a); fmprb_init(b); fmprb_init(c); fmprb_randtest_precise(a, state, 1 + n_randint(state, 1000), 100); fmprb_acos(b, a, prec1); fmprb_acos(c, a, prec2); if (!fmprb_overlaps(b, c)) { printf("FAIL: overlap\n\n"); printf("a = "); fmprb_print(a); printf("\n\n"); printf("b = "); fmprb_print(b); printf("\n\n"); printf("c = "); fmprb_print(c); printf("\n\n"); abort(); } /* check sin(asin(x)) = x */ fmprb_cos(c, b, prec1); if (!fmprb_contains(c, a)) { printf("FAIL: functional equation\n\n"); printf("a = "); fmprb_print(a); printf("\n\n"); printf("b = "); fmprb_print(b); printf("\n\n"); printf("c = "); fmprb_print(c); printf("\n\n"); abort(); } fmprb_clear(a); fmprb_clear(b); fmprb_clear(c); } flint_randclear(state); flint_cleanup(); printf("PASS\n"); return EXIT_SUCCESS; }
int main() { slong iter; flint_rand_t state; flint_printf("agm...."); fflush(stdout); flint_randinit(state); for (iter = 0; iter < 100000; iter++) { fmprb_t a, b, c; fmpq_t q, r; mpfr_t t, u; slong prec = 2 + n_randint(state, 200); fmprb_init(a); fmprb_init(b); fmprb_init(c); fmpq_init(q); fmpq_init(r); mpfr_init2(t, prec + 100); mpfr_init2(u, prec + 100); fmprb_randtest(a, state, 1 + n_randint(state, 200), 3); fmprb_randtest(b, state, 1 + n_randint(state, 200), 3); fmprb_randtest(c, state, 1 + n_randint(state, 200), 3); fmprb_agm(c, a, b, prec); if (fmprb_equal(a, b)) { if (!fmprb_contains(c, a)) { flint_printf("FAIL: containment (identity)\n\n"); flint_printf("a = "); fmprb_print(a); flint_printf("\n\n"); flint_printf("b = "); fmprb_print(b); flint_printf("\n\n"); flint_printf("c = "); fmprb_print(c); flint_printf("\n\n"); abort(); } } else { fmprb_get_rand_fmpq(q, state, a, 1 + n_randint(state, 200)); fmprb_get_rand_fmpq(r, state, b, 1 + n_randint(state, 200)); fmpq_get_mpfr(t, q, MPFR_RNDN); fmpq_get_mpfr(u, r, MPFR_RNDN); mpfr_agm(t, t, u, MPFR_RNDN); if (!fmprb_contains_mpfr(c, t)) { flint_printf("FAIL: containment\n\n"); flint_printf("a = "); fmprb_print(a); flint_printf("\n\n"); flint_printf("b = "); fmprb_print(b); flint_printf("\n\n"); flint_printf("c = "); fmprb_print(c); flint_printf("\n\n"); abort(); } } fmprb_clear(a); fmprb_clear(b); fmprb_clear(c); fmpq_clear(q); fmpq_clear(r); mpfr_clear(t); mpfr_clear(u); } flint_randclear(state); flint_cleanup(); flint_printf("PASS\n"); return EXIT_SUCCESS; }