Dissertation/Thesis Abstract

The author has requested that access to this graduate work be delayed until 2019-05-15. After this date, this graduate work will be available on an open access basis.
A New Block Krylov Subspace Framework with Applications to Functions of Matrices Acting on Multiple Vectors
by Lund, Kathryn, Ph.D., Temple University, 2018, 212; 10792061
Abstract (Summary)

We propose a new framework for understanding block Krylov subspace methods, which hinges on a matrix-valued inner product. We can recast the ``classical" block Krylov methods, such as O'Leary's block conjugate gradients, global methods, and loop-interchange methods, within this framework. Leveraging the generality of the framework, we develop an efficient restart procedure and error bounds for the shifted block full orthogonalization method (Sh-BFOM( m)). Regarding BFOM as the prototypical block Krylov subspace method, we propose another formalism, which we call modified BFOM, and show that block GMRES and the new block Radau-Lanczos method can be regarded as modified BFOM. In analogy to Sh-BFOM(m), we develop an efficient restart procedure for shifted BGMRES with restarts (Sh-BGMRES(m)), as well as error bounds.

Using this framework and shifted block Krylov methods with restarts as a foundation, we formulate block Krylov subspace methods with restarts for matrix functions acting on multiple vectors f(A) B. We obtain convergence bounds for B(FOM)2 (BFOM for Functions Of Matrices) and block harmonic methods (i.e., BGMRES-like methods) for matrix functions.

With various numerical examples, we illustrate our theoretical results on Sh-BFOM and Sh-BGMRES. We also analyze the matrix polynomials associated to the residuals of these methods. Through a variety of real-life applications, we demonstrate the robustness and versatility of B(FOM)2 and block harmonic methods for matrix functions. A particularly interesting example is the tensor t-function, our proposed definition for the function of a tensor in the tensor t-product formalism. Despite the lack of convergence theory, we also show that the block Radau-Lanczos modification can reduce the number of cycles required to converge for both linear systems and matrix functions.

Indexing (document details)
Advisor: Szyld, Daniel B., Frommer, Andreas
Commitee: Bolten, Matthias, Queisser, Gillian, Seibold, Benjamin
School: Temple University
Department: Mathematics
School Location: United States -- Pennsylvania
Source: DAI-B 79/09(E), Dissertation Abstracts International
Subjects: Applied Mathematics
Keywords: Block Krylov Subspace Methods, Functions of matrices, Functions of tensors, Matrix polynomials, Shifted linear systems, Tensor t-product
Publication Number: 10792061
ISBN: 9780355954814
Copyright © 2019 ProQuest LLC. All rights reserved. Terms and Conditions Privacy Policy Cookie Policy