SciPy

numpy.ndarray.__ipow__

ndarray.__ipow__

x.__ipow__(y) <==> x**y