The scipy.special.zeta and scipy.special.zetac don't take complex inputs.

A previously closed issue on this argued that the complex use case is not supported.

However, there isn't really much of a use case for zeta() without complex variables.

This limitation is like having a euclidian distance function that only supported one-dimensional inputs.