void test_scalar_div_si(char * in, slong x, char * out) { int ans; fmpz_poly_q_t rop, op; char * res; fmpz_poly_q_init(op); fmpz_poly_q_set_str(op, in); fmpz_poly_q_init(rop); fmpz_poly_q_scalar_div_si(rop, op, x); res = fmpz_poly_q_get_str(rop); ans = !strcmp(out, res); if (!ans) { flint_printf("test_scalar_div_si: failed\n"); flint_printf(" Expected \"%s\", got \"%s\"\n", out, res); abort(); } fmpz_poly_q_clear(op); fmpz_poly_q_clear(rop); flint_free(res); }
int main(void) { int i, result; flint_rand_t state; printf("scalar_div_si... "); fflush(stdout); flint_randinit(state); /* Check aliasing of a and b */ for (i = 0; i < 100; i++) { fmpz_poly_q_t a, b; long x; fmpz_poly_q_init(a); fmpz_poly_q_init(b); fmpz_poly_q_randtest(b, state, n_randint(state, 50), 50, n_randint(state, 50), 50); x = z_randtest_not_zero(state); fmpz_poly_q_scalar_div_si(a, b, x); fmpz_poly_q_scalar_div_si(b, b, x); result = fmpz_poly_q_equal(a, b); if (!result) { printf("FAIL:\n"); fmpz_poly_q_print(a), printf("\n\n"); fmpz_poly_q_print(b), printf("\n\n"); abort(); } fmpz_poly_q_clear(a); fmpz_poly_q_clear(b); } /* Check that (a + b) / x == a / x + b / x */ for (i = 0; i < 100; i++) { fmpz_poly_q_t a, b, c, d; long x; fmpz_poly_q_init(a); fmpz_poly_q_init(b); fmpz_poly_q_init(c); fmpz_poly_q_init(d); fmpz_poly_q_randtest(a, state, n_randint(state, 50), 50, n_randint(state, 50), 50); fmpz_poly_q_randtest(b, state, n_randint(state, 50), 50, n_randint(state, 50), 50); x = z_randtest_not_zero(state); fmpz_poly_q_scalar_div_si(c, a, x); fmpz_poly_q_scalar_div_si(d, b, x); fmpz_poly_q_add(d, c, d); fmpz_poly_q_add(c, a, b); fmpz_poly_q_scalar_div_si(c, c, x); result = fmpz_poly_q_equal(c, d) && fmpz_poly_q_is_canonical(c); if (!result) { printf("FAIL:\n"); printf("a = "), fmpz_poly_q_print(a), printf("\n\n"); printf("b = "), fmpz_poly_q_print(b), printf("\n\n"); printf("c = "), fmpz_poly_q_print(c), printf("\n\n"); printf("d = "), fmpz_poly_q_print(d), printf("\n\n"); printf("x = %ld\n\n", x); abort(); } fmpz_poly_q_clear(a); fmpz_poly_q_clear(b); fmpz_poly_q_clear(c); fmpz_poly_q_clear(d); } flint_randclear(state); _fmpz_cleanup(); printf("PASS\n"); return EXIT_SUCCESS; }