{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,8,21]],"date-time":"2026-08-21T17:34:32Z","timestamp":1787333672262,"version":"build-2736575974"},"reference-count":37,"publisher":"Society for Industrial & Applied Mathematics (SIAM)","issue":"4","funder":[{"DOI":"10.13039\/100014600","name":"MathWorks","doi-asserted-by":"crossref","id":[{"id":"10.13039\/100014600","id-type":"DOI","asserted-by":"crossref"}]},{"DOI":"10.13039\/501100000266","name":"Engineering and Physical Sciences Research Council","doi-asserted-by":"publisher","award":["EP\/P020720\/1"],"award-info":[{"award-number":["EP\/P020720\/1"]}],"id":[{"id":"10.13039\/501100000266","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100000288","name":"Royal Society","doi-asserted-by":"publisher","id":[{"id":"10.13039\/501100000288","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["SIAM J. Matrix Anal. Appl."],"published-print":{"date-parts":[[2019,1]]},"abstract":"<jats:p>The most popular algorithms for computing the matrix exponential are those based on the scaling and squaring technique. For optimal efficiency these are usually tuned to a particular precision of floating-point arithmetic. We design a new scaling and squaring algorithm that takes the unit roundoff of the arithmetic as input and chooses the algorithmic parameters in order to keep the forward error in the underlying Pad\u00e9 approximation below the unit roundoff. To do so, we derive an explicit expression for all the coefficients in an error expansion for Pad\u00e9 approximants to the exponential and use it to obtain a new bound for the truncation error. We also derive a new technique for selecting the internal parameters used by the algorithm, which at each step decides whether to scale or to increase the degree of the approximant. The algorithm can employ diagonal Pad\u00e9 approximants or Taylor approximants and can be used with a Schur decomposition or in transformation-free form. Our numerical experiments show that the new algorithm performs in a forward stable way for a wide range of precisions and that the most accurate of our implementations, the Taylor-based transformation-free variant, is superior to existing alternatives.<\/jats:p>","DOI":"10.1137\/18m1228876","type":"journal-article","created":{"date-parts":[[2019,10,2]],"date-time":"2019-10-02T09:40:40Z","timestamp":1570009240000},"page":"1233-1256","source":"Crossref","is-referenced-by-count":26,"title":["An Arbitrary Precision Scaling and Squaring Algorithm for the Matrix Exponential"],"prefix":"10.1137","volume":"40","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-6015-391X","authenticated-orcid":true,"given":"Massimiliano","family":"Fasi","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-5956-4976","authenticated-orcid":true,"given":"Nicholas J.","family":"Higham","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"351","published-online":{"date-parts":[[2019,10,1]]},"reference":[{"key":"atypb1","doi-asserted-by":"publisher","DOI":"10.1137\/17M1145227"},{"key":"atypb2","doi-asserted-by":"publisher","DOI":"10.1137\/09074721X"},{"key":"atypb3","doi-asserted-by":"publisher","DOI":"10.1137\/110852553"},{"key":"atypb4","doi-asserted-by":"publisher","DOI":"10.1137\/140973979"},{"key":"atypb5","doi-asserted-by":"publisher","DOI":"10.1137\/16M1057577"},{"key":"atypb6","doi-asserted-by":"publisher","DOI":"10.1137\/141000671"},{"key":"atypb7","doi-asserted-by":"publisher","DOI":"10.1016\/j.cam.2018.07.042"},{"key":"atypb8","doi-asserted-by":"publisher","DOI":"10.1016\/j.apm.2015.09.050"},{"key":"atypb9","doi-asserted-by":"publisher","DOI":"10.1145\/2491491.2491494"},{"key":"atypb10","doi-asserted-by":"publisher","DOI":"10.1007\/s101070100263"},{"key":"atypb11","doi-asserted-by":"publisher","DOI":"10.1016\/j.laa.2019.04.001"},{"key":"atypb12","doi-asserted-by":"publisher","DOI":"10.1137\/17M1129866"},{"key":"atypb13","doi-asserted-by":"publisher","DOI":"10.1137\/140997610"},{"key":"atypb14","doi-asserted-by":"crossref","unstructured":"W. Gautschi,\n                      Numerical Analysis\n                      , Birkh\u00e4user, Basel, 2011.","DOI":"10.1007\/978-0-8176-8259-0"},{"key":"atypb15","doi-asserted-by":"publisher","DOI":"10.1007\/s11075-005-8141-0"},{"key":"atypb17","doi-asserted-by":"publisher","DOI":"10.1137\/04061101X"},{"key":"atypb18","doi-asserted-by":"crossref","unstructured":"N. J. Higham,\n                      Functions of Matrices: Theory and Computation\n                      , SIAM, Philadelphia, 2008,https:\/\/doi.org\/10.1137\/1.9780898717778.","DOI":"10.1137\/1.9780898717778"},{"key":"atypb19","doi-asserted-by":"publisher","DOI":"10.1137\/090768539"},{"key":"atypb20","doi-asserted-by":"publisher","DOI":"10.1137\/16M1084225"},{"key":"atypb21","doi-asserted-by":"publisher","DOI":"10.1137\/10081232X"},{"key":"atypb22","doi-asserted-by":"publisher","DOI":"10.1137\/130906118"},{"key":"atypb23","doi-asserted-by":"publisher","DOI":"10.1023\/A:1026152731904"},{"key":"atypb24","doi-asserted-by":"publisher","DOI":"10.1137\/S0895479899356080"},{"key":"atypb25","doi-asserted-by":"publisher","DOI":"10.1109\/TC.2017.2690633"},{"key":"atypb27","doi-asserted-by":"publisher","DOI":"10.1016\/j.anucene.2016.09.004"},{"key":"atypb28","first-page":"221","volume":"1898","author":"Laguerre E. N.","journal-title":"France"},{"key":"atypb29","doi-asserted-by":"publisher","DOI":"10.1137\/0704033"},{"key":"atypb32","doi-asserted-by":"publisher","DOI":"10.1007\/BF01935650"},{"key":"atypb33","doi-asserted-by":"publisher","DOI":"10.1137\/1020098"},{"key":"atypb34","doi-asserted-by":"publisher","DOI":"10.1137\/S00361445024180"},{"key":"atypb36","doi-asserted-by":"publisher","DOI":"10.1137\/18M1170352"},{"key":"atypb37","doi-asserted-by":"publisher","DOI":"10.1137\/0202007"},{"key":"atypb38","doi-asserted-by":"publisher","DOI":"10.13182\/NSE10-81"},{"key":"atypb40","doi-asserted-by":"publisher","DOI":"10.1080\/00207160.2013.791392"},{"key":"atypb42","doi-asserted-by":"publisher","DOI":"10.1137\/0714065"},{"key":"atypb43","unstructured":"S. Zhao,\n                      Matrix Exponential Approximation for Burnup Equation\n                      , M.Sc. thesis, The University of Manchester, Manchester, UK, 2017."},{"key":"atypb44","doi-asserted-by":"crossref","unstructured":"P. Zimmermann, A. Casamayou, N. Cohen, G. Connan, T. Dumont, L. Fousse, F. Maltey, M. Meulien, M. Mezzarobba, C. Pernet, N. M. Thi\u00e9ry, E. Bray, J. Cremona, M. Forets, A. Ghitza, and H. Thomas,\n                      Computational Mathematics with SageMath\n                      , SIAM, Philadelphia, 2018,https:\/\/doi.org\/10.1137\/1.9781611975468.","DOI":"10.1137\/1.9781611975468"}],"container-title":["SIAM Journal on Matrix Analysis and Applications"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/epubs.siam.org\/doi\/pdf\/10.1137\/18M1228876","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,8,21]],"date-time":"2026-08-21T17:17:07Z","timestamp":1787332627000},"score":1,"resource":{"primary":{"URL":"https:\/\/epubs.siam.org\/doi\/10.1137\/18M1228876"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2019,1]]},"references-count":37,"journal-issue":{"issue":"4","published-print":{"date-parts":[[2019,1]]}},"alternative-id":["10.1137\/18M1228876"],"URL":"https:\/\/doi.org\/10.1137\/18m1228876","relation":{},"ISSN":["0895-4798","1095-7162"],"issn-type":[{"value":"0895-4798","type":"print"},{"value":"1095-7162","type":"electronic"}],"subject":[],"published":{"date-parts":[[2019,1]]}}}