Skip to main content
2 de 6
se eliminaron 4 caracteres en el cuerpo
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68

El problema se encuentra en:

for i in range(5): 
    linha_aux = linha_aux + 1
    if pivos_indices[i] > pivos_indices[linha_aux] and linha_aux < 5 and i < 5:

Tienes un array de 5 elementos, el valor máximo que alcanza i en tu for es 4, pero para i = 4 linha_aux vale 5 y pivos_indices[5] no es válido porque ese índice esta fuera de rango para un array de 5 elementos (en el que el indice máximo es 4). and linha_aux < 5 no sirve de protección porque no llega a evaluarse, ya que la indexación se produce antes de que se evalúe.

Lo lógico es hacer que el for itere hasta len(array) - 1, en este caso for i in range(4).

Hay una forma más simple de ordenar el array en función del numero de ceros consecutivos iniciales. Valiendonos de tu idea de usar np.nonzero podemos usar el array resultante con el indice de los pivotes. Después lo ordenamos y obtenemos los indices de los elementos ordenados (numpy.argsort) y por último hacemos un indexado de mat con él:

import numpy as np

def buscar_pivo(L):
    p = np.nonzero(L)[1]
    if p.size:
        return p[0]
    return L.size

mat = np.matrix([[0,  5,  2,  7,  8],
                 [0,  0,  0, 14, 16],
                 [0,  0,  0,  0,  0],
                 [0,  0, 10, 16, 22],
                 [0,  0,  0,  0,  8]]).astype(float)

res = mat[np.argsort(np.apply_along_axis(buscar_pivo, axis=1, arr=mat))]
print(res)

Salida:

[[  0.   5.   2.   7.   8.]
 [  0.   0.  10.  16.  22.]
 [  0.   0.   0.  14.  16.]
 [  0.   0.   0.   0.   8.]
 [  0.   0.   0.   0.   0.]]

Para tu ejemplo original obtenemos:

>>> mat = np.matrix([[0,  5,  2,  7,  8],
                     [8, 10,  4, 14, 16],
                     [0,  0,  0,  0,  0],
                     [2,  6, 10, 16, 22],
                     [3,  5,  8,  9, 15]]).astype(float)
>>> res = mat[np.argsort(np.apply_along_axis(buscar_pivo, axis=1, arr=mat))]
>>> res  
[[  8.  10.   4.  14.  16.]
 [  2.   6.  10.  16.  22.]
 [  3.   5.   8.   9.  15.]
 [  0.   5.   2.   7.   8.]
 [  0.   0.   0.   0.   0.]]
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68