Декартово произведение точек массива x и y на один массив 2D-точек

У меня есть два массива numpy, которые определяют оси x и y сетки. Например:

x = numpy.array([1,2,3])
y = numpy.array([4,5])

Я хотел бы сгенерировать декартово произведение этих массивов, чтобы сгенерировать:

array([[1,4],[2,4],[3,4],[1,5],[2,5],[3,5]])

Таким образом, что это не так уж неэффективно, поскольку мне нужно делать это много раз в цикле. Я предполагаю, что преобразование их в список Python и использование itertools.productи обратно в массив numpy — не самая эффективная форма.