NumCosmoMathBinSplit

Binary splitting of rational series in arbitrary precision.

A NcmBinSplit holds the partial sum $$S(n_1, n_2) = \sum_{n=n_1}^{n_2-1} \frac{a(n)}{b(n)} \prod_{k=n_1}^{n} \frac{p(k)}{q(k)}$$ of a series with integer $p$, $q$, $a$ and $b$ as the GMP integers $T$, $B$, $Q$ and $P$, with $P = \prod p(k)$, $Q = \prod q(k)$, $B = \prod b(k)$ and $S = T / (BQ)$. The evaluators of a series, NcmBinSplitEval, are generated by including ncm_binsplit_eval.c, see there. ncm_binsplit_join() combines consecutive ranges, and ncm_binsplit_eval_prec() extends a sum until the next term is below a relative precision.