<<

. 34
( 83 .)



>>

Many further numerical examples are given in the literature [39,41“ 44,50,84].


Appendix A. Stieltjes series and functions

A Stieltjes series is a formal expansion

(’1) j j z j
f(z) = (A.1)
j=0

with partial sums
n
(’1) j j z j :
fn (z) = (A.2)
j=0

The coe cients n are the moments of an uniquely given positive measure (t) that has inÿnitely
many di erent values on 06t ¡ ∞ [4, p. 159]:

t n d (t);
= n ∈ N0 : (A.3)
n
0
142 H.H.H. Homeier / Journal of Computational and Applied Mathematics 122 (2000) 81“147


Formally, the Stieltjes series can be identiÿed with a Stieltjes integral

d (t)
f(z) = ; |arg(z)| ¡ : (A.4)
1 + zt
0

If such an integral exists for a function f then the function is called a Stieltjes function. For every
Stieltjes function there exist a unique asymptotical Stieltjes series (A.1), uniformly in every sector
|arg(z)| ¡ ‚ for all ‚ ¡ . For any Stieltjes series, however, several di erent corresponding Stieltjes
functions may exist. To ensure uniqueness, additional criteria are necessary [88, Section 4.3].
In the context of convergence acceleration and summation of divergent series, it is important that
for given z the tails f(z) ’ fn (z) of a Stieltjes series are bounded in absolute value by the next
term of the series,
n+1
|f(z) ’ fn (z)|6 n+1 z z¿0: (A.5)
Hence, for Stieltjes series the remainder estimates may be chosen as
!n = (’1)n+1 n+1
n+1 z : (A.6)
fn (z), i.e., to a t˜ variant.
This corresponds to !n =


Appendix B. Derivation of the recursive scheme (148)
(k) (k)
We show that for the divided di erence operator = n [{{x n }}] the identity
n
(k) (k)
(x n+k+1 + ˜) n+1 ((x)˜ g(x)) ’ (x n + ˜) n ((x)˜ g(x))
(k+1)
((x)˜+1 g(x)) = (B.1)
n
x n+k+1 ’ x n
holds. The proof is based on the Leibniz formula for divided di erences (see, e.g., [69, p. 50]) that
(k)
yields upon use of (x)˜+1 = (x + ˜)(x)˜ and n (x) = x n k; 0 + k; 1
k+1
(k+1’j)
(k+1) (k+1) ( j)
((x)˜+1 g(x)) = ˜ n ((x)˜ g(x)) + n (x) n+j ((x)˜ g(x))
n
j=0
(k)
(k+1)
= (x n + ˜) ((x)˜ g(x)) + n+1 ((x)˜ g(x)): (B.2)
n

Using the recursion relation of the divided di erences, one obtains
(k) (k)
n+1 ((x)˜ g(x)) ’ n ((x)˜ g(x)) (k)
(k+1)
((x)˜+1 g(x)) = (x n + ˜) + n+1 ((x)˜ g(x)): (B.3)
n
x n+k+1 ’ x n
Simple algebra then yields Eq. (B.1).
Comparison with Eq. (140) shows that using the interpolation conditions gn = g(x n ) = sn =!n and
˜ = k ’ 1 in Eq. (B.1) yields the recursion for the numerators in Eq. (148), while the recursion
for the denominators in Eq. (148) follows for ˜ = k ’ 1 and using the interpolation conditions
gn = g(x n ) = 1=!n . In each case, the initial conditions follow directly from Eq. (140) in combination
with the deÿnition of the divided di erence operator: For k = 0, we use (a)’1 = 1=(a ’ 1) and obtain
(k)
n (x n )k’1 gn = (x n )’1 gn = gn =(x n ’ 1).
H.H.H. Homeier / Journal of Computational and Applied Mathematics 122 (2000) 81“147 143


Appendix C. Two lemmata

Lemma C.1. Deÿne
k n+j
—¦ (k)
A= ; (C.1)
j
(n + j)r+1
j=0

—¦ —¦ (k)
k
(k) j
where is a zero of multiplicity m of (z) = j=0 j z . Then
—¦
(’1)m d m (k)
r+m
n+m
A∼ () (n ’ ∞): (C.2)
r nr+m+1 d xm

Proof. Use

1 1
exp(’at)t r dt;
= a¿0 (C.3)
ar+1 r! 0

to obtain
k
∞ n ∞
1 —¦
—¦ (k) n+j r (k)
( exp(’t))t r dt:
A= exp(’(n + j)t)t dt = exp(’nt) (C.4)
j
r! r!
0 0
j=0

Taylor expansion of the polynomial yields due to the zero at
—¦
(’ )m d m (k) (x)
—¦
(k)
t m (1 + O(t)):
( exp(’t)) = (C.5)
m! d xm
x=

Invoking Watson™s lemma [6, p. 263 ] completes the proof.

—¦ (k)
(k)
Lemma C.2. Assume that assumption (C-3 ) of Theorem 13 holds. Further assume ’ for
n; j j
n ’ ∞. Then; Eq. (280) holds.

Proof. We have
j’1
!n+j n+t
j j
∼ exp ∼ exp(j n ) (C.6)
n
!n n
t=0

for large n. Hence,
k k
!n —¦
—¦ (k)
(k)
exp( n ))’j = (k)
∼ j( (1= + n) (C.7)
n; j
!n+j
j=0 j=0

—¦
(k)
Since the characteristic polynomial (z) has a zero of order at z = 1= according to the
assumptions, Eq. (280) follows using Taylor expansion.
144 H.H.H. Homeier / Journal of Computational and Applied Mathematics 122 (2000) 81“147


References

[1] M. Abramowitz, I. Stegun, Handbook of Mathematical Functions, Dover Publications, New York, 1970.
[2] A.C. Aitken, On Bernoulli™s numerical solution of algebraic equations, Proc. Roy. Soc. Edinburgh 46 (1926)
289“305.
[3] G.A. Baker Jr., Essentials of Padà Approximants, Academic Press, New York, 1975.
e
[4] G.A. Baker Jr., P. Graves-Morris, Padà Approximants, Part I: Basic Theory, Addison-Wesley, Reading, MA, 1981.
e
[5] G.A. Baker Jr., P. Graves-Morris, Padà Approximants, 2nd Edition, Cambridge Unievrsity Press, Cambridge, GB,
e
1996.
[6] C.M. Bender, S.A. Orszag, Advanced Mathematical Methods for Scientists and Engineers, McGraw-Hill, Singapore,
1987.
[7] C. Brezinski, Accà là ration de la Convergence en Analyse Numà rique, Springer, Berlin, 1977.
ee e
à Ã
[8] C. Brezinski, Algorithmes d™Accà là ration de la Convergence “ Etude Numà rique, Editions Technip, Paris, 1978.
ee e
[9] C. Brezinski, A general extrapolation algorithm, Numer. Math. 35 (1980) 175“180.
[10] C. Brezinski, Padà -Type Approximation and General Orthogonal Polynomials, Birkhauser, Basel, 1980.
e
[11] C. Brezinski, A Bibliography on Continued Fractions, Padà Approximation, Extrapolation and Related Subjects,
e
Prensas Universitarias de Zaragoza, Zaragoza, 1991.
[12] C. Brezinski (Ed.), Continued Fractions and Padà Approximants, North-Holland, Amsterdam, 1991.
e
[13] C. Brezinski, A.C. Matos, A derivation of extrapolation algorithms based on error estimates, J. Comput. Appl.
Math. 66 (1“2) (1996) 5“26.
[14] C. Brezinski, M. Redivo Zaglia, Extrapolation Methods, Theory and Practice, North-Holland, Amsterdam, 1991.
[15] C. Brezinski, M. Redivo Zaglia, A general extrapolation algorithm revisited, Adv. Comput. Math. 2 (1994) 461“477.
[16] J. Cioslowski, E.J. Weniger, Bulk properties from ÿnite cluster calculations, VIII Benchmark calculations on the
e ciency of extrapolation methods for the HF and MP2 energies of polyacenes, J. Comput. Chem. 14 (1993)
1468“1481.
„„
[17] J. CÃzek, F. Vinette, E.J. Weniger, Examples on the use of symbolic computation in physics and chemistry:
applications of the inner projection technique and of a new summation method for divergent series, Int. J. Quantum
Chem. Symp. 25 (1991) 209“223.
„„
[18] J. CÃzek, F. Vinette, E.J. Weniger, On the use of the symbolic language Maple in physics and chemistry:
several examples, in: de Groot R.A, Nadrchal J. (Eds.), Proceedings of the Fourth International Conference on
Computational Physics PHYSICS COMPUTING ™92, World Scientiÿc, Singapore, 1993, pp. 31“44.
[19] J.E. Drummond, A formula for accelerating the convergence of a general series, Bull. Austral. Math. Soc. 6 (1972)
69“74.
[20] A. Erdà lyi, W. Magnus, F. Oberhettinger, F.G. Tricomi, Higher Transcendental Functions, Vol. I, McGraw-Hill,
e
New York, 1953.
[21] T. Fessler, W.F. Ford, D.A. Smith, HURRY: an acceleration algorithm for scalar sequences and series, ACM Trans.
Math. Software 9 (1983) 346“354.
[22] W.F. Ford, A. Sidi, An algorithm for a generalization of the Richardson extrapolation process, SIAM J. Numer.
Anal. 24 (5) (1987) 1212“1232.
[23] B. Germain-Bonne, Transformations de suites, Rev. Francaise Automat. Inform. Rech. Oper. 7 (R-1) (1973) 84“90.
[24] P.R. Graves-Morris (Ed.), Padà Approximants, The Institute of Physics, London, 1972.
e
[25] P.R. Graves-Morris (Ed.), Padà Approximants and their Applications, Academic Press, London, 1973.
e
[26] P.R. Graves-Morris, E.B. Sa , R.S. Varga (Ed.), Rational Approximation and Interpolation, Springer, Berlin, 1984.
[27] J. Grotendorst, A Maple package for transforming series, sequences and functions, Comput. Phys. Comm. 67 (1991)
325“342.
[28] J. Grotendorst, E.O. Steinborn, Use of nonlinear convergence accelerators for the e cient evaluation of GTO
molecular integrals, J. Chem. Phys. 84 (1986) 5617“5623.
[29] J. Grotendorst, E.J. Weniger, E.O. Steinborn, E cient evaluation of inÿnite-series representations for overlap,
two-center nuclear attraction and Coulomb integrals using nonlinear convergence accelerators, Phys. Rev. A 33
(1986) 3706“3726.
[30] E.R. Hansen, A Table of Series and Products, Prentice-Hall, Englewood-Cli s, NJ, 1975.
H.H.H. Homeier / Journal of Computational and Applied Mathematics 122 (2000) 81“147 145


[31] T. Hasegawa, A. Sidi, An automatic integration procedure for inÿnite range integrals involving oscillatory kernels,
Numer. Algorithms 13 (1996) 1“19.
[32] T. Havie, Generalized Neville type extrapolation schemes, BIT 19 (1979) 204“213.
[33] H.H.H. Homeier, Integraltransformationsmethoden und Quadraturverfahren fur Molekulintegrale mit B-Funktionen,
Theorie und Forschung, Vol. 121, S. Roderer Verlag, Regensburg, 1990, also Doctoral dissertation, Universitat
Regensburg.
[34] H.H.H. Homeier, A Levin-type algorithm for accelerating the convergence of Fourier series, Numer. Algorithms 3
(1992) 245“254.
[35] H.H.H. Homeier, Some applications of nonlinear convergence accelerators, Int. J. Quantum Chem. 45 (1993) 545“
562.
[36] H.H.H. Homeier, A hierarchically consistent, iterative sequence transformation, Numer. Algorithms 8 (1994) 47“81.
[37] H.H.H. Homeier, Nonlinear convergence acceleration for orthogonal series, in: R. Gruber, M. Tomassini (Eds.),
Proceedings of the sixth Joint EPS“APS International Conference on Physics Computing, Physics Computing ™94,
European Physical Society, Boite Postale 69, CH-1213 Petit-Lancy, Genf, Schweiz, 1994, pp. 47“50.
[38] H.H.H. Homeier, Determinantal representations for the J transformation, Numer. Math. 71 (3) (1995) 275“288.
[39] H.H.H. Homeier, Analytical and numerical studies of the convergence behavior of the J transformation, J. Comput.
Appl. Math. 69 (1996) 81“112.
[40] H.H.H. Homeier, Extrapolationsverfahren fur Zahlen-, Vektor- und Matrizenfolgen und ihre Anwendung in der
Theoretischen und Physikalischen Chemie, Habilitation Thesis, Universitat Regensburg, 1996.
[41] H.H.H. Homeier, Extended complex series methods for the convergence acceleration of Fourier series, Technical
Report TC-NA-97-3, Institut fur Physikalische und Theoretische Chemie, Universitat Regensburg, D-93040
Regensburg, 1997.
[42] H.H.H. Homeier, On an extension of the complex series method for the convergence acceleration of orthogonal
expansions. Technical Report TC-NA-97-4, Institut fur Physikalische und Theoretische Chemie, Universitat
Regensburg, D-93040 Regensburg, 1997.

<<

. 34
( 83 .)



>>