Skip to main content
cuerpo editado
Origen Enlace
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68
good_rowsgood_cols = numpy.any(matriz_v.T != 0,  axis = 1)
matriz_v = matriz_v[:, good_rows][good_rowsgood_cols][good_cols, :]
matriz = [[0.003,     0,     0,     -1,     0,    0,    0], 
          [    0, 0.004,     0,     -1,     0,    0,    0],
          [    0,     0,     0.1,   -1,     0,    0,    0],
          [    1,     1,     1,      0,     0,    0,    0], 
          [    0,     0,     0,      0,     0,    0,    0],
          [    0,     0,     0,      0,     0,    0,    0],
          [    0,     0,     0,      1,     0,    0,    0]]
          
        

matriz_v = numpy.array(matriz)
good_rowsgood_cols = numpy.any(matriz_v.T != 0,  axis = 1)
matriz_v = matriz_v[:, good_rows][good_rowsgood_cols][good_cols, :]
good_rows = numpy.any(matriz_v.T != 0,  axis = 1)
matriz_v = matriz_v[:, good_rows][good_rows, :]
matriz = [[0.003,     0,     0,     -1,     0,    0,    0], 
          [    0, 0.004,     0,     -1,     0,    0,    0],
          [    0,     0,     0.1,   -1,     0,    0,    0],
          [    1,     1,     1,      0,     0,    0,    0], 
          [    0,     0,     0,      0,     0,    0,    0],
          [    0,     0,     0,      0,     0,    0,    0],
          [    0,     0,     0,      1,     0,    0,    0]]
          
        

matriz_v = numpy.array(matriz)
good_rows = numpy.any(matriz_v.T != 0,  axis = 1)
matriz_v = matriz_v[:, good_rows][good_rows, :]
good_cols = numpy.any(matriz_v.T != 0,  axis = 1)
matriz_v = matriz_v[:, good_cols][good_cols, :]
matriz = [[0.003,     0,     0,     -1,     0,    0,    0], 
          [    0, 0.004,     0,     -1,     0,    0,    0],
          [    0,     0,     0.1,   -1,     0,    0,    0],
          [    1,     1,     1,      0,     0,    0,    0], 
          [    0,     0,     0,      0,     0,    0,    0],
          [    0,     0,     0,      0,     0,    0,    0],
          [    0,     0,     0,      1,     0,    0,    0]]
          
        

matriz_v = numpy.array(matriz)
good_cols = numpy.any(matriz_v.T != 0,  axis = 1)
matriz_v = matriz_v[:, good_cols][good_cols, :]
Origen Enlace
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68

El problema es que matriz_v_c y matriz_v_c_f se reddefinen por cada ciclo del primer for, esto causa que al final solo se aplique el resultado de la última iteración.

Se puede hacer de esta forma pero es complicarse la vida sin necesidad, recurda que NumPy y Python hacen 'magia':

  • Primero buscamos una manera de crear un array de boleanos en el que cada elemento será True si esa columna es buena (tiene algún elemento que no sea 0) y False en caso contrario. En tu ejemplo el array sería:

      [ True  True  True  True False False False]
    
  • Despues aplicamos slicing sobre la matriz original usando ese array de voleanos. Si, resulta que Numpy permite hacer cosas como:

      >>> import numpy as np
      >>> m = np.array([[1,2,3],
                       [4,5,6]])
      >>> filtro = np.array([True, False, True])
      >>> m[:, filtro]
      array([[1, 3],
             [4, 6]])
    

En tu caso el codigo sería simplemente:

good_rows = numpy.any(matriz_v.T != 0,  axis = 1)
matriz_v = matriz_v[:, good_rows][good_rows, :]

Se usa numpy.any que retorna True en cuanto algún elemento de la columna no es 0 (Short-circuit evaluation), lo que aumenta la eficiencia sobre sum o all.

Por si se quiere reproducir el código con tu ejemplo:

matriz = [[0.003,     0,     0,     -1,     0,    0,    0], 
          [    0, 0.004,     0,     -1,     0,    0,    0],
          [    0,     0,     0.1,   -1,     0,    0,    0],
          [    1,     1,     1,      0,     0,    0,    0], 
          [    0,     0,     0,      0,     0,    0,    0],
          [    0,     0,     0,      0,     0,    0,    0],
          [    0,     0,     0,      1,     0,    0,    0]]
          
        

matriz_v = numpy.array(matriz)
good_rows = numpy.any(matriz_v.T != 0,  axis = 1)
matriz_v = matriz_v[:, good_rows][good_rows, :]

Salida:

[[ 0.003  0.     0.    -1.   ]
 [ 0.     0.004  0.    -1.   ]
 [ 0.     0.     0.1   -1.   ]
 [ 1.     1.     1.     0.   ]]