int
f_elliptic_p_laurent_n(acb_ptr res, const acb_t z, void * param, slong order, slong prec)
{
slong n;
acb_t tau;
if (order > 1)
flint_abort(); /* Would be needed for Taylor method. */
n = ((slong *)(param))[0];
acb_init(tau);
acb_onei(tau);
acb_modular_elliptic_p(res, z, tau, prec);
acb_pow_si(tau, z, -n - 1, prec);
acb_mul(res, res, tau, prec);
acb_clear(tau);
return 0;
}
int main()
{
slong iter;
flint_rand_t state;
flint_printf("elliptic_p....");
fflush(stdout);
flint_randinit(state);
/* check test values */
for (iter = 0; iter < 100; iter++)
{
slong i;
acb_t z, tau, p1, p2;
acb_init(z);
acb_init(tau);
acb_init(p1);
acb_init(p2);
for (i = 0; i < NUM_TESTS; i++)
{
acb_set_dddd(z, testdata[i][0], 0.0, testdata[i][1], 0.0);
acb_set_dddd(tau, testdata[i][2], 0.0, testdata[i][3], 0.0);
acb_set_dddd(p2, testdata[i][4], EPS, testdata[i][5], EPS);
acb_modular_elliptic_p(p1, z, tau, 2 + n_randint(state, 1000));
if (!acb_overlaps(p1, p2))
{
flint_printf("FAIL (test value)\n");
flint_printf("tau = "); acb_printd(tau, 15); flint_printf("\n\n");
flint_printf("z = "); acb_printd(z, 15); flint_printf("\n\n");
flint_printf("p1 = "); acb_printd(p1, 15); flint_printf("\n\n");
flint_printf("p2 = "); acb_printd(p2, 15); flint_printf("\n\n");
abort();
}
}
acb_clear(z);
acb_clear(tau);
acb_clear(p1);
acb_clear(p2);
}
/* Test periods */
for (iter = 0; iter < 2000; iter++)
{
acb_t tau, z1, z2, p1, p2;
slong m, n, e0, prec0, prec1, prec2;
acb_init(tau);
acb_init(z1);
acb_init(z2);
acb_init(p1);
acb_init(p2);
e0 = 1 + n_randint(state, 10);
prec0 = 2 + n_randint(state, 1000);
prec1 = 2 + n_randint(state, 1000);
prec2 = 2 + n_randint(state, 1000);
acb_randtest(tau, state, prec0, e0);
if (arf_sgn(arb_midref(acb_imagref(tau))) < 0)
acb_neg(tau, tau);
acb_randtest(z1, state, prec0, e0);
acb_randtest(p1, state, prec0, e0);
acb_randtest(p2, state, prec0, e0);
/* z2 = z1 + m + n*tau */
m = n_randint(state, 10);
n = n_randint(state, 10);
acb_add_ui(z2, z1, m, prec0);
acb_addmul_ui(z2, tau, n, prec0);
acb_modular_elliptic_p(p1, z1, tau, prec1);
acb_modular_elliptic_p(p2, z2, tau, prec2);
if (!acb_overlaps(p1, p2))
{
flint_printf("FAIL (overlap)\n");
flint_printf("tau = "); acb_printd(tau, 15); flint_printf("\n\n");
flint_printf("z1 = "); acb_printd(z1, 15); flint_printf("\n\n");
flint_printf("z2 = "); acb_printd(z2, 15); flint_printf("\n\n");
flint_printf("p1 = "); acb_printd(p1, 15); flint_printf("\n\n");
flint_printf("p2 = "); acb_printd(p2, 15); flint_printf("\n\n");
abort();
}
acb_modular_elliptic_p(z1, z1, tau, prec1);
if (!acb_overlaps(z1, p1))
{
flint_printf("FAIL (aliasing)\n");
flint_printf("tau = "); acb_printd(tau, 15); flint_printf("\n\n");
flint_printf("z1 = "); acb_printd(z1, 15); flint_printf("\n\n");
flint_printf("p1 = "); acb_printd(p1, 15); flint_printf("\n\n");
abort();
}
acb_clear(tau);
acb_clear(z1);
acb_clear(z2);
acb_clear(p1);
acb_clear(p2);
}
flint_randclear(state);
flint_cleanup();
flint_printf("PASS\n");
return EXIT_SUCCESS;
}
