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. ]]