x = (1 + np.linspace(0, 5, np.prod(xs))).astype(dtype).reshape(xs)