def divide(x, y): with np.errstate(divide='ignore', invalid='ignore'): z = np.true_divide(x, y) z[~ np.isfinite(z)] = 0 return z