An efficient algorithm for regularization of Laplace transform inversion in real case

We address design of a numerical algorithm for solving the linear system arising in numerical inversion of Laplace transforms in real case [L. D’Amore, A. Murli, Regularization of a Fourier series method for the Laplace transform inversion with real data, Inverse Problems 18 (2002) 1185–1205]. The matrix has a condition number that grows almost exponentially and the singular values decay gradually towards zero. In such a case, because of this intrinsic strong instability, the main difficulty of any numerical computation is the ability of discovering at run time, only using data, what is the maximum attainable accuracy on the solution.

In this paper, we use GMRES with the aim of relating the current residuals to the maximum attainable accuracy of the approximate solution by using a suitable stopping rule. We prove that GMRES stops after, at most, as many iterations as the number of the largest eigenvalues (compared to the machine epsilon). We use a split preconditioner that symmetrically precondition the initial system. By this way, the largest eigenvalue dynamically provides the estimate of the condition number of the matrix.