Discovery of the double exponential transformation and its developments

This article is mainly concerned how the double exponential formula for numerical integration I was discovered and how it has been developed thereafter. For evaluation of I = integral(1)(-1) f(x)dx H. Takahasi and M. Mori took advantage of the optimality of the trapezoidal formula over (-infinity, infinity) with all equal mesh size and proposed in 1974 x = phi(t) = tanh(pi/2 sinh t) as all optimal choice of variable transformation which I = integral(infinity)(-infinity) f(phi(t))phi'(t)dt. If the trapezoidal formula transforms the original integral is applied to the transformed integral and the resulting infinite summation is properly truncated a quadrature formula I approximate to h Sigma(n)(k=-n) f(phi(kh))phi'(kh) is obtained. Its error is expressed as O(exp(-CN/log N)) as a function of the number N (= 2n + 1) of function evaluations. Since the integrand decays double exponentially after the transformation it is called the double exponential (DE) formula. It is also shown that the formula is optimal in the sense that there is no quadrature formula obtained by variable transformation whose error decays faster than O(exp(-CN/log N)) as N becomes large. Since the paper by Takahasi and Mori was published the DE formula has gradually come to be used widely in various fields of science and engineering. In fact we call find papers in which the DE formula is successfully used in the fields of molecular physics, fluid dynamics, statistics, civil engineering, financial engineering, in particular in the field of the boundary element method, and so on. The DE transformation has turned out to be also useful for evaluation of indefinite integrals, for solution of integral equations and for solution of ordinary differential equations, so that the scope of its applications is expected to spread also in the future.