# Efficient and Stable Arnoldi Restarts for Matrix Functions Based on Quadrature

## Abstract

### Keywords

### MSC codes

## Get full access to this article

View all available purchase options and get full access to this article.

## References

*A generalization of the steepest descent method for matrix functions*, Electron. Trans. Numer. Anal., 28 (2008), pp. 206--222.

*Implementation of a restarted Krylov subspace method for the evaluation of matrix functions*, Linear Algebra Appl., 429 (2008), pp. 229--314.

*Padé Approximants*, Cambridge University Press, Cambridge, UK, 1996.

*Error estimation and evaluation of matrix functions via the Faber transform*, SIAM J. Numer. Anal., 47 (2009), pp. 3849--3883.

*Potential Theory on Locally Compact Abelian Groups*, Springer-Verlag, Heidelberg, 1975.

*An iterative method to compute the sign function of a non-Hermitian matrix and its application to the overlap Dirac operator at nonzero chemical potential*, Comput. Phys. Comm., 177 (2007), pp. 933--943.

*An efficient implicit FEM scheme for fractional-in-space reaction--diffusion equations*, SIAM J. Sci. Comput., 34 (2012), pp. A2145--A2172.

*Computation of the matrix $p$th root and its Fréchet derivative by integrals*, Electron. Trans. Numer. Anal., 39 (2012), pp. 414--436.

*A Schur--Parlett algorithm for computing matrix functions*, SIAM J. Matrix Anal. Appl., 25 (2003), pp. 464--485.

*Methods of Numerical Integration*, Academic Press, New York, 1975.

*Divided differences*, Surv. Approx. Theory, 1 (2005), pp. 46--69.

*Two polynomial methods of calculating functions of symmetric matrices*, Comput. Math. Math. Physi., 29 (1989), pp. 112--121.

*Extended Krylov subspaces: Approximation of the matrix square root and related functions*, SIAM J. Matrix Anal. Appl., 19 (1998), pp. 775--771.

*Spectral deferred correction methods for ordinary differential equations*, BIT, 40 (2000), pp. 241--266.

*A restarted Krylov subspace method for the evaluation of matrix functions*, SIAM J. Numer. Anal., 44 (2006), pp. 2481--2504.

*Deflated restarting for matrix functions*, SIAM J. Matrix Anal. Appl., 32 (2011), pp. 621--641.

*Computing Functions of Matrices Using Krylov Subspace Methods*, Technical report, Chalmers University of Technology, Gothenburg, Sweden, 1990.

*FUNM_QUAD: An Implementation of a Stable, Quadrature-Based Restarted Arnoldi Method for Matrix Functions*, Technical report, IMACM Preprint 14/04, 2014.

*Marix functions*, in Model Order Reduction: Theory, Research Aspects and Applications, H. A. van der Vorst, W. H. A. Schilders, and J. Rommes, eds., Mathematics in Industry, Springer-Verlag, Berlin, 2008.

*Efficient solution of parabolic equations by Krylov approximation methods*, SIAM J. Sci. Statist. Comput., 13 (1992), pp. 1236--1264.

*Quadrature formulae on half-infinite intervals*, BIT, 31 (1991), pp. 437--446.

*Rational Krylov Methods for Operator Functions*, Ph.D. thesis, Fakultät für Mathematik und Informatik der Technischen Universität Bergakademie Freiberg, 2010.

*Rational Krylov approximation of matrix functions: Numerical methods and optimal pole selection*, GAMM-Mitt., 36 (2013), pp. 8--31.

*A black-box rational Arnoldi variant for Cauchy--Stieltjes matrix functions*, BIT, 24 (2013), pp. 595--616.

*Applied and Computational Complex Analysis, vol. 2*, John Wiley & Sons, New York, 1977.

*Functions of Matrices: Theory and Computation*, SIAM, Philadelphia, 2008.

*Computing matrix functions*, Acta Numer., 19 (2010), pp. 159--208.

*On Krylov subspace approximations to the matrix exponential operator*, SIAM J. Numer. Anal., 34 (1997), pp. 1911--1925.

*Exponential integrators*, Acta Numer., 19 (2010), pp. 209--286.

*A restarted Lanczos approximation to functions of a symmetric matrix*, IMA J. Numer. Anal., 30 (2010), pp. 1044--1061.

*Calculation of functions of unsymmetric matrices using Arnoldi's method*, Comput. Math. Math. Phys., 31 (1991), pp. 1--9.

*A new investigation of the extended Krylov subspace method for matrix function evaluations*, Numer. Linear Algebra Appl., 17 (2010), pp. 615--638.

*An iteration method for the solution of the eigenvalue problem of linear differential and integral operators*, J. Res. National Bureau of Standards, 45 (1950), pp. 255--282.

*Infinite series for logarithm of matrix, applied to identification of linear continuous-time multivariable systems from discrete-time models*, Electron. Lett., 27 (1991), pp. 1468--1470.

*RD-rational approximations of the matrix exponential.*, BIT, 44 (2004), pp. 595--615.

*Exactly massless quarks on the lattice*, Phys. Lett. B, 417 (1998), pp. 141--144.

*Approximate solutions and eigenvalue bounds from Krylov subspaces*, Numer. Linear Algebra Appl., 1 (1993), pp. 1--7.

*Analysis of some Krylov subspace approximations to the exponential operator*, SIAM J. Numer. Anal., 29 (1992), pp. 209--228.

*Iterative Methods for Sparse Linear Systems, 2nd ed.*, SIAM, Philadelphia, 2003.

*GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems*, SIAM J. Sci. Statist. Comput., 7 (1986), pp. 856--869.

*Vectorized adaptive quadrature in Matlab*, J. Comput. Appl. Math., 211 (2008), pp. 131--140.

*Elementary Real and Complex Analysis*, MIT Press, Cambridge, MA, 1973.

*The representation of social processes by Markov models*, Amer. J. Sociol., 82 (1976), pp. 1--54.

*On restart and error estimation for Krylov approximation of $w=f(A)v$*, SIAM J. Sci. Comput., 29 (2007), pp. 2426--2441.

*Talbot quadratures and rational approximations*, BIT, 46 (2006), pp. 653--670.

*Numerical methods for the QCD overlap operator, I. Sign-function and error bounds*, Comput. Phys. Comm., 146 (2002), pp. 203--224.

*Preconditioning Lanczos approximations to the matrix exponential*, SIAM J. Sci. Comput., 27 (2006), pp. 1438--1457.

*Interpolation and Approximation by Rational Functions in the Complex Domain*, AMS, Providence, RI, 1969.

*Optimizing Talbot's contours for the inversion of the Laplace transform*, SIAM J. Numer. Anal., 44 (2006), pp. 2342--2362.

*Parabolic and hyperbolic contours for computing the Bromwich integral*, Math. Comp., 76 (2007), pp. 1341--1356.

*Application of elliptic functions to the question of functions deviating least and most from zero*, Zap. Imp. Akad. Nauk. St. Petersburg, 30 (1877), pp. 1--59.

## Information & Authors

### Information

#### Published In

#### Copyright

#### History

**Submitted**: 29 August 2013

**Accepted**: 12 February 2014

**Published online**: 22 May 2014

#### Keywords

#### MSC codes

### Authors

## Metrics & Citations

### Metrics

### Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

#### Cited By

- Randomized Sketching for Krylov Approximations of Large-Scale Matrix FunctionsSIAM Journal on Matrix Analysis and Applications, Vol. 44, No. 3 | 25 July 2023
- Low‐rank updates of matrix square rootsNumerical Linear Algebra with Applications, Vol. 22 | 9 August 2023
- Krylov Subspace Residual and Restarting for Certain Second Order Differential EquationsSIAM Journal on Scientific Computing, Vol. 0, No. 0 | 19 July 2023
- Low-Memory Krylov Subspace Methods for Optimal Rational Matrix Function ApproximationSIAM Journal on Matrix Analysis and Applications, Vol. 44, No. 2 | 19 May 2023
- Krylov Subspace Restarting for Matrix Laplace TransformsSIAM Journal on Matrix Analysis and Applications, Vol. 44, No. 2 | 19 May 2023
- Coarse Grid Corrections in Krylov Subspace Evaluations of the Matrix ExponentialSIAM Journal on Scientific Computing, Vol. 45, No. 3 | 27 June 2023
- Integral representations for higher-order Fréchet derivatives of matrix functions: Quadrature algorithms and new results on the level-2 condition numberLinear Algebra and its Applications, Vol. 656 | 1 Jan 2023
- Error Bounds for Lanczos-Based Matrix Function ApproximationSIAM Journal on Matrix Analysis and Applications, Vol. 43, No. 2 | 12 May 2022
- Some methods for solving equations with an operator function and applications for problems with a fractional power of an operatorJournal of Computational and Applied Mathematics, Vol. 407 | 1 Jun 2022
- A study of defect-based error estimates for the Krylov approximation of φ-functionsNumerical Algorithms, Vol. 90, No. 1 | 8 November 2021
- Block Gram-Schmidt algorithms and their stability propertiesLinear Algebra and its Applications, Vol. 638 | 1 Apr 2022
- Factorized Schemes for First and Second Order Evolution Equations with Fractional Powers of OperatorsComputational Methods in Applied Mathematics, Vol. 22, No. 2 | 20 October 2021
- 3-D Large-Scale TEM Modeling Using Restarting Polynomial Krylov MethodIEEE Transactions on Geoscience and Remote Sensing, Vol. 60 | 1 Jan 2022
- An approximate representation of a solution to fractional elliptical BVP via solution of parabolic IVPJournal of Computational and Applied Mathematics, Vol. 391 | 1 Aug 2021
- Splitting schemes for non-stationary problems with a rational approximation for fractional powers of the operatorApplied Numerical Mathematics, Vol. 165 | 1 Jul 2021
- An Accurate Restarting for Shift-and-Invert Krylov Subspaces Computing Matrix Exponential Actions of Nonsymmetric MatricesComputational Mathematics and Mathematical Physics, Vol. 61, No. 5 | 1 July 2021
- Rational Krylov for Stieltjes matrix functions: convergence and pole selectionBIT Numerical Mathematics, Vol. 61, No. 1 | 4 August 2020
- A Comparison of Limited-memory Krylov Methods for Stieltjes Functions of Hermitian MatricesSIAM Journal on Matrix Analysis and Applications, Vol. 42, No. 1 | 21 January 2021
- Residual and Restarting in Krylov Subspace Evaluation of the $\varphi$ FunctionSIAM Journal on Scientific Computing, Vol. 43, No. 6 | 4 November 2021
- Fast Rational Lanczos Method for the Toeplitz Symmetric Positive Semidefinite Matrix FunctionsSimulation Tools and Techniques | 27 April 2021
- Exponential Time Integrators for Unsteady Advection–Diffusion Problems on Refined MeshesNumerical Geometry, Grid Generation and Scientific Computing | 7 May 2021
- Approximation of the matrix exponential for matrices with a skinny field of valuesBIT Numerical Mathematics, Vol. 60, No. 4 | 15 May 2020
- Limited‐memory polynomial methods for large‐scale matrix functionsGAMM-Mitteilungen, Vol. 43, No. 3 | 10 September 2020
- High order numerical schemes for solving fractional powers of elliptic operatorsJournal of Computational and Applied Mathematics, Vol. 372 | 1 Jul 2020
- Approximation of a fractional power of an elliptic operatorNumerical Linear Algebra with Applications, Vol. 27, No. 3 | 31 January 2020
- The tensor t‐function: A definition for functions of third‐order tensorsNumerical Linear Algebra with Applications, Vol. 27, No. 3 | 20 February 2020
- Computing performability measures in Markov chains by means of matrix functionsJournal of Computational and Applied Mathematics, Vol. 368 | 1 Apr 2020
- Computable upper error bounds for Krylov approximations to matrix exponentials and associated $${\varvec{\varphi }}$$-functionsBIT Numerical Mathematics, Vol. 60, No. 1 | 11 September 2019
- Block Krylov Subspace Methods for Functions of Matrices II: Modified Block FOMSIAM Journal on Matrix Analysis and Applications, Vol. 41, No. 2 | 20 May 2020
- ART: Adaptive residual-time restarting for Krylov subspace matrix exponential evaluationsJournal of Computational and Applied Mathematics, Vol. 364 | 1 Jan 2020
- Comparison of exponential and implicit time integrators for unsteady advection-diffusion problems on refined meshesKeldysh Institute Preprints, No. 119 | 1 January 2020
- Inexact Arnoldi residual estimates and decay properties for functions of non-Hermitian matricesBIT Numerical Mathematics, Vol. 59, No. 4 | 15 June 2019
- Rational approximations to fractional powers of self-adjoint positive operatorsNumerische Mathematik, Vol. 143, No. 1 | 11 May 2019
- A block rational Krylov method for 3-D time-domain marine controlled-source electromagnetic modellingGeophysical Journal International, Vol. 218, No. 1 | 8 March 2019
- Multigrid approach in shifted linear systems for the non-degenerated twisted mass operatorComputer Physics Communications, Vol. 236 | 1 Mar 2019
- Adaptive residual-time restarting for Krylov subspace matrix exponential evaluationsKeldysh Institute Preprints, No. 127 | 1 January 2019
- Exponential Krylov time integration for modeling multi-frequency optical response with monochromatic sourcesJournal of Computational and Applied Mathematics, Vol. 340 | 1 Oct 2018
- Efficient Implementation of Rational Approximations to Fractional Differential OperatorsJournal of Scientific Computing, Vol. 76, No. 1 | 22 December 2017
- Enhanced Dependability Evaluation Through Krylov Methods and Matrix Functions: The Case of Load-Sharing Systems2018 48th Annual IEEE/IFIP International Conference on Dependable Systems and Networks Workshops (DSN-W) | 1 Jun 2018
- A shifted block FOM algorithm with deflated restarting for matrix exponential computationsApplied Numerical Mathematics, Vol. 127 | 1 May 2018
- Scabies in residential care homes: Modelling, inference and interventions for well-connected population sub-unitsPLOS Computational Biology, Vol. 14, No. 3 | 26 March 2018
- Computing the Weighted Geometric Mean of Two Large-Scale Matrices and Its Inverse Times a VectorSIAM Journal on Matrix Analysis and Applications, Vol. 39, No. 1 | 1 February 2018
- Numerical Solution of Time-Dependent Problems with Fractional Power Elliptic OperatorComputational Methods in Applied Mathematics, Vol. 18, No. 1 | 17 August 2017
- Acceleration of contour integration techniques by rational Krylov subspace methodsJournal of Computational and Applied Mathematics, Vol. 316 | 1 May 2017
- Rational Approximation to the Fractional Laplacian Operator in Reaction-Diffusion ProblemsSIAM Journal on Scientific Computing, Vol. 39, No. 1 | 17 January 2017
- The Radau--Lanczos Method for Matrix FunctionsSIAM Journal on Matrix Analysis and Applications, Vol. 38, No. 3 | 1 August 2017
- Handling Fermions on the LatticeLattice Quantum Chromodynamics | 21 October 2016
- Approximation of functions of large matrices with Kronecker structureNumerische Mathematik, Vol. 135, No. 1 | 22 March 2016
- APPROXIMATION OF THE LINEAR COMBINATION OF <i>φ</i>-FUNCTIONS USING THE BLOCK SHIFT-AND-INVERT KRYLOV SUBSPACE METHODJournal of Applied Analysis & Computation, Vol. 7, No. 4 | 1 Jan 2017
- Error bounds and estimates for Krylov subspace approximations of Stieltjes matrix functionsBIT Numerical Mathematics, Vol. 56, No. 3 | 17 December 2015
- Block Krylov subspace methods for approximating the linear combination of φ -functions arising in exponential integratorsComputers & Mathematics with Applications, Vol. 72, No. 4 | 1 Aug 2016
- Parallel matrix function evaluation via initial value ODE modelingComputers & Mathematics with Applications, Vol. 72, No. 1 | 1 Jul 2016
- A framework of the harmonic Arnoldi method for evaluating φ-functions with applications to exponential integratorsAdvances in Computational Mathematics, Vol. 42, No. 3 | 1 December 2015
- Matching exponential-based and resolvent-based centrality measuresJournal of Complex Networks, Vol. 4, No. 2 | 29 June 2015
- Multigrid preconditioning for the overlap operator in lattice QCDNumerische Mathematik, Vol. 132, No. 3 | 27 May 2015
- Estimating the Largest Elements of a MatrixSIAM Journal on Scientific Computing, Vol. 38, No. 5 | 19 October 2016
- Convergence of Restarted Krylov Subspace Methods for Stieltjes Functions of MatricesSIAM Journal on Matrix Analysis and Applications, Vol. 35, No. 4 | 18 December 2014