pub fn zeta_core<P, E, V, const ZETAC: bool, const DERIV: bool>(s: V) -> (V, V)where
E: FloatElement + ZetaConsts + BernoulliNumbers,
V: FloatVector<Element = E> + SpecializedTranscendentalMath<E> + SpecializedSpecialMath<E>,
P: Policy,Expand description
zeta(s) (ZETAC = false) or zeta(s) - 1 (ZETAC = true), and its derivative when
DERIV is set. zeta and zetac differ by a constant, so the one derivative serves both.