×

zbMATH — the first resource for mathematics

Efficient mixed rational and polynomial approximation of matrix functions. (English) Zbl 1283.65047
Summary: This paper presents an efficient method for computing approximations for general matrix functions based on mixed rational and polynomial approximations. A method to obtain this kind of approximation from rational approximations is given, reaching the highest efficiency when transforming nondiagonal rational approximations with a higher numerator degree than the denominator degree. Then, the proposed mixed rational and polynomial approximation can be successfully applied for matrix functions which have any type of rational approximation, such as Padé, Chebyshev, etc., with maximum efficiency for higher numerator degrees than the denominator degrees. The efficiency of the mixed rational and polynomial approximation is compared with the best existing evaluating schemes for general polynomial and rational approximations, providing greater theoretical accuracy with the same cost in terms of matrix multiplications. It is well known that diagonal rational approximants are generally more accurate than the corresponding nondiagonal rational approximants which have the same computational cost. Using the proposed mixed approximation we show that the above statement is no longer true, and nondiagonal rational approximants are in fact generally more accurate than the corresponding diagonal rational approximants with the same cost.

MSC:
65F60 Numerical computation of matrix exponential and similar matrix functions
41A10 Approximation by polynomials
41A20 Approximation by rational functions
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] N.J. Higham, Functions of Matrices: Theory and Computation, Society for Industrial and Applied Mathematics, Philadelphia, PA, USA, 2008. · Zbl 1167.15001
[2] Golub, G.H.; Loan, C.V., Matrix computations, Johns hopkins studies in math. sci., (1996), The Johns Hopkins University Press
[3] Paterson, M.S.; Stockmeyer, L.J., On the number of nonscalar multiplications necessary to evaluate polynomials, SIAM J. comput., 2, 1, 60-66, (1973) · Zbl 0262.65033
[4] Sastre, J.; Ibáñez, J.; Defez, E.; Ruiz, P., Efficient orthogonal matrix polynomial based method for computing matrix exponential, Appl. math. comput., 217, 14, 6451-6463, (2011) · Zbl 1211.65052
[5] Sastre, J.; Ibáñez, J.; Defez, E.; Ruiz, P., Accurate matrix exponential computation to solve coupled differential models in engineering, Math. comput. model., 54, 1835-1840, (2011) · Zbl 1235.65042
[6] Al-Mohy, A.H.; Higham, N.J., A new scaling and squaring algorithm for the matrix exponential, SIAM J. matrix anal. appl., 31, 3, 970-989, (2009) · Zbl 1194.15021
[7] Defez, E.; Sastre, J.; Ibáñez, J.; Ruiz, P., Computing matrix functions solving coupled differential models, Math. comput. model., 50, 5-6, 831-839, (2009) · Zbl 1185.65078
[8] Defez, E.; Jódar, L., Some applications of Hermite matrix polynomials series expansions, J. comput. appl. math., 99, 105-117, (1998) · Zbl 0929.33006
[9] Moler, C.B.; Loan, C.V., Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later, SIAM rev., 45, 3-49, (2003) · Zbl 1030.65029
[10] S. Blackford, J. Dongarra, Installation guide for LAPACK, LAPACK Working Note 411, Department of Computer Science University of Tenessee, 1999.
[11] Higham, N.J., The scaling and squaring method for the matrix exponential revisited, SIAM J. matrix anal. appl., 26, 4, 1179-1193, (2005) · Zbl 1081.65037
[12] Calvetti, D.; Gallopoulos, E.; Reichel, L., Incomplete partial fractions for parallel evaluation of rational matrix functions, J. comput. appl. math., 59, 349-380, (1995) · Zbl 0839.65054
[13] Calvetti, D.; Reichel, L., On the evaluation of polynomial coefficients, Numer. alg., 33, 153-161, (2003) · Zbl 1035.65156
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.