Solution of linear equations with skyline-stored symmetric matrix

Abstract

Fortran IV subroutines for the in-core solution of linear algebraic systems with a sparse, symmetric, skyline-stored coefficient matrix are presented. Such systems arise in a variety of applications, notably the numerical discretization of conservative physical systems by finite differences or finite element techniques. The routines can be used for processing constrained systems without need for prearranging equations. The application to ‘superelement’ condensation of large-scale systems is discussed.