Large-scale linear discrete ill-posed problems are generally solved by Krylov subspace iterative methods. However, these methods can be difficult to implement so that they execute efficiently in a multiprocessor environment, because some of the computations have to be carried out sequentially. This is due to the fact that only one new basis vector of the Krylov solution subspace is generated in each iteration. It is therefore interesting to investigate the performance of other solution methods that use a solution subspace basis that can be generated in parallel and, therefore, more efficiently on many computers. This paper proposes solution methods that use a solution subspace basis that is made up of discretized Chebyshev polynomials. It compares their performance to a Krylov subspace method that is based on partial Golub-Kahan bidiagonalization of the system matrix, and to a randomized method. The application of a solution subspace basis made up of discretized Chebyshev polynomial is found to be competitive when solving linear discrete ill-posed problems in one space-dimension and for some problems in higher space-dimensions.

Solution of linear discrete ill-posed problems by discretized Chebyshev expansion

Bai X.
;
Buccini A.;Reichel L.
2021

Abstract

Large-scale linear discrete ill-posed problems are generally solved by Krylov subspace iterative methods. However, these methods can be difficult to implement so that they execute efficiently in a multiprocessor environment, because some of the computations have to be carried out sequentially. This is due to the fact that only one new basis vector of the Krylov solution subspace is generated in each iteration. It is therefore interesting to investigate the performance of other solution methods that use a solution subspace basis that can be generated in parallel and, therefore, more efficiently on many computers. This paper proposes solution methods that use a solution subspace basis that is made up of discretized Chebyshev polynomials. It compares their performance to a Krylov subspace method that is based on partial Golub-Kahan bidiagonalization of the system matrix, and to a randomized method. The application of a solution subspace basis made up of discretized Chebyshev polynomial is found to be competitive when solving linear discrete ill-posed problems in one space-dimension and for some problems in higher space-dimensions.
978-1-6654-5843-6
Chebyshev expansion method; Linear discrete ill-posed problems; Tikhonov regularization
File in questo prodotto:
File Dimensione Formato  
Chebyshev.pdf

Solo gestori archivio

Tipologia: versione pre-print
Dimensione 480.89 kB
Formato Adobe PDF
480.89 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11584/332832
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 0
  • ???jsp.display-item.citation.isi??? 0
social impact