Dyadic diagonalization of positive definite band matrices and efficient B-spline orthogonalization
A dyadic algorithm for diagonalizing an arbitrary positive definite band matrix, referred to as a band Gramian, is obtained to efficiently orthogonalize the B-splines. The algorithm can be also used as a fast inversion method for a band Gramian characterized by remarkable sparsity of the diagonalizing matrix. There are two versions of the algorithm: the first one is more efficient and is applicabl
