void test_submul(char * in1, char * in2, char * in3, char * out) { int ans; fmpz_poly_q_t rop, op1, op2; char * res; fmpz_poly_q_init(rop); fmpz_poly_q_set_str(rop, in1); fmpz_poly_q_init(op1); fmpz_poly_q_set_str(op1, in2); fmpz_poly_q_init(op2); fmpz_poly_q_set_str(op2, in3); fmpz_poly_q_submul(rop, op1, op2); res = fmpz_poly_q_get_str(rop); ans = !strcmp(out, res); if (!ans) { flint_printf("test_submul: failed\n"); flint_printf(" Expected \"%s\", got \"%s\"\n", out, res); abort(); } fmpz_poly_q_clear(rop); fmpz_poly_q_clear(op1); fmpz_poly_q_clear(op2); flint_free(res); }
int main(void) { int i, result; flint_rand_t state; printf("submul... "); fflush(stdout); flint_randinit(state); /* Check aliasing of a and c */ for (i = 0; i < 100; i++) { fmpz_poly_q_t a, b, c; fmpz_poly_q_init(a); fmpz_poly_q_init(b); fmpz_poly_q_init(c); 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); fmpz_poly_q_set(c, a); fmpz_poly_q_submul(c, a, b); fmpz_poly_q_submul(a, a, b); result = fmpz_poly_q_equal(a, c) && fmpz_poly_q_is_canonical(a); if (!result) { printf("FAIL:\n"); fmpz_poly_q_print(a), printf("\n\n"); fmpz_poly_q_print(b), printf("\n\n"); fmpz_poly_q_print(c), printf("\n\n"); abort(); } fmpz_poly_q_clear(a); fmpz_poly_q_clear(b); fmpz_poly_q_clear(c); } /* Check aliasing of b and c */ for (i = 0; i < 100; i++) { fmpz_poly_q_t a, b, c; fmpz_poly_q_init(a); fmpz_poly_q_init(b); fmpz_poly_q_init(c); 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); fmpz_poly_q_set(c, b); fmpz_poly_q_submul(c, a, b); fmpz_poly_q_submul(b, a, b); result = fmpz_poly_q_equal(b, c); if (!result) { printf("FAIL:\n"); fmpz_poly_q_print(a), printf("\n\n"); fmpz_poly_q_print(b), printf("\n\n"); fmpz_poly_q_print(c), printf("\n\n"); abort(); } fmpz_poly_q_clear(a); fmpz_poly_q_clear(b); fmpz_poly_q_clear(c); } flint_randclear(state); _fmpz_cleanup(); printf("PASS\n"); return EXIT_SUCCESS; }