static __inline__ void
zeta_coeff_k(zeta_bsplit_t S, slong k, slong n, slong s)
{
arb_set_si(S->D, 2 * (n + k));
arb_mul_si(S->D, S->D, n - k, ARF_PREC_EXACT);
arb_set_si(S->Q1, k + 1);
arb_mul_si(S->Q1, S->Q1, 2*k + 1, ARF_PREC_EXACT);
if (k == 0)
{
arb_zero(S->A);
arb_one(S->Q2);
}
else
{
arb_set_si(S->A, k % 2 ? 1 : -1);
arb_mul(S->A, S->A, S->Q1, ARF_PREC_EXACT);
arb_ui_pow_ui(S->Q2, k, s, ARF_PREC_EXACT);
}
arb_mul(S->Q3, S->Q1, S->Q2, ARF_PREC_EXACT);
arb_zero(S->B);
arb_set(S->C, S->Q1);
}
void Lib_Arb_Mul_Si(ArbPtr f, ArbPtr g, int32_t x, int32_t prec)
{
arb_mul_si( (arb_ptr) f, (arb_ptr) g, x, prec);
}
