dist/src/eta.c

0 /* eta -- Functions for computing the Dedekind eta function
26 eta_series (mpcb_ptr eta, mpcb_srcptr q, mpfr_exp_t expq, int N)
27    /* Evaluate 2N+1 terms of the Dedekind eta function without the q^(1/24)
47    mpcb_set_ui_ui (eta, 1, 0, p);
54    mpcb_add (eta, eta, q3nm1);
55    mpcb_add (eta, eta, q3np1);
63       mpcb_add (eta, eta, q3nm1);
64       mpcb_add (eta, eta, q3np1);
74    mpcr_mul (r2, r, eta->r);
75    mpcr_add (eta->r, eta->r, r);
76    mpcr_add (eta->r, eta->r, r2);
87 mpcb_eta_q24 (mpcb_ptr eta, mpcb_srcptr q24)
90       the function computes eta (z).
111       mpcb_set_inf (eta);
116          mpcb_set_inf (eta);
121          eta_series (eta, q, expq, N);
122          mpcb_mul (eta, eta, q24);
226 mpcb_eta_err (mpcb_ptr eta, mpc_srcptr z, unsigned long int err_re,
228    /* Given z=x+i*y in the fundamental domain, compute eta (z).
243    mpcb_eta_q24 (eta, q24);
252       |Re z| <= 1/2 and |z| >= 1, compute Dedekind eta (z).
258    mpcb_t eta;
262    mpcb_init (eta);
269       mpcb_eta_err (eta, zl, 0, 0);
272          ok = mpcb_can_round (eta, MPC_PREC_RE (rop), MPC_PREC_IM (rop),
289          ok = mpfr_zero_p (mpc_imagref (eta->c));
290          mpcb_add (eta, eta, fuzzb);
291          ok &= mpcb_can_round (eta, MPC_PREC_RE (rop), 2, rnd);
300       inex = mpcb_round (rop, eta, rnd);
302       inex = MPC_INEX (mpfr_set (mpc_realref (rop), mpc_realref (eta->c),
307    mpcb_clear (eta);