y = numpy.arange(3 * 4).reshape(3, 4)
z = numpy.arange(4 * 4).reshape(4, 4)
(x.reshape(2, 1, 1, 4) * y.reshape(1, 3, 1, 4) * z.reshape(1, 1, 4, 4)).sum(
    axis=3)
Как я могу векторизовать этот тройной цикл над двумерными массивами в numpy?