{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,8,21]],"date-time":"2026-08-21T18:34:50Z","timestamp":1787337290053,"version":"build-2736575974"},"reference-count":15,"publisher":"Society for Industrial & Applied Mathematics (SIAM)","issue":"4","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["SIAM J. Sci. Comput."],"published-print":{"date-parts":[[2006,1]]},"abstract":"<jats:p>Special functions are pervasive in all fields of science. The most well\u2010known application areas are in physics, engineering, chemistry, computer science, and statistics. Because of their importance, several books and a large collection of papers have been devoted to the numerical computation of these functions. The technique for providing a floating\u2010point implementation of a function differs substantially when going from a fixed finite precision context to a finite multiprecision context. In the former, the aim is to provide an optimal mathematical model, valid on a reduced argument range and requiring as few operations as possible. Here optimal means that, in relation to the model\u2019s complexity, the truncation error is as small as it can get. The total relative error, including round\u2010off error and possible argument reduction effect, should not exceed a prescribed threshold. In a finite multiprecision context, the goal is to provide a more generic technique, from which an approximant yielding the user\u2010defined accuracy can be deduced at runtime. Hence best approximants are not an option since these models have to be recomputed every time the precision is altered and a function evaluation is requested. At the same time the generic technique should generate an approximant of as low complexity as possible. In the current approach we point out how continued fraction representations of functions can be helpful in the multiprecision context. The newly developed generic technique is based mainly on the use of sharpened a priori truncation error estimates for real continued fraction representations of a real variable, developed in section 3. As illustrated in section 4, the technique is very efficient and even quite competitive when compared to the traditional fixed precision implementations. The implementation is reliable in the sense that it allows one to return a sharp interval enclosure for the requested function evaluation, at the same cost. The paper follows a recipe style. In section 2 we gather the ingredients for the new results. In section 3 we construct or prepare, for a general function $f(x)$, a continued fraction approximant satisfying all requirements of a proper implementation. In section 4 the procedure is illustrated with results obtained for several specific functions $f(x)$.<\/jats:p>","DOI":"10.1137\/050629203","type":"journal-article","created":{"date-parts":[[2006,9,8]],"date-time":"2006-09-08T18:06:03Z","timestamp":1157738763000},"page":"1437-1462","source":"Crossref","is-referenced-by-count":9,"title":["Efficient and Reliable Multiprecision Implementation of Elementary and Special Functions"],"prefix":"10.1137","volume":"28","author":[{"given":"Annie","family":"Cuyt","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Brigitte","family":"Verdonk","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Haakon","family":"Waadeland","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"351","published-online":{"date-parts":[[2006,9,5]]},"reference":[{"key":"R1","doi-asserted-by":"crossref","unstructured":"G.A. Baker, Jr., and P. Graves\u2010Morris,\n                      Pad\u00e9 Approximants\n                      , 2nd ed., Cambridge University Press, Cambridge, UK, 1996.","DOI":"10.1017\/CBO9780511530074"},{"key":"R2","doi-asserted-by":"publisher","DOI":"10.1007\/BF01390055"},{"key":"R3","doi-asserted-by":"crossref","unstructured":"B.A. Cipra,\n                      A new testament for special functions?\n                      , SIAM News, 31 (2) (1998), p. 20.","DOI":"10.1126\/science.279.5347.31"},{"key":"R4","unstructured":"W.J. Cody and W. Waite,\n                      Software Manual for the Elementary Functions\n                      , Prentice\u2013Hall, Englewood Cliffs, NJ, 1980."},{"key":"R5","doi-asserted-by":"publisher","DOI":"10.1137\/1009002"},{"key":"R6","doi-asserted-by":"publisher","DOI":"10.1137\/0720088"},{"key":"R7","doi-asserted-by":"publisher","DOI":"10.1007\/BF02166031"},{"key":"R8","unstructured":"L. Jacobsen,\n                      Convergence acceleration and analytic continuation by means of modifications of continued fractions\n                      , in Pad\u00e9 Approximants and Continued Fractions, Kongelige Norske Videnskabers Selskab 1, H. Waadeland and H. Wallin, eds., Universitetsforlaget, Trondheim\u2010Oslo\u2010Bergen\u2010Troms\u00f8, Norway, 1983, pp. 19\u201333."},{"key":"R9","doi-asserted-by":"publisher","DOI":"10.1090\/S0025-5718-1974-0373265-5"},{"key":"R10","doi-asserted-by":"publisher","DOI":"10.1216\/rmjm\/1181069962"},{"key":"R11","unstructured":"L. Lorentzen and H. Waadeland,\n                      Continued Fractions with Applications\n                      , North\u2013Holland, Amsterdam, 1992."},{"key":"R12","unstructured":"Sun Microsystems,\n                      fdlibm\n                      , www.netlib.org\/fdlibm (1995)."},{"key":"R13","unstructured":"C. Moler,\n                      The tetragamma function and numerical craftsmanship\n                      , MATLAB News & Notes, February (2002)."},{"key":"R14","doi-asserted-by":"crossref","unstructured":"J.\u2010M. Muller,\n                      Elementary Functions. Algorithms and Implementation\n                      , Birkh\u00e4user Boston, Boston, 1997.","DOI":"10.1007\/978-1-4757-2646-6"},{"key":"R15","doi-asserted-by":"publisher","DOI":"10.1007\/BF01396057"}],"container-title":["SIAM Journal on Scientific Computing"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/epubs.siam.org\/doi\/pdf\/10.1137\/050629203","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,8,21]],"date-time":"2026-08-21T17:54:54Z","timestamp":1787334894000},"score":1,"resource":{"primary":{"URL":"https:\/\/epubs.siam.org\/doi\/10.1137\/050629203"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2006,1]]},"references-count":15,"journal-issue":{"issue":"4","published-print":{"date-parts":[[2006,1]]}},"alternative-id":["10.1137\/050629203"],"URL":"https:\/\/doi.org\/10.1137\/050629203","relation":{},"ISSN":["1064-8275","1095-7197"],"issn-type":[{"value":"1064-8275","type":"print"},{"value":"1095-7197","type":"electronic"}],"subject":[],"published":{"date-parts":[[2006,1]]}}}