[pyar] Matriz linealmente independiente

[Mauricio Cleveland]

Hola, escribo para consultar si alguien ha trabajado con numpy, necesito
encontrar una sub-matriz linealmente independiente.

Actualmente uso esto:

import itertools
import numpy.linalg

r = numpy.array([[1,2,2,6],[2,5,4,17],[3,9,6,1.15]])

def busca_li(A):
    N,M = A.shape
    for colset in itertools.combinations(xrange(M), N):
        B = A[:,colset]
        if numpy.linalg.det(m) != 0:
            return B, colset



print busca_li(r)

Pero la matriz que ingreso es de 250 columnas y 400 filas, se demora una
eternidad.

Conocen algún metodo mas eficiente?

Gracias de antemano!

-- 
Mauricio Cleveland
Web: http://www.universodigital.cl
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <http://listas.python.org.ar/pipermail/pyar/attachments/20101029/2a586a26/attachment.html>