Cálculo de matrices:

Tengo que resolver el cálculo de grandes matrices y desearía saber si
alguien puede orientarme en la búsqueda de bibliografía y alguna subrutina
en "C" que me permita abordar la solución.

La matriz a calcular tiene las siguientes características:

* Matriz cuadrada.
* El orden "n" puede ser de hasta 10.000.
* La diagonal principal tiene coeficientes distintos de cero y su
valor absoluto mayor a 1.
* Cada una de las filas tiene un número muy reducido de coeficientes
diferente de cero (entre 3 a 6).
* Los coeficientes diferentes de cero de las filas no se agrupan
alrededor de la diagonal, ni tampoco en sus extremos, por lo general se
encuentran distribuidos en forma muy separada dentro de la misma fila.
*
* Dado que el cálculo del problema es recurrente, la rutina de cálculo
para resolver la matriz debería permitir:
*
* Obtener los coeficientes independientes en un número reducido de
iteraciones.
* Dentro de lo posible solo debería considerar los coeficientes
diferentes de cero para ahorrar memoria.

A la fecha solo encontré como mejor método el de Gauss-Seidel, pero este
tiene la desventaja del gran número de iteraciones necesarias para obtener
los coeficientes de los términos independientes.

Mi correo es garciahugomario@...