[SOLUCIONADO] Invertir una matriz usando CVXOPT

[ser_saso]

Que tal, pues el tema en cuestión es simple, como invertir una matriz usando específicamente ese paquete de python.
Puedo invertir una matriz usando el paquete numpy, sin embargo como uso matrices dispersas en mis códigos, mi asesor me dió la tarea de hacerlo usando CVXOPT.

El comando que se usa es:

cvxopt.lapack.getri(A, ipiv)

Sin embargo se me pide como argumento la matriz a ser invertida (A) y un argumento llamado ipiv el cual es descrito como:

The argument ipiv is an integer matrix of length at least min{m, n}. ///(El argumento ipiv es una matrix entera de longitud al menos min{m, n}.)

y como verán no tengo idea a que se refiere este último argumento.
Si alguien ha usado este paquete o conoce otra manera de invertir una matriz le agradecería me comente como hacerlo.

Saludos y gracias
Sergio Salazar

[razpeitia]

Lee la documentación.

cvxopt.lapack.getri(A, ipiv)
Computes the inverse of a matrix.

On entry, A and ipiv must contain the factorization as computed by gesv or getrf. On exit, A contains the matrix inverse.

Ejemplo con una matriz simple de 2x2

Código Python:

from cvxopt import matrix, normal
from cvxopt.lapack import getrf, getri
 
n = 2
#A = normal(n, n)
A = matrix([4., 3., 3., 2.], (n,n))
old = matrix(A)
ipiv = matrix(0, (n,1))
getrf(A, ipiv)
getri(A, ipiv)
print old
print A
print old * A
 

Lo que hace primero es factorizar, y eso lo pone en ipiv, después con A e ipiv sacas la inversa.

[ser_saso]

mmm ok ok, pero aún tengo la duda si ipiv significa algo en específico... es decir porque la necesidad imperiosa de pedirnos ese argumento...
De cualquier manera gracias...

[razpeitia]

Si, ipiv son las matrices U y L (combinadas) resultado de la factorización LU.

Lo hacen porque necesitan esa matriz para obtener la matriz inversa de manera eficiente. De lo contrario seria muy lenta.

[ser_saso]

ahora todo tiene sentido...
ipiv la incializamos en cero para después de la factorización usar la matriz de permutación guardada en ipiv en las operaciones para el calculo de la inversa.

Muchas gracias!!