int main() { long iter; flint_rand_t state; printf("det...."); fflush(stdout); flint_randinit(state); for (iter = 0; iter < 10000; iter++) { fmpq_mat_t Q; fmpq_t Qdet; acb_mat_t A; acb_t Adet, imagunit; long n, qbits, prec; int imaginary; n = n_randint(state, 8); qbits = 1 + n_randint(state, 100); prec = 2 + n_randint(state, 200); imaginary = n_randint(state, 2); fmpq_mat_init(Q, n, n); fmpq_init(Qdet); acb_mat_init(A, n, n); acb_init(Adet); acb_init(imagunit); fmpq_mat_randtest(Q, state, qbits); fmpq_mat_det(Qdet, Q); acb_mat_set_fmpq_mat(A, Q, prec); if (imaginary) { acb_onei(imagunit); acb_mat_scalar_mul_acb(A, A, imagunit, prec); } acb_mat_det(Adet, A, prec); if (imaginary) { acb_onei(imagunit); acb_inv(imagunit, imagunit, prec); acb_pow_ui(imagunit, imagunit, n, prec); acb_mul(Adet, Adet, imagunit, prec); } if (!acb_contains_fmpq(Adet, Qdet)) { printf("FAIL (containment, iter = %ld)\n", iter); printf("n = %ld, prec = %ld\n", n, prec); printf("\n"); printf("Q = \n"); fmpq_mat_print(Q); printf("\n\n"); printf("Qdet = \n"); fmpq_print(Qdet); printf("\n\n"); printf("A = \n"); acb_mat_printd(A, 15); printf("\n\n"); printf("Adet = \n"); acb_printd(Adet, 15); printf("\n\n"); printf("Adet = \n"); acb_print(Adet); printf("\n\n"); abort(); } fmpq_mat_clear(Q); fmpq_clear(Qdet); acb_mat_clear(A); acb_clear(Adet); acb_clear(imagunit); } flint_randclear(state); flint_cleanup(); printf("PASS\n"); return EXIT_SUCCESS; }
int main() { slong iter; flint_rand_t state; flint_printf("exp...."); fflush(stdout); flint_randinit(state); /* check exp(A)*exp(c*A) = exp((1+c)*A) */ for (iter = 0; iter < 500 * arb_test_multiplier(); iter++) { acb_mat_t A, E, F, EF, G; fmpq_mat_t Q; acb_t c, d; slong n, qbits, prec; n = n_randint(state, 5); qbits = 2 + n_randint(state, 300); prec = 2 + n_randint(state, 300); acb_init(c); acb_init(d); fmpq_mat_init(Q, n, n); acb_mat_init(A, n, n); acb_mat_init(E, n, n); acb_mat_init(F, n, n); acb_mat_init(EF, n, n); acb_mat_init(G, n, n); _fmpq_mat_randtest_for_exp(Q, state, qbits); acb_mat_set_fmpq_mat(A, Q, prec); acb_mat_exp(E, A, prec); acb_randtest(c, state, prec, 10); acb_mat_scalar_mul_acb(F, A, c, prec); acb_mat_exp(F, F, prec); acb_add_ui(d, c, 1, prec); acb_mat_scalar_mul_acb(G, A, d, prec); acb_mat_exp(G, G, prec); acb_mat_mul(EF, E, F, prec); if (!acb_mat_overlaps(EF, G)) { flint_printf("FAIL\n\n"); flint_printf("n = %wd, prec = %wd\n", n, prec); flint_printf("c = \n"); acb_printd(c, 15); flint_printf("\n\n"); flint_printf("A = \n"); acb_mat_printd(A, 15); flint_printf("\n\n"); flint_printf("E = \n"); acb_mat_printd(E, 15); flint_printf("\n\n"); flint_printf("F = \n"); acb_mat_printd(F, 15); flint_printf("\n\n"); flint_printf("E*F = \n"); acb_mat_printd(EF, 15); flint_printf("\n\n"); flint_printf("G = \n"); acb_mat_printd(G, 15); flint_printf("\n\n"); flint_abort(); } acb_clear(c); acb_clear(d); fmpq_mat_clear(Q); acb_mat_clear(A); acb_mat_clear(E); acb_mat_clear(F); acb_mat_clear(EF); acb_mat_clear(G); } flint_randclear(state); flint_cleanup(); flint_printf("PASS\n"); return EXIT_SUCCESS; }