wi, th, hi = np.broadcast_arrays(w, theta, h)