Python Numpy中的组合学

Question

我想编写一个单个功能，以输出2个矩阵的所有可能组合：

def combine(*args):
    return np.array(np.meshgrid(args)).T.reshape(-1, len(args)+1)

但是当传递时：

print(combine(np.array([1,2,3]), np.array([4,5,6])))

输出：

[[1 2 3]
 [4 5 6]]

如何使其工作？我想保持它自动，而不是简单地通过（args [0]，args [1]）

Answer 1

Python Itertools的直接使用：

In [134]: import itertools
In [135]: a,b = [1,2,3], [4,5,6]
In [137]: list(itertools.product(a,b))
Out[137]: [(1, 4), (1, 5), (1, 6), (2, 4), (2, 5), (2, 6), (3, 4), (3, 5), (3, 6)]

或作为数组：

In [145]: np.array(list(itertools.product(a,b)))
Out[145]: 
array([[1, 4],
       [1, 5],
       [1, 6],
       [2, 4],
       [2, 5],
       [2, 6],
       [3, 4],
       [3, 5],
       [3, 6]])