|
|
Reorthogonalization strategies in Golub-Kahan iterative bidiagonalization
Šmelík, Martin ; Hnětynková, Iveta (advisor) ; Kučera, Václav (referee)
The main goal of this thesis is to describe Golub-Kahan iterative bidiagonalization and its connection with Lanczos tridiagonalization and Krylov space theory. The Golub-Kahan iterative bidiagonalization is based on short recurrencies and when computing in finite precision arithmetics, the loss of orthogonality often occurs. Consequently, with the aim to reduce the loss of orthogonality, we focus on various reorthogonalization strategies. We compare them in numerical experiments on testing matrices available in the MATLAB environment. We study the dependency of the loss of orthogonalization and computational time on the choice of the method or the attributes of the matrix.
|
|
|
Regularization properties of Krylov subspace methods
Kučerová, Andrea ; Hnětynková, Iveta (advisor) ; Kučera, Václav (referee)
The aim of this thesis is to study and describe regularizing properties of iterative Krylov subspace methods for finding a solution of linear algebraic ill- posed problems contaminated by white noise. First we explain properties of this kind of problems, especially their sensitivity to small perturbations in data. It is shown that classical methods for solving approximation problems (such as the least squares method) fail here. Thus we turn to explanation of regularizing pro- perties of projections onto Krylov subspaces. Basic Krylov regularizing methods are considered, namely RRGMRES, CGLS, and LSQR. The results are illustrated on model problems from Regularization toolbox in MATLAB. 1
|
|
|
Reorthogonalization strategies in Golub-Kahan iterative bidiagonalization
Šmelík, Martin ; Hnětynková, Iveta (advisor) ; Kučera, Václav (referee)
The main goal of this thesis is to describe Golub-Kahan iterative bidiagonalization and its connection with Lanczos tridiagonalization and Krylov space theory. The Golub-Kahan iterative bidiagonalization is based on short recurrencies and when computing in finite precision arithmetics, the loss of orthogonality often occurs. Consequently, with the aim to reduce the loss of orthogonality, we focus on various reorthogonalization strategies. We compare them in numerical experiments on testing matrices available in the MATLAB environment. We study the dependency of the loss of orthogonalization and computational time on the choice of the method or the attributes of the matrix.
|
guest ::
