numpy.meshgrid

x, y = np.meshgrid([1,2],[3,4,5,6])

x =

array([[1, 2],
           [1, 2],
           [1, ②],
           [1, 2]])

y =

array([[3, 3],
           [4, 4],
           [5, ⑤],
           [6, 6]])

可能な格子点が得られるといった感じ？（x,yの同じ位置の組が格子点）