void divide(ElementType& result, const ElementType& a, const ElementType& b) const { // We need to handle the case when result is a, or b. // This is why we use the temporary value 'c'. ElementType c; init(c); #if 0 printf("entering divide\n"); printf(" a = "); fq_zech_print_pretty(&a, mContext); printf("\n b = "); fq_zech_print_pretty(&b, mContext); #endif assert(not is_zero(b)); invert(c, b); #if 0 printf("\n 1/b = "); fq_zech_print_pretty(&c, mContext); #endif mult(result, c, a); #if 0 printf("\n a/b = "); fq_zech_print_pretty(&result, mContext); printf("\n"); #endif clear(c); }
int main(void) { int j, i, result; fq_zech_ctx_t ctx; FLINT_TEST_INIT(state); flint_printf("mul_fmpz... "); fflush(stdout); for (j = 0; j < 50; j++) { fq_zech_ctx_randtest(ctx, state); for (i = 0; i < 200; i++) { fmpz_t x; fq_nmod_t aa, bb; fq_zech_t a, b, c; fq_nmod_init(aa, ctx->fq_nmod_ctx); fq_nmod_init(bb, ctx->fq_nmod_ctx); fmpz_init(x); fmpz_randtest_mod_signed(x, state, fq_zech_ctx_prime(ctx)); fq_nmod_randtest(aa, state, ctx->fq_nmod_ctx); fq_zech_set_fq_nmod(a, aa, ctx); fq_nmod_mul_fmpz(bb, aa, x, ctx->fq_nmod_ctx); fq_zech_set_fq_nmod(b, bb, ctx); fq_zech_mul_fmpz(c, a, x, ctx); result = (fq_zech_equal(b, c, ctx)); if (!result) { flint_printf("FAIL:\n\n"); fq_zech_ctx_print(ctx); flint_printf("\n"); flint_printf("aa = "); fq_nmod_print_pretty(aa, ctx->fq_nmod_ctx); flint_printf("\n"); flint_printf("a = "); fq_zech_print_pretty(a, ctx); flint_printf("\n"); flint_printf("b = "); fq_zech_print_pretty(b, ctx); flint_printf("\n"); flint_printf("c = "); fq_zech_print_pretty(c, ctx); flint_printf("\n"); abort(); } fmpz_clear(x); fq_nmod_clear(bb, ctx->fq_nmod_ctx); fq_nmod_clear(aa, ctx->fq_nmod_ctx); } for (i = 0; i < 200; i++) { fmpz_t x; fq_nmod_t aa, bb; fq_zech_t a, b; fq_nmod_init(aa, ctx->fq_nmod_ctx); fq_nmod_init(bb, ctx->fq_nmod_ctx); fmpz_init(x); fmpz_randtest_mod_signed(x, state, fq_zech_ctx_prime(ctx)); fq_nmod_randtest(aa, state, ctx->fq_nmod_ctx); fq_zech_set_fq_nmod(a, aa, ctx); fq_nmod_mul_fmpz(bb, aa, x, ctx->fq_nmod_ctx); fq_zech_set_fq_nmod(b, bb, ctx); fq_zech_mul_fmpz(a, a, x, ctx); result = (fq_zech_equal(b, a, ctx)); if (!result) { flint_printf("FAIL:\n\n"); fq_zech_ctx_print(ctx); flint_printf("\n"); flint_printf("aa = "); fq_nmod_print_pretty(aa, ctx->fq_nmod_ctx); flint_printf("\n"); flint_printf("a = "); fq_zech_print_pretty(a, ctx); flint_printf("\n"); flint_printf("b = "); fq_zech_print_pretty(b, ctx); flint_printf("\n"); abort(); } fmpz_clear(x); fq_nmod_clear(bb, ctx->fq_nmod_ctx); fq_nmod_clear(aa, ctx->fq_nmod_ctx); } fq_zech_ctx_clear(ctx); } FLINT_TEST_CLEANUP(state); flint_printf("PASS\n"); return EXIT_SUCCESS; }