Nonumpy as dimensões da matriz resultante variam em tempo de execução. Muitas vezes há confusão entre uma matriz 1d e uma matriz 2d com 1 coluna. Em um caso, eu posso iterar sobre as colunas; no outro, não posso.

Como você resolve elegantemente esse problema? Para evitar desarrumar meu código comif declarações checando a dimensionalidade, eu uso esta função:

def reshape_to_vect(ar):
    if len(ar.shape) == 1:
      return ar.reshape(ar.shape[0],1)
    return ar

No entanto, isso parece deselegante e caro. Existe uma solução melhor?