Algorithms 52-53. Determination of an interpolating quintic spline function with equally spaced and double knots
Inversion of band matrices is the simplest, direct way of computing interpolating splines. On the other hand, when studying the convergence problems, the bounds of related inverse matrices are useful. This paper contains algorithms of inversion, based on LR decomposition, of tri- and five-diagonal matrices appearing in spline fitting problems. (In this case LR decomposition is possible and unique if one of the diagonals is fixed.) For example, inversion of a tri-diagonal matrix of dimension n requires...
Page 1