liu.seSearch for publications in DiVA
Change search
ReferencesLink to record
Permanent link

Direct link
A numerical solution of a Cauchy problem for an elliptic equation by Krylov subspaces
Linköping University, Department of Mathematics, Scientific Computing. Linköping University, The Institute of Technology.ORCID iD: 0000-0003-2281-856X
University of Bologna.
2009 (English)In: INVERSE PROBLEMS, ISSN 0266-5611 , Vol. 25, no 6, 065002- p.Article in journal (Refereed) Published
Abstract [en]

We study the numerical solution of a Cauchy problem for a self-adjoint elliptic partial differential equation u(zz) - L-u = 0 in three space dimensions (x, y, z), where the domain is cylindrical in z. Cauchy data are given on the lower boundary and the boundary values on the upper boundary are sought. The problem is severely ill-posed. The formal solution is written as a hyperbolic cosine function in terms of the two-dimensional elliptic operator L (via its eigenfunction expansion), and it is shown that the solution is stabilized (regularized) if the large eigenvalues are cut off. We suggest a numerical procedure based on the rational Krylov method, where the solution is projected onto a subspace generated using the operator L-1. This means that in each Krylov step, a well-posed two-dimensional elliptic problem involving L is solved. Furthermore, the hyperbolic cosine is evaluated explicitly only for a small symmetric matrix. A stopping criterion for the Krylov recursion is suggested based on the relative change of an approximate residual, which can be computed very cheaply. Two numerical examples are given that demonstrate the accuracy of the method and the efficiency of the stopping criterion.

Place, publisher, year, edition, pages
2009. Vol. 25, no 6, 065002- p.
National Category
URN: urn:nbn:se:liu:diva-18557DOI: 10.1088/0266-5611/25/6/065002OAI: diva2:220589
Original Publication: Lars Eldén and Valeria Simoncini, A numerical solution of a Cauchy problem for an elliptic equation by Krylov subspaces, 2009, INVERSE PROBLEMS, (25), 6, 065002. Copyright: Iop Publishing Ltd Available from: 2009-06-10 Created: 2009-06-01 Last updated: 2013-08-30Bibliographically approved

Open Access in DiVA

fulltext(3119 kB)634 downloads
File information
File name FULLTEXT01.pdfFile size 3119 kBChecksum SHA-512
Type fulltextMimetype application/pdf

Other links

Publisher's full text

Search in DiVA

By author/editor
Eldén, Lars
By organisation
Scientific ComputingThe Institute of Technology

Search outside of DiVA

GoogleGoogle Scholar
Total: 634 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Altmetric score

Total: 480 hits
ReferencesLink to record
Permanent link

Direct link