SciPy

numpy.ndarray.__imul__

ndarray.__imul__

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