zeta

→Tensor

zeta(x: Tensor, q: Tensor)

source edit

Hurwitz zeta function $\zeta(x, q)$ .

Computes the Hurwitz zeta function, a two-argument generalisation of the Riemann zeta that arises in moment-generating functions of discrete distributions, in expansions of polygamma functions, and as a regulariser in analytic number theory.

Parameters

xTensor

Exponent; must satisfy

\Re(x) > 1

for the series to converge. Any floating-point dtype.

qTensor

Shift parameter; broadcast-compatible with x. The non-positive integers are poles.

Returns

Tensor

$\zeta(x, q)$ element-wise, broadcast to the common shape of x and q.

Notes

Series definition:

\zeta(x, q) = \sum_{k=0}^\infty \frac{1}{(k + q)^x}.

The implementation accumulates an explicit prefix of 12 terms and closes the tail with the Euler–Maclaurin correction

\sum_{k=N}^\infty (k + q)^{-x} \approx \frac{a^{1-x}}{x - 1} + \frac{a^{-x}}{2} + \sum_j \frac{B_{2j}}{(2j)!} \prod_{i=0}^{2j-2}(x + i)\, a^{-x-2j+1},

where $a = q + N$ and the $B_{2j}$ are Bernoulli numbers; that is accurate to float64 rounding for $x > 1$ and $q > 0$ . A negative non-integer q is real only for an integer x and is reduced to positive arguments; otherwise it is NaN. The non-positive integers q are poles (+inf), x = 1 is +inf and x < 1 is NaN. The Riemann zeta is the special case $\zeta(x) = \zeta(x, 1)$ .

Examples

>>> import lucid
>>> from lucid.special import zeta
>>> x = lucid.tensor([2.0, 3.0])
>>> q = lucid.tensor([1.0, 1.0])
>>> zeta(x, q)
tensor([1.645, 1.202])

Used by 1

lucid.special