Computing with NumPy arrays