A Levinson-like algorithm for symmetric strongly nonsingular higher order semiseparable plus band matrices (Articolo in rivista)

Type
Label
  • A Levinson-like algorithm for symmetric strongly nonsingular higher order semiseparable plus band matrices (Articolo in rivista) (literal)
Anno
  • 2007-01-01T00:00:00+01:00 (literal)
Alternative label
  • Vandebril R., Mastronardi N., Van Barel M. (2007)
    A Levinson-like algorithm for symmetric strongly nonsingular higher order semiseparable plus band matrices
    in Journal of computational and applied mathematics
    (literal)
Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#autori
  • Vandebril R., Mastronardi N., Van Barel M. (literal)
Pagina inizio
  • 75 (literal)
Pagina fine
  • 97 (literal)
Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#numeroVolume
  • 198 (literal)
Rivista
Note
  • ISI Web of Science (WOS) (literal)
Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#affiliazioni
  • K.U.Leuven, Dept. Computerwetenschappen, Celestijnenlaan 200A, 3000 Leuven (Heverlee), Belgium. Istituto per le Applicazioni del Calcolo ”M.Picone” sez. Bari, CNR, via G. Amendola 122/D, I-70126 Bari, Italy (literal)
Titolo
  • A Levinson-like algorithm for symmetric strongly nonsingular higher order semiseparable plus band matrices (literal)
Abstract
  • In this paper, we will derive a solver for a symmetric strongly nonsingular higher order generator representable semiseparable plus band matrix. The solver we will derive is based on the Levinson algorithm, which is used for solving strongly nonsingular Toeplitz systems. In the first part an $O(p^2n)$ solver for a semiseparable matrix of semiseparability rank $ p$ is derived, and in a second part we derive an $O(l^2n) $solver for a band matrix with bandwidth $ 2l + 1.$ Both solvers are constructed in a similar way: firstly a Yule–Walker-like equation needs to be solved, and secondly this solution is used for solving a linear equation with an arbitrary right-hand side. Finally, a combination of the above methods is presented to solve linear systems with semiseparable plus band coefficient matrices. The overall complexity of this solver is (l + p)^2 n $ plus lower order terms. In the final section numerical experiments are performed. Attention is paid to the timing and the accuracy of the described methods. (literal)
Prodotto di
Autore CNR
Insieme di parole chiave

Incoming links:


Autore CNR di
Prodotto
Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#rivistaDi
Insieme di parole chiave di
data.CNR.it