import numpy as np litres = np.array([12., 18., 24.]) hours = np.array([2., 0., 4.]) rates = np.full(3, np.nan) np.divide(litres, hours, out=rates, where=hours > 0) print(rates[[0, 2]].tolist()) print(np.isnan(rates[1]))