No complex arguments in scipy.special.zeta or scipy.special.zetac

This is a feature request. Complex numbers are not supported by the scipy.special.zeta and scipy.special.zetac functions.

Reproducing code example:

zeta(2+1j)

Error message:

---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-29-07480cb9365e> in <module>()
----> 1 zeta(2+1j)

/usr/local/lib/python3.6/dist-packages/scipy/special/basic.py in zeta(x, q, out)
   2461     if q is None:
   2462         q = 1
-> 2463     return _zeta(x, q, out)
   2464 

TypeError: ufunc '_zeta' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''

Scipy/Numpy/Python version information:

0.19.1 1.14.5 sys.version_info(major=3, minor=6, micro=3, releaselevel='final', serial=0)

[person142]

@sohambhattacharyya an important question here is what application you have in mind when requesting complex arguments. The “standard” application for complex arguments of the Riemann zeta function is of course investigating the zeros along the critical strip, and this is not at all a procedure for which SciPy is suitable. Our special functions typically only accept double precision arguments, and properly investigating zeros requires using arbitrary precision. Arbitrary precision arguments are, at least for now, out of scope for SciPy.

[rgommers]

Okay looks like that’s a no to the feature request for the time being. So closing.