SOLVEBLOK: A Package for Solving Almost Block Diagonal Linear Systems
Carl de Boor, Richard G. Weiss · ACM Transactions on Mathematical Software · 1980
In spline approximation and m the numerical solution of (systems of) orchnary differential equations, one IS eventually faced with linear systems of a rather characteristic, almost block diagonal, form which reflects the underlying geometry of the problems.Whde these systems could be treated as banded systems, it is more efficient, both m the generation of the systems and m their numerical solution, to work explicitly with their block structure.The paper describes a package of Fortran routines for the solutmn of an almost block diagonal system by Gaussian eliminatmn with scaled row pivoting In addition, alternative approaches to the efficient solution of such systems are discussed Key Words and Phrases" almost block diagonal systems, Gaussian eliminatmn, splme approximation, ordinary differential equations CR Categories: 5.13, 5.14, 5.17 The Algorithm.SOLVEBLOK.