Código de computación para la resolución directa de sistemas de ecuaciones algebraicas lineales con matriz esparcida, simetrica y definida positiva: codigo SPFACT

Abstract

Una variedad de métodos numéricos utilizados en la práctica conduce a la generación de grandes sistemas de ecuaciones algebraicas lineales. En muchos casos la matriz del sistema resulta esparcida o sea que la mayor parte de sus elementos son iguales a cero. En tales casos se deben emplear procedimientos de cómputos especiales para almacenar sólo los no-ceros en la memoria de la computadora, para realizar operaciones sólo con los no-ceros, y para reducir la generaicón de nuevos no-ceros durante los cálculos tanto como sea posible. Así se asegura un empleo eficiente del tiempo y memoria de computadora disponibles. Durante los últimos años los algoritmos empleados para cálculo con matrices esparcidas se han vuelto muy sofisticados, y los cómputos se pueden realizar ahora de una manera muy densa y eficiente. Aquí presentamos un paquete de subrutina en Fortran que incorpora los últimos adelantos en lo que ha empezado a llmarse la Tecnología de Matrices Esparcidas, para resolver grandes sistemas de ecuaciones lineales con matriz simétrica y definida positiva, del tipo que aparece en el método de Elementos Finitos.