The browser you are using is not supported by this website. All versions of Internet Explorer are no longer supported, either by us or Microsoft (read more here: https://www.microsoft.com/en-us/microsoft-365/windows/end-of-ie-support).

Please use a modern browser to fully experience our website, such as the newest versions of Edge, Chrome, Firefox or Safari etc.

Portrait of Krzysztof Podgórski. Photo.

Krzysztof Podgórski

Professor, Head of the Department of Statistics

Portrait of Krzysztof Podgórski. Photo.

Dyadic diagonalization of positive definite band matrices and efficient B-spline orthogonalization

Author

  • Xijia Liu
  • Hiba Nassar
  • Krzysztof Podgórski

Summary, in English

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 applicable to a Toeplitz band Gramian while the second one is more general, works with any Gramian matrix, but is more computationally intensive. In the context of the B-splines, these two cases result in new symmetric orthogonalization procedures and correspond to equally and arbitrarily spaced knots, respectively. In the algorithm, the sparsity of a band Gramian is utilized to produce a natural dyadic net of orthogonal splines, rather than a sequence of them. Such a net is thus naturally referred to as a splinet. The splinets exploit “near-orthogonalization” of the B-splines and feature locality expressed through a small size of the total support set and computational efficiency that is a result of a small number of inner product evaluations needed for their construction. These and other efficiencies are formally quantified by upper bounds and asymptotic rates with respect to the number of splines in a splinet. An additional assessment is provided through numerical experiments. They suggest that the theoretical bounds are rather conservative and the method is even more efficient than the bounds indicate. The dyadic net-like structures and the locality bear some resemblance to wavelets but in fact, the splinets are fundamentally different because they do not aim at capturing the resolution scales. The orthogonalization method together with efficient spline algebra and calculus has been implemented in R-package Splinets available on CRAN.

Department/s

  • Department of Statistics

Publishing year

2022-11

Language

English

Publication/Series

Journal of Computational and Applied Mathematics

Volume

414

Document type

Journal article

Publisher

Elsevier

Topic

  • Control Engineering

Keywords

  • B-splines
  • Band matrices
  • Dyadic structure
  • Matrix inversion
  • Orthogonalization

Status

Published

ISBN/ISSN/Other

  • ISSN: 0377-0427