A restarted Lanczos approximation to functions of a symmetric matrix

IMA Journal of Numerical Analysis

Volumn 30, Issue 4, 2010, Pages 1044-1061

  Ilic M ^a   Turner, I W ^a   Simpson, D P ^a  

^a QUEENSLAND UNIVERSITY OF TECHNOLOGY   (Australia)

Author keywords

error bounds; fractional Poisson equation; Gaussian Markov random fields; Krylov subspace methods; matrix functions; Stieltjes transforms

Indexed keywords

ERROR ANALYSIS; MARKOV PROCESSES;

ERROR BOUND; GAUSSIAN MARKOV RANDOM FIELD; KRYLOV SUBSPACE METHOD; MATRIX FUNCTIONS; STIELTJES TRANSFORM;

POISSON EQUATION;

EID: 77958182839     PISSN: 02724979     EISSN: 14643642     Source Type: Journal    
DOI: 10.1093/imanum/drp003     Document Type: Article

References (27)

1. BERG, C., BOYADZHIEV, K. & DELAUBENFELS, R. (1993) Generation of generators of holomorphic semigroups. J. Aust. Math. Soc. Ser. A, 55, 246-269.
2. BERG, C. & FORST, G. (1975) Potential Theory on Locally Compact Abelian Groups. Berlin: Springer.
3. BERG, C. & PEDERSEN, H. (2001) A completely monotone function related to the gamma function. J. Comput. Appl. Math., 133, 219-230.
4. BESAG, J. E. & HIGDON, D. (1999) Bayesian analysis of agricultural field experiments (with discussion). J. R. Stat. Soc. B, 61, 691-746.
5. BURDEN, R. & FAIRES, J. (2001) Numerical Analysis, 7th edn., Belmont, CA: Brooks/Cole.
6. CHENG, S. H., HIGHAM, N. J., KENNEY, C. S. & LAUB, A. J. (2001) Approximating the logarithm of a matrix to specified accuracy. SIAM J. Matrix Anal. Appl., 22, 1112-1125.
7. CRESSIE, N. (1991) Statistics For Spacial Data. Wiley Series in Probability and Mathematical Statistics. New York: John Wiley.
8. DE BOOR, C. (2005) Divided differences. Surv. Approx. Theory, 1, 46-69.
9. DRUSKIN, V. & KNIZHNERMAN, L. (1995) Krylov subspace approximations of eigenpairs and matrix functions in exact and computer arithmetic. Numer. Linear Algebra Appl., 2, 205-217.
10. EIERMANN, M. & ERNST, O. G. (2006) A restarted Krylov subspace method for the evaluation of matrix functions. SIAM J. Numer. Anal., 44, 2481-2504.
11. GILKS, W. R., RICHARDSON, S. & SPIEGELHALTER, D. J. (eds) (1998) Markov Chain Monte Carlo in Practice. New York: Chapman & Hall.
12. HIGHAM, N. J. (2008) Functions of Matrices: Theory and Computation. Philadelphia, PA: SIAM.
13. HOCHBRUCK, M. & LUBICH, C. (1997) On Krylov subspace approximations to the matrix exponential operator. SIAM J. Numer. Anal., 34, 1911-1925.
14. ILIĆ , M., LIU, F., TURNER, I. W. & ANH, V. (2005) Numerical approximation of a fractional-in-space diffusion equation I. Fract. Calc. Appl. Anal., 8, 323-341.
15. ILIĆ , M., LIU, F., TURNER, I. W. & ANH, V. (2006) Numerical approximation of a fractional-in-space diffusion equation II with nonhomogeneous boundary conditions. Fract. Calc. Appl. Anal., 9, 333-349.
16. ILIĆ , M., TURNER, I. W. & ANH, V. (2008) Numerical solution of the fractional Poisson equations using an adaptively preconditioned Lanczos method. J. Appl. Math. Stoch. Anal. doi:10.1155/2008/104525.
17. ILIĆ , M., TURNER, I. W. & PETTITT, A. N. (2004) Bayesian computations and efficient algorithms for computing functions of large, sparse matrices. ANZIAM J., 45, C504-C518.
18. LOPEZ, L. & SIMONCINI, V. (2006) Analysis of projection methods for rational function approximation to the matrix exponential. SIAM J. Numer. Anal., 44, 613-635.
19. PETTITT, A. N., WEIR, I. S. & HART, A. G. (2002) A conditional autoregressive Gaussian process for irregularly spaced multivariate data with application to modelling large sets of binary data. Stat. Comput., 12, 353-367.
20. POPOLIZIO, M. & SIMONCINI, V. (2008) Acceleration techniques for approximating the matrix exponential operator. SIAM. J. Matrix Anal. Appl., 30, 657-683.
21. RUE, H. & HELD, L. (1975) Gaussian Markov Random Fields: Theory and Applications. London: Chapman & Hall.
22. SAAD, Y. (1992) Analysis of some Krylov subspace approximations to the matrix exponential operator. SIAM J. Numer. Anal., 29, 209-228.
23. SAAD, Y. (1993) Iterative Methods for Sparse Linear Systems. Boston, MA: PWS.
24. SIDJE, R. B. (1998) Expokit: a software package for computing matrix exponentials. ACM Trans. Math. Softw., 24, 130-156.
25. SIMPSON, D. P., TURNER, I. W. & PETTITT, A. N. (2009) Fast sampling from a Gaussian Markov random field using Krylov subspace approaches. (in preparation).
26. VAN DEN ESHOF, J., FROMMER, A., LIPPERT, T., SCHILLING, K. & VAN DER VORST, H. (2002) Numerical methods for the QCD overlap operator. I. Sign-function and error bounds. Comput. Phys. Commun., 146, 203-224. (Pubitemid 34624104)
27. VAN DER VORST, H. A. (1987) An iterative solution method for solving f (A)x = b using Krylov subspace information obtained for the symmetric positive definite matrix A. J. Comput. Appl. Math., 18, 249-263.