{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,5,5]],"date-time":"2026-05-05T05:48:39Z","timestamp":1777960119393,"version":"3.51.4"},"reference-count":36,"publisher":"Springer Science and Business Media LLC","issue":"5","license":[{"start":{"date-parts":[[2024,8,13]],"date-time":"2024-08-13T00:00:00Z","timestamp":1723507200000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2024,8,13]],"date-time":"2024-08-13T00:00:00Z","timestamp":1723507200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"funder":[{"name":"Instituto Polit\u00e9cnico de Coimbra"}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["Fract Calc Appl Anal"],"published-print":{"date-parts":[[2024,10]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:p>It is well-known that the two-parameter Mittag-Leffler (ML) function plays a key role in Fractional Calculus. In this paper, we address the problem of computing this function, when its argument is a square matrix. Effective methods for solving this problem involve the computation of higher order derivatives or require the use of mixed precision arithmetic. In this paper, we provide an alternative method that is derivative-free and works entirely using IEEE standard double precision arithmetic. If certain conditions are satisfied, our method uses a Taylor series representation for the ML function; if not, it switches to a Schur-Parlett technique that will be combined with the Cauchy integral formula. A detailed discussion on the choice of a convenient contour is included. Theoretical and numerical issues regarding the performance of the proposed algorithm are discussed. A set of numerical experiments shows that our novel approach is competitive with the state-of-the-art method for IEEE double precision arithmetic, in terms of accuracy and CPU time. For matrices whose Schur decomposition has large blocks with clustered eigenvalues, our method far outperforms the other. Since our method does not require the efficient computation of higher order derivatives, it has the additional advantage of being easily extended to other matrix functions (e.g., special functions).<\/jats:p>","DOI":"10.1007\/s13540-024-00326-9","type":"journal-article","created":{"date-parts":[[2024,8,13]],"date-time":"2024-08-13T22:01:55Z","timestamp":1723586515000},"page":"2248-2274","update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":7,"title":["Computing the Mittag-Leffler function of a matrix argument"],"prefix":"10.1007","volume":"27","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-8132-520X","authenticated-orcid":false,"given":"Jo\u00e3o R.","family":"Cardoso","sequence":"first","affiliation":[]}],"member":"297","published-online":{"date-parts":[[2024,8,13]]},"reference":[{"key":"326_CR1","first-page":"2031","volume":"236","author":"R Agarwal","year":"1953","unstructured":"Agarwal, R.: A propos d\u2019une note de M. Pierre Humbert. C. R. Acad. Sci. Paris 236, 2031\u20132032 (1953)","journal-title":"C. R. Acad. Sci. Paris"},{"key":"326_CR2","first-page":"414","volume":"39","author":"J Cardoso","year":"2012","unstructured":"Cardoso, J.: Computation of the matrix $$p$$th root and its Fr\u00e9chet derivative by integrals. Electron. Trans. Numer. Anal. 39, 414\u2013436 (2012)","journal-title":"Electron. Trans. Numer. Anal."},{"key":"326_CR3","unstructured":"Cardoso, J.; Matrix Mittag-Leffler Function https:\/\/www.mathworks.com\/matlabcentral\/fileexchange\/168776-matrix-mittag-leffler-function, MATLAB Central File Exchange. Retrieved June 25, (2024)"},{"key":"326_CR4","doi-asserted-by":"publisher","first-page":"464","DOI":"10.1137\/S0895479802410815","volume":"25","author":"P Davies","year":"2003","unstructured":"Davies, P., Higham, N.: A Schur-Parlett algorithm for computing matrix functions. SIAM J. Matrix Anal. Appl. 25, 464\u2013485 (2003)","journal-title":"SIAM J. Matrix Anal. Appl."},{"key":"326_CR5","volume-title":"Higher Transcendental Functions","author":"A Erd\u00e9lyi","year":"1955","unstructured":"Erd\u00e9lyi, A., Magnus, W., Oberhettinger, F., Tricomi, F.: Higher Transcendental Functions, vol. 3. McGraw-Hill, New York (1955)"},{"key":"326_CR6","doi-asserted-by":"publisher","first-page":"182","DOI":"10.1016\/j.laa.2019.04.001","volume":"574","author":"M Fasi","year":"2019","unstructured":"Fasi, M.: Optimality of the Paterson-Stockmeyer method for evaluating matrix polynomials and rational matrix functions. Lin. Algebr. Appl. 574, 182\u2013200 (2019)","journal-title":"Lin. Algebr. Appl."},{"issue":"5","key":"326_CR7","doi-asserted-by":"publisher","first-page":"1045","DOI":"10.1016\/j.matcom.2010.10.009","volume":"81","author":"R Garrappa","year":"2011","unstructured":"Garrappa, R., Popolizio, M.: On the use of matrix functions for fractional partial differential equations. Math. Comput. Simul. 81(5), 1045\u20131056 (2011)","journal-title":"Math. Comput. Simul."},{"key":"326_CR8","doi-asserted-by":"publisher","first-page":"135","DOI":"10.1137\/140971191","volume":"53","author":"R Garrappa","year":"2015","unstructured":"Garrappa, R.: Numerical evaluation of two and three parameter Mittag-Leffler function. SIAM J. Numer. Anal. 53, 135\u2013169 (2015)","journal-title":"SIAM J. Numer. Anal."},{"key":"326_CR9","unstructured":"Garrappa, R.: The Mittag-Leffler Function. https:\/\/www.mathworks.com\/matlabcentral\/fileexchange\/48154-the-mittag-leffler-Function. MATLAB Central File Exchange. Retrieved July 4, (2023)"},{"issue":"5","key":"326_CR10","doi-asserted-by":"publisher","first-page":"977","DOI":"10.1016\/j.camwa.2016.11.028","volume":"74","author":"R Garrappa","year":"2017","unstructured":"Garrappa, R., Moret, I., Popolizio, M.: On the time-fractional Schr\u00f6dinger equation: theoretical analysis and numerical solution by matrix Mittag-Leffler functions. Comput. Math. Appl. 74(5), 977\u2013992 (2017)","journal-title":"Comput. Math. Appl."},{"key":"326_CR11","doi-asserted-by":"publisher","first-page":"129","DOI":"10.1007\/s10915-018-0699-5","volume":"77","author":"R Garrappa","year":"2018","unstructured":"Garrappa, R., Popolizio, M.: Computing the matrix Mittag-Leffler function with applications to fractional calculus. J. Sci. Comput. 77, 129\u2013153 (2018)","journal-title":"J. Sci. Comput."},{"key":"326_CR12","doi-asserted-by":"publisher","DOI":"10.56021\/9781421407944","volume-title":"Matrix Computations","author":"G Golub","year":"2013","unstructured":"Golub, G., Van Loan, C.: Matrix Computations, 4th edn. Johns Hopkins University Press, Baltimore (2013)","edition":"4"},{"key":"326_CR13","doi-asserted-by":"crossref","unstructured":"Gorenflo, R., Kilbas, A., Mainardi, F., Rogosin, S.: Mittag-Leffler Functions, Related Topics and Applications. Springer-Verlag, 2nd ed. (2020)","DOI":"10.1007\/978-3-662-61550-8"},{"key":"326_CR14","unstructured":"Gorenflo, R., Loutchko, J., Luchko, Y.: Computation of the Mittag-Leffler function $$E_{\\alpha ,\\beta } (z)$$ and its derivatives. Fract. Calc. Appl. Anal. 5(4), 491\u2013518 (2002). Erratum: Fract. Calc. Appl. Anal. 6 (2003)"},{"key":"326_CR15","doi-asserted-by":"publisher","DOI":"10.1137\/1.9780898717778","volume-title":"Functions of Matrices: Theory and Computation","author":"N Higham","year":"2008","unstructured":"Higham, N.: Functions of Matrices: Theory and Computation. Society for Industrial and Applied Mathematics, Philadelphia (2008)"},{"issue":"3","key":"326_CR16","doi-asserted-by":"publisher","first-page":"1401","DOI":"10.1137\/20M1365326","volume":"42","author":"N Higham","year":"2021","unstructured":"Higham, N., Liu, X.: A Multiprecision derivative-free Schur-Parlett algorithm for computing matrix functions. SIAM J. Matrix Anal. Appl. 42(3), 1401\u20131422 (2021)","journal-title":"SIAM J. Matrix Anal. Appl."},{"key":"326_CR17","unstructured":"Higham, N.: The Matrix Function Toolbox. https:\/\/www.mathworks.com\/matlabcentral\/fileexchange\/20820-the-matrix-function-toolbox, MATLAB Central File Exchange. Retrieved August 7, (2023)"},{"key":"326_CR18","volume-title":"Topics in Matrix Analysis","author":"R Horn","year":"1994","unstructured":"Horn, R., Johnson, C.: Topics in Matrix Analysis. Cambridge Univ. Press, Cambridge, Paperback Edition (1994)"},{"key":"326_CR19","unstructured":"Humbert, P., Agarwal, R.: Sur la fonction de Mittag-Leffler et quelquenes de ses g\u00e9n\u00e9ralisationes. Bull. Sci. Math. (Ser. II). 77, 180\u2013185 (1953)"},{"key":"326_CR20","first-page":"937","volume":"135","author":"M Mittag-Leffler","year":"1902","unstructured":"Mittag-Leffler, M.: Sur l\u2019int\u00e9grale de Laplace-Abel. Comp. Rend. Acad. Sci. Paris 135, 937\u2013939 (1902)","journal-title":"Comp. Rend. Acad. Sci. Paris"},{"key":"326_CR21","first-page":"537","volume":"136","author":"M Mittag-Leffler","year":"1903","unstructured":"Mittag-Leffler, M.: Une g\u00e9n\u00e9ralization de l\u2019int\u00e9grale de Laplace-Abel. Comp. Rend. Acad. Sci. Paris 136, 537\u2013539 (1903)","journal-title":"Comp. Rend. Acad. Sci. Paris"},{"key":"326_CR22","first-page":"554","volume":"137","author":"M Mittag-Leffler","year":"1903","unstructured":"Mittag-Leffler, M.: Sur la nouvelle fonction $$E_\\alpha (x)$$. Comp. Rend. Acad. Sci. Paris 137, 554\u2013558 (1903)","journal-title":"Comp. Rend. Acad. Sci. Paris"},{"key":"326_CR23","unstructured":"Mittag-Leffler, M.: Sopra la funzione $$E_\\alpha (x)$$. Rend. R. Acc. Lincei, (Ser. 5) 13, 3\u20135 (1904)"},{"key":"326_CR24","doi-asserted-by":"crossref","unstructured":"Mittag-Leffler, M.: Sur la representation analytique d\u2019une branche uniforme d\u2019une fonction monog\u00e8ne (cinqui\u00e8me note) Acta Math. 29, 101\u2013181 (1905)","DOI":"10.1007\/BF02403200"},{"issue":"5","key":"326_CR25","doi-asserted-by":"publisher","first-page":"2144","DOI":"10.1137\/080738374","volume":"49","author":"I Moret","year":"2011","unstructured":"Moret, I., Novati, P.: On the convergence of Krylov subspace methods for matrix Mittag-Leffler functions. SIAM J. Numer. Anal. 49(5), 2144\u20132164 (2011)","journal-title":"SIAM J. Numer. Anal."},{"issue":"6","key":"326_CR26","doi-asserted-by":"publisher","first-page":"725","DOI":"10.1515\/ijnsns-2018-0358","volume":"20","author":"M Ortigueira","year":"2019","unstructured":"Ortigueira, M., Lopes, A., Machado, J.: On the numerical computation of the Mittag-Leffler function. Int. J. Nonlinear Sci. Numer. Simul. 20(6), 725\u2013736 (2019)","journal-title":"Int. J. Nonlinear Sci. Numer. Simul."},{"key":"326_CR27","doi-asserted-by":"publisher","first-page":"117","DOI":"10.1016\/0024-3795(76)90018-5","volume":"14","author":"B Parlett","year":"1976","unstructured":"Parlett, B.: A recurrence among the elements of functions of triangular matrices. Linear Algebr. Appl. 14, 117\u2013121 (1976)","journal-title":"Linear Algebr. Appl."},{"issue":"1","key":"326_CR28","doi-asserted-by":"publisher","first-page":"60","DOI":"10.1137\/0202007","volume":"2","author":"M Paterson","year":"1973","unstructured":"Paterson, M., Stockmeyer, L.: On the number of nonscalar multiplications necessary to evaluate polynomials. SIAM J. Comput. 2(1), 60\u201366 (1973)","journal-title":"SIAM J. Comput."},{"key":"326_CR29","unstructured":"Podlubny, I.: Fractional Differential Equations, An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications. Math. Sci. Engrg. 198, Academic Press Inc., San Diego, CA (1999)"},{"key":"326_CR30","unstructured":"Podlubny, I.: Mittag-Leffler Function. https:\/\/www.mathworks.com\/matlabcentral\/fileexchange\/8738-mittag-leffler-function, MATLAB Central File Exchange. Retrieved July 4, (2023)"},{"issue":"1","key":"326_CR31","doi-asserted-by":"publisher","first-page":"69","DOI":"10.1137\/070700280","volume":"47","author":"H Seybold","year":"2008","unstructured":"Seybold, H., Hilfer, R.: Numerical algorithm for calculating the generalized Mittag-Leffler function. SIAM J. Numer. Anal. 47(1), 69\u201388 (2008)","journal-title":"SIAM J. Numer. Anal."},{"key":"326_CR32","doi-asserted-by":"publisher","first-page":"22","DOI":"10.1016\/j.apnum.2014.11.007","volume":"90","author":"E Sousa","year":"2015","unstructured":"Sousa, E., Li, C.: A weighted finite difference method for the fractional diffusion equation based on the Riemann-Liouville derivative. Appl. Numer. Math. 90, 22\u201337 (2015)","journal-title":"Appl. Numer. Math."},{"key":"326_CR33","doi-asserted-by":"publisher","first-page":"68","DOI":"10.1007\/s10444-022-09984-w","volume":"48","author":"E Sousa","year":"2022","unstructured":"Sousa, E.: Convergence of consistent and inconsistent schemes for fractional diffusion problems with boundaries. Adv. Comput. Math. 48, 68 (2022)","journal-title":"Adv. Comput. Math."},{"key":"326_CR34","doi-asserted-by":"publisher","first-page":"558","DOI":"10.1553\/etna_vol54s558","volume":"54","author":"F Tatsuoka","year":"2021","unstructured":"Tatsuoka, F., Sogabe, T., Miyatake, Y., Kemmochi, T., Zhang, S.-L.: Computing the matrix fractional power with the double exponential formula. Electron. Trans. Numer. Anal. 54, 558\u2013580 (2021)","journal-title":"Electron. Trans. Numer. Anal."},{"issue":"3","key":"326_CR35","doi-asserted-by":"publisher","first-page":"385","DOI":"10.1137\/130932132","volume":"56","author":"L Trefethen","year":"2014","unstructured":"Trefethen, L., Weideman, J.: The exponentially convergent trapezoidal rule. SIAM Rev. 56(3), 385\u2013458 (2014)","journal-title":"SIAM Rev."},{"key":"326_CR36","doi-asserted-by":"publisher","first-page":"191","DOI":"10.1007\/BF02403202","volume":"29","author":"A Wiman","year":"1905","unstructured":"Wiman, A.: \u00dcber den fundamentalsatz der theorie der funkntionen $$E_\\alpha (x)$$. Acta Math. 29, 191\u2013201 (1905)","journal-title":"Acta Math."}],"container-title":["Fractional Calculus and Applied Analysis"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1007\/s13540-024-00326-9.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/article\/10.1007\/s13540-024-00326-9\/fulltext.html","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1007\/s13540-024-00326-9.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,9,30]],"date-time":"2024-09-30T20:35:38Z","timestamp":1727728538000},"score":1,"resource":{"primary":{"URL":"https:\/\/link.springer.com\/10.1007\/s13540-024-00326-9"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2024,8,13]]},"references-count":36,"journal-issue":{"issue":"5","published-print":{"date-parts":[[2024,10]]}},"alternative-id":["326"],"URL":"https:\/\/doi.org\/10.1007\/s13540-024-00326-9","relation":{},"ISSN":["1311-0454","1314-2224"],"issn-type":[{"value":"1311-0454","type":"print"},{"value":"1314-2224","type":"electronic"}],"subject":[],"published":{"date-parts":[[2024,8,13]]},"assertion":[{"value":"27 February 2024","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"30 July 2024","order":2,"name":"revised","label":"Revised","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"31 July 2024","order":3,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"13 August 2024","order":4,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}},{"order":1,"name":"Ethics","group":{"name":"EthicsHeading","label":"Declarations"}},{"value":"The author declares no conflict of interest.","order":2,"name":"Ethics","group":{"name":"EthicsHeading","label":"Conflicts of interest"}}]}}