是否有沿轴执行numpy.outer的矢量化方法？

Question

>>> x = np.array([['a0', 'a1'],['b0','b1']])

>>> y = np.array([['x0', 'x1'],['y0','y1']])

>>> iterable = [np.outer(x[i],y[i]) for i in xrange(x.shape[0])]

>>> elbareti = np.asarray(iterable)

>>> elbareti

array([[[ 'a0'*'x0', 'a0'*'x1' ],

        [ 'a1'*'x0', 'a1'*'x1' ]],

       [[ 'b0'*'y0', 'b0'*'y1' ],

        [ 'b1'*'y0', 'b1'*'y1' ]]])

由于我正在计划使用大型阵列，因此是否有更多类似numpy的版本？我觉得答案就在我的鼻子底下，我在想这与有关reduce，但是numpy的版本仅适用于ufuncs，而不适用于函数。甚至一个提示将不胜感激。