def get_meshgrid(xmin, xmax, nx, xdim):
    x1d = np.linspace(xmin, xmax, nx)
    vals = ([x1d] * xdim)
    xds = np.meshgrid(*vals)
    xs = np.concatenate([xd.reshape((- 1), 1) for xd in xds], axis=1)
    return xs
def get_meshgrid(xmin, xmax, nx, xdim):
    x1d = np.linspace(xmin, xmax, nx)
    vals = ([x1d] * xdim)
    xds = np.meshgrid(*vals)
    xs = np.concatenate([xd.reshape((- 1), 1) for xd in xds], axis=1)
    return xs