In this paper we analyze the behavior of the LSQR algorithm for the solution of compact operator equations in Hilbert spaces. We present results concerning existence of Krylov solutions and the rate of convergence in terms of an ℓp sequence where p depends on the summability of the singular values of the operator. Under stronger regularity requirements we also consider the decay of the error. Finally we study the approximation of the dominant singular values of the operator attainable with the bidiagonal matrices generated by the Lanczos bidiagonalization and the arising low rank approximations. Some numerical experiments on classical test problems are presented.

%B Linear Algebra and its Applications %V 583 %P 146-164 %G eng %U https://www.sciencedirect.com/science/article/pii/S0024379519303714 %R https://doi.org/10.1016/j.laa.2019.08.024 %0 Journal Article %J Calcolo %D 2019 %T On Krylov solutions to infinite-dimensional inverse linear problems %A Noe Caruso %A Alessandro Michelangeli %A Paolo Novati %XWe discuss, in the context of inverse linear problems in Hilbert space, the notion of the associated infinite-dimensional Krylov subspace and we produce necessary and sufficient conditions for the Krylov-solvability of a given inverse problem, together with a series of model examples and numerical experiments.

%B Calcolo %V 56 %P 1–25 %G eng %0 Report %D 2018 %T Truncation and convergence issues for bounded linear inverse problems in Hilbert space %A Noe Caruso %A Alessandro Michelangeli %A Paolo Novati %X We present a general discussion of the main features and issues that (bounded) inverse linear problems in Hilbert space exhibit when the dimension of the space is infinite. This includes the set-up of a consistent notation for inverse problems that are genuinely infinite-dimensional, the analysis of the finite-dimensional truncations, a discussion of the mechanisms why the error or the residual generically fail to vanish in norm, and the identification of practically plausible sufficient conditions for such indicators to be small in some weaker sense. The presentation is based on theoretical results together with a series of model examples and numerical tests. %I SISSA %G en %U http://preprints.sissa.it/handle/1963/35326 %1 35637 %2 Mathematics %4 1 %$ Submitted by Maria Pia Calandra (calapia@sissa.it) on 2018-11-20T13:02:23Z No. of bitstreams: 1 infinite_dim_truncation_sissapreprint.pdf: 1326956 bytes, checksum: 75a6af69b0bca0c5b9b5283e640d89be (MD5)