Thanks for looking into it!

`cpow` is indeed complicated.  We had some discussion in NumPy about `0**y` branch cuts (I did yet not finish it, because thinking about the branch cuts is tricky).

It is almost reassuring that the C standard also hedges out :).  Although, it means that there is no reliance on math libraries if we want to treat the special cases more carefully :(.