It would be nice to have a pure imaginary exp implemented like npy_cexp, but only for pure imaginary arguments (maybe npy_iexp or npy_ciexp). It would take a real value $\theta$ and evaluate $e^{i\cdot\theta}$. The hope is this would be a little bit faster in cases where it is know there is no real part to the complex number.