SciPy

numpy.ndarray.__iand__

ndarray.__iand__

x.__iand__(y) <==> x&=y