Algorithm 894: On a block Schur - Parlett algorithm for φ-functions based on the sep-inverse estimate

ACM Transactions on Mathematical Software

Volumn 36, Issue 2, 2009, Pages

  Koikari, Souji ^a  

^a KYOTO UNIVERSITY   (Japan)

Author keywords

functions; Block Schur Parlett algorithm; Exponential integrators; Matrix functions

Indexed keywords

BLOCK SCHUR - PARLETT ALGORITHM; EXPONENTIAL INTEGRATORS; FORTRAN 95; INVERSE FUNCTIONS; MATRIX; MATRIX FUNCTIONS; NUMERICAL EXPERIMENTS; SCALING AND SQUARING METHOD;

ALGORITHMS; FUNCTION EVALUATION; SUBROUTINES;

RECURSIVE FUNCTIONS;

EID: 65849282797     PISSN: 00983500     EISSN: 15577295     Source Type: Journal    
DOI: 10.1145/1499096.1499101     Document Type: Article

References (26)

1. ANDERSON, E., BAI, Z., BISHOF, C. H., BLACKFORD, S., DEMMEL, J. W., DONGARRA, J. J., DU CROZ, J. J., GREENBAUM, A., HAMMARLING, S. J., MCKENNEY, A., AND SORENSEN, D. C. 1999. LAPACK Users' Guide, 3rd ed. Society for Industrial and Applied Mathematics, Philadelphia, PA.
2. BAI, Z. AND DEMMEL, J. W. 1993. On swapping diagonal blocks in real Schur form. Lin. Alg. Appl. 186, 73-95.
3. BARTELS, R. H. AND STEWART, G. W. 1972. Algorithm 432: Solution of the matrix equation AX + X B = C. Commun. ACM 15, 9, 820-826.
4. BERLAND, H., SKAFLESTAD, B., AND WRIGHT, W. M. 2007. EXPINT - a MATLAB package for exponential integrators. ACM Trans. Math. Softw. 33, 1, Article 4. (Pubitemid 46494536)
5. BLACKFORD, L. S., DEMMEL, J., DONGARRA, J., DUFF, I., HAMMARLING, S., HENRY, G., HEROUX, M., KAUFMAN, L., LUMSDAINE, A., PETITET, A., POZO, R., REMINGTON, K., AND WHALEY, R. C. 2002. An updated set of basic linear algebra subprograms (BLAS). ACM Trans. Math. Softw. 28, 2, 135-151. (Pubitemid 135701673)
6. DAVIES, P. I. AND HIGHAM, N. J. 2003. A Schur-Parlett algorithm for computing matrix functions. SIAM J. Matrix Anal. Appl. 25, 2, 464-485.
7. FRANK, W. L. 1958. Computing eigenvalues of complex matrices by determinant evaluation and by methods of Danilewski and Wielandt. J. Soc. Indust. Appl. Math. 6, 4, 378-392.
8. HIGHAM, N. J. 1988. FORTRAN codes for estimating the one-norm of a real or complex matrix, with applications to condition estimation. ACM Trans. Math. Softw. 14, 4, 381-396.
9. HIGHAM, N. J. 2002. Accuracy and Stability of Numerical Algorithm, 2nd ed. Society for Industrial and Applied Mathematics, Philadelphia, PA.
10. HIGHAM, N. J. 2005. The scaling and squaring method for the matrix exponential revisited. SIAM J. Matrix Anal. Appl. 26, 4, 1179-1193.
11. HIGHAM, N. J. AND TISSEUR, F. 2000. A block algorithm for matrix 1-norm estimation with an application to 1-norm pseudospectra. SIAM J. Matrix Anal. Appl. 21, 4, 1185-1201. (Pubitemid 33217064)
12. INTEL CORPORATION. 2008. Intel Fortran compiler. http://developer.intel. com/software/products/.
13. JONSSON, I. AND KåGSTRöM, B. 2002. Recursive blocked algorithms for solving triangular systems. I. One-sided and coupled Sylvester-type matrix equations. ACM Trans. Mathe. Softw. 28, (4), 392-415. (Pubitemid 135688491)
14. KåAGSTRöM, B. AND POROMAA, P. 1992. Distributed and shared memory block algorithms for the triangular Sylvester equations with sep-1 estimators. SIAM J. Matrix Anal. Appl. 13, 1, 90-101.
15. KåGSTRöM, B. AND WESTIN, L. 1989. Generalized Schur methods with condition estimators for solving the generalized Sylvester equation. IEEE Trans. Autom. Cont. 34, 745-751.
16. KASSAM, A. K. AND TREFETHEN, L. N. 2005. Fourth-order time stepping for stiff PDEs. SIAM J. Sci. Comput. 26, 4, 1214-1233.
17. KOIKARI, S. 2007. An error analysis of the modified scaling and squaring method. Comput. Math. Appl. 53, 1293-1305. (Pubitemid 46755552)
18. MINCHEV, B. AND WRIGHT, W. M. 2005. A review of exponential integrators for first order semilinear problem. Tech. rep. No. 2/2005. The Norwegian University of Science and Technology. Trondhein, Norway.
19. MOLER, C. B. AND VAN LOAN, C. F. 2003. Nineteen dubious ways to compute the exponential of a matrix, twenty five years later. SIAM Rev. 45, 1, 3-49.
20. PARLETT, B. N. 1976. A recurrence among the elements of functions of triangular matrices. Lin. Alg. Appl. 14, 117-121.
21. SIDJE, R.B. 1998. EXPOKIT: A software package for computing matrix exponentials. ACM Trans. Math. Softw. 24, 1, 130-156. (Pubitemid 128447702)
22. SUN MICROSYSTEMS. 2008. Sun Fortran 95 compiler. http://developers.sun. com/sunstudio/.
23. TOMMILA, M. 2005. apfloat, a C++ high performance arbitrary precision arithmetic package, Version 2.41. http://www.apfloat.org/.
24. TORCZON, V. 1991. On the convergence of the multidirectional search algorithm. SIAM J. Optimiz. 1, 1, 123-145.
25. VAUGHT, A. 2006. G95 Manual. http://www.g95.org/.
26. WARD, R. C. 1977. Numerical computation of the matrix exponential with accuracy estimate. SIAM J. Numer. Anal. 14, 4, 600-610.