Abstract
In this paper we use the Lanczos process for preconditioning discrete ill-posed problems. We show that by few steps of this process one can obtain a well qualified and efficient preconditioner. This is a general method in the sense that it is not limited only to special structured matrices and the matrix–vector multiplications can be carried out in O(n) operations. Also even in problems with structured matrices this preconditioner performs more efficiently than the circulant and Kronecker product approximate preconditioners.
Similar content being viewed by others
References
Bjorck A (1996) Numerical methods for least squares problems. SIAM, Philadelphia
Calvetti D, Golub GH, Reichel L (1999) Estimation of the L-curve via Lanczos bidiagonalization. BIT 39: 603–619
Calvetti D, Reichel L, Shuibi A (2004) L-curve and curvature bounds for Tikhonov regularization. Numer Algorithms 35(2–4): 301–314
Engl HW, Hanke M, Neubauer A (1996) Regularization of inverse problems. Kluwer, Dordrecht
Golub GH, Luk FT, Overton ML (1981) A block Lanczos method for computing the singular values and corresponding singular vectors of a matrix. AMC Trans Math Softw 7(2): 149–169
Golub GH, Reinsch C (1970) Singular value decomposition and least squares solutions. Numer Math 14(5): 403–420
Golub GH, Van Loan CF (1996) Matrix computations, 3rd edn. The Johns Hopkins University Press, Baltimore
Hanke M (1995) Conjugate gradient type methods for ill-posed problems. In: Pitman research notes in mathematics, no 327. Longman, Harlow
Hanke M, Hansen PC (1993) Regularization methods for large-scale problems. Surv Math Ind 3(4): 253–315
Hanke M, Nagy J (1998) Inverse Toeplitz preconditioners for ill-posed problems. Linear Algebra Appl 284: 137–156
Hanke M, Nagy J, Plemmons R (1993) Preconditioned iterative regularization methods for ill-posed problems. In: Reichel L, Ruttan A, Varga RS (eds) Numerical linear algebra. de Gruyter, Berlin, pp 141–163
Hansen PC (1994) Regularization tools: a MATLAB package for analysis and solution of discrete ill-posed problems. Numer Algorithms 6:1–35. Software is available in Netlib at http://www.netlib.org
Hansen PC (1998) Rank-deficient and discrete ill-posed problems. SIAM, Philadelphia
Hansen PC, Jensen TK, Rodriguez G (2007) An adaptive pruning algorithm for the discrete L-curve criterion. J Comput Appl Math 198: 483–492
Kamm J, Nagy JG (1998) Kronecker product and SVD approximations in image restoration. Linear Algebra Appl 284: 177–192
Kamm J, Nagy JG (2000) Optimal Kronecker product approximation of block Toeplitz matrices. SIAM J Matrix Anal Appl 22: 155–172
Kilmer ME, Oleary DP (2001) Choosing regularization parameter in iterative methods for ill-posed problems. SIAM J Matrix Anal Appl 22: 1204–1221
Nagy JG, Kilmer ME (2006) Kronecker product approximation for preconditioning in three-dimensional imaging applications. IEEE Trans Image Process 15: 604–614
Paige CC (1974) Bidiagonalization of matrices and solution of linear equations. SIAM J Numer Anal 11: 197–209
RestoreTools (2002) An object oriented Matlab package for image restoration. http://www.mathcs.emory.edu/nagy/RestoreTools
Saad Y (2004) Iterative methods for sparse linear systems, 2nd edn. SIAM, Philadelphia
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by X. Chen.
Rights and permissions
About this article
Cite this article
Rezghi, M., Hosseini, S.M. Lanczos based preconditioner for discrete ill-posed problems. Computing 88, 79–96 (2010). https://doi.org/10.1007/s00607-010-0090-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00607-010-0090-3