Estoy usando uBLAS de Boost en un código numérico y tengo un solucionador 'pesado' en su lugar:

http://www.crystalclearsoftware.com/cgi-bin/boost_wiki/wiki.pl?LU_Matrix_Inversion

El código funciona excelentemente, sin embargo, es dolorosamente lento. Después de un poco de investigación, encontréUMFPACK, que es un solucionador de matriz disperso (entre otras cosas). Mi código genera grandes matrices dispersas que necesito invertir con mucha frecuencia (resuelva más correctamente, el valor de la matriz inversa es irrelevante), por lo que UMFPACk y la clase Sparse_Matrix de BOOST parece ser un matrimonio feliz.

UMFPACK solicita la matriz dispersa especificada por tres vectores: un recuento de entradas, índices de filas y las entradas. (Ver ejemplo)

Mi pregunta se reduce a: ¿puedo obtener estos tres vectores de manera eficiente de la clase de matriz dispersa de BOOST?