Quercus Robur's questions

Tenho dois arrays numpy ae boutros que a.shape[:-1]e b.shapesão transmissíveis. Com apenas essa restrição, quero calcular um array cde acordo com o seguinte:

c = numpy.empty(numpy.broadcast_shapes(a.shape[:-1],b.shape),a.dtype)
for i in range(a.shape[-1]):
    c[...,i] = a[...,i] * b

O código acima certamente funciona, mas gostaria de saber se existe uma maneira mais elegante (e idiomática) de fazê-lo.

Suponha que eu receba alguns arrays numpy, digamos a, be c, que são assumidos como transmissíveis. Existe um padrão ou, de outra forma, uma maneira elegante de encontrar o formato do array após a transmissão? Claro, algo como (a+b+c).shapefuncionaria, mas é muito ineficiente se eu estiver interessado apenas no formato do resultado.