Andreas Lauschke Consulting


Convergence Acceleration with Canonical Contractions of Continued Fractions


Exponential Function



The following compares three function approximations to the exponential function. One is the well-known Taylor series (green curve), the second is a continued fraction expansion that was developed as an even contraction of a continued fraction derived from a power series using Euler's method to convert a series into an S-fraction (gray curve). This continued fraction expansion is already the fastest-converging one from among the ones presented in: Cuyt, Petersen, et al.: "Handbook of Continued Fractions for Special Functions" (Springer, 2008), which is the reason I chose it as the starting point for my own contraction method. The third (blue-dashed curve) is a canonical even contraction I have developed from the second.


Comparison: Taylor Series, Continued Fraction Expansion, and Contracted Continued Fraction Expansion of exp(z)




The original continued fraction expansion for exp(z) is (gray curve above)

original

The contracted continued fraction expansion for exp(z) I have developed from this is (blue-dashed curve above)

contracted

with
partial numerators

and
partial denominators

Using the backwards-recursion algorithm (see General page) this can be implemented in a procedural manner fairly straight-forward using a simple loop. Here in Java:

Java code

Notice that all polynomials use the Horner form (see General page) for much better efficiency.

It is now up to the individual computing platform and tweaking of compiler and runtime options to find out which one of the two continued fraction expansions is faster for a prescribed accuracy. In the case of the contracted continued fraction the terms to be computed are more complicated, which costs computing time, but fewer steps are required to attain a prescribed accuracy. This is a question that cannot be answered a-priori as compiled computing platforms are too different and the values chosen for compiler and runtime settings are too diverse to make a definitive statement here. However, already for n=0 or n=1 the convergence acceleration of the contracted continued fraction is so significant that the number of additional computations required in every step by the contracted continued fraction is not of importance. On today's computing platforms the standard arithmetic operations +,-,*, and / are "cheap" and fast, which ties in with another advantage of continued fraction expansions: they simply don't need any further arithmetic operations than these four!

For example, picking n=0, the comparison is between

Original

and

Contracted

Both yield the same accuracy, the latter is the contraction of the former, and which one is calculated faster depends entirely on the computing platform and compiler and runtime settings. Using simpler terms means having to go "deeper", and going "shallow" means using more complicated terms.

This website is powered by webMathematica.
Please visit www.webMathematica.com to find out how you can empower your website with Mathematica.

webMathematica banner

Copyright (c) 2009 Andreas Lauschke Consulting. All Rights Reserved.