SciPy

numpy.ndarray.__and__

ndarray.__and__

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