Convergence rates for inverse-free rational approximation of matrix functions
dc.contributor.author | Jagels, Carl | |
dc.contributor.author | Mach, Thomas | |
dc.contributor.author | Reichel, Lothar | |
dc.contributor.author | Vandebril, Raf | |
dc.date.accessioned | 2017-01-09T03:33:16Z | |
dc.date.available | 2017-01-09T03:33:16Z | |
dc.date.issued | 2016-12-01 | |
dc.description.abstract | This article deduces geometric convergence rates for approximating matrix functions via inverse-free rational Krylov methods. In applications one frequently encounters matrix functions such as the matrix exponential or matrix logarithm; often the matrix under consideration is too large to compute the matrix function directly, and Krylov subspace methods are used to determine a reduced problem. If many evaluations of a matrix function of the form f(A)v with a large matrix A are required, then it may be advantageous to determine a reduced problem using rational Krylov subspaces. These methods may give more accurate approximations of f(A)v with subspaces of smaller dimension than standard Krylov subspace methods. Unfortunately, the system solves required to construct an orthogonal basis for a rational Krylov subspace may create numerical difficulties and/or require excessive computing time. This paper investigates a novel approach to determine an orthogonal basis of an approximation of a rational Krylov subspace of (small) dimension from a standard orthogonal Krylov subspace basis of larger dimension. The approximation error will depend on properties of the matrix A and on the dimension of the original standard Krylov subspace. We show that our inverse-free method for approximating the rational Krylov subspace converges geometrically (for increasing dimension of the standard Krylov subspace) to a rational Krylov subspace. The convergence rate may be used to predict the dimension of the standard Krylov subspace necessary to obtain a certain accuracy in the approximation. Computed examples illustrate the theory developed. | ru_RU |
dc.identifier.citation | Jagels, C., Mach, T., Reichel, L., & Vandebril, R. (2016). Convergence rates for inverse-free rational approximation of matrix functions. Linear Algebra and Its Applications, 510, 291-310. DOI: 10.1016/j.laa.2016.08.029 | ru_RU |
dc.identifier.uri | http://nur.nu.edu.kz/handle/123456789/2202 | |
dc.language.iso | en | ru_RU |
dc.publisher | Linear Algebra and Its Applications | ru_RU |
dc.rights | Attribution-NonCommercial-ShareAlike 3.0 United States | * |
dc.rights.uri | http://creativecommons.org/licenses/by-nc-sa/3.0/us/ | * |
dc.subject | approximation | ru_RU |
dc.subject | convergence rate | ru_RU |
dc.subject | iterative method | ru_RU |
dc.subject | matrix function | ru_RU |
dc.subject | rational Krylov | ru_RU |
dc.title | Convergence rates for inverse-free rational approximation of matrix functions | ru_RU |
dc.type | Article | ru_RU |
Files
Original bundle
1 - 1 of 1
Loading...
- Name:
- Carl Jagels, Thomas Mach, Lothar Reichel, Raf Vandebril.pdf
- Size:
- 16.76 KB
- Format:
- Adobe Portable Document Format
- Description:
License bundle
1 - 1 of 1
No Thumbnail Available
- Name:
- license.txt
- Size:
- 6.22 KB
- Format:
- Item-specific license agreed upon to submission
- Description: