About the Project

Wynn%20epsilon%20algorithm

AdvancedHelp

(0.003 seconds)

1—10 of 222 matching pages

1: 22.16 Related Functions
§22.16(ii) Jacobi’s Epsilon Function
Integral Representations
See Figure 22.16.2. …
Quasi-Addition and Quasi-Periodic Formulas
Relation to Theta Functions
2: 18.40 Methods of Computation
For applications in which the OP’s appear only as terms in series expansions (compare §18.18(i)) the need to compute them can be avoided altogether by use instead of Clenshaw’s algorithm3.11(ii)) and its straightforward generalization to OP’s other than Chebyshev. … A simple set of choices is spelled out in Gordon (1968) which gives a numerically stable algorithm for direct computation of the recursion coefficients in terms of the moments, followed by construction of the J-matrix and quadrature weights and abscissas, and we will follow this approach: Let N be a positive integer and define … It is now necessary to take the limit ε 0 + of F ( x + i ε ) , and the imaginary part is the required Stieltjes–Perron inversion: …Results of low ( 2 to 3 decimal digits) precision for w ( x ) are easily obtained for N 10 to 20 . Gautschi (2004, p. 119–120) has explored the ε 0 + limit via the Wynn ε -algorithm, (3.9.11) to accelerate convergence, finding four to eight digits of precision in w ( x ) , depending smoothly on x , for N 4000 , for an example involving first numerator Legendre OP’s. …
3: 33.24 Tables
§33.24 Tables
  • Abramowitz and Stegun (1964, Chapter 14) tabulates F 0 ( η , ρ ) , G 0 ( η , ρ ) , F 0 ( η , ρ ) , and G 0 ( η , ρ ) for η = 0.5 ( .5 ) 20 and ρ = 1 ( 1 ) 20 , 5S; C 0 ( η ) for η = 0 ( .05 ) 3 , 6S.

  • Curtis (1964a) tabulates P ( ϵ , r ) , Q ( ϵ , r ) 33.1), and related functions for = 0 , 1 , 2 and ϵ = 2 ( .2 ) 2 , with x = 0 ( .1 ) 4 for ϵ < 0 and x = 0 ( .1 ) 10 for ϵ 0 ; 6D.

  • 4: 20 Theta Functions
    Chapter 20 Theta Functions
    5: 3.9 Acceleration of Convergence
    The ratio of the Hankel determinants in (3.9.9) can be computed recursively by Wynn’s epsilon algorithm:
    ε 1 ( n ) = 0 ,
    Then t n , 2 k = ε 2 k ( n ) . … If s n is the n th partial sum of a power series f , then t n , 2 k = ε 2 k ( n ) is the Padé approximant [ ( n + k ) / k ] f 3.11(iv)). For further information on the epsilon algorithm see Brezinski and Redivo Zaglia (1991, pp. 78–95). …
    6: 3.11 Approximation Techniques
    The Padé approximants can be computed by Wynn’s cross rule. Any five approximants arranged in the Padé table as … For the recursive computation of [ n + k / k ] f by Wynn’s epsilon algorithm, see (3.9.11) and the subsequent text. … is of fundamental importance in the FFT algorithm. …For further details and algorithms, see Van Loan (1992). …
    7: Bibliography W
  • R. S. Ward (1987) The Nahm equations, finite-gap potentials and Lamé functions. J. Phys. A 20 (10), pp. 2679–2683.
  • T. Weider (1999) Algorithm 794: Numerical Hankel transform by the Fortran program HANKEL. ACM Trans. Math. Software 25 (2), pp. 240–250.
  • E. J. Weniger (2003) A rational approximant for the digamma function. Numer. Algorithms 33 (1-4), pp. 499–507.
  • M. E. Wojcicki (1961) Algorithm 44: Bessel functions computed recursively. Comm. ACM 4 (4), pp. 177–178.
  • P. Wynn (1966) Upon systems of recursions which obtain among the quotients of the Padé table. Numer. Math. 8 (3), pp. 264–269.
  • 8: 2.11 Remainder Terms; Stokes Phenomenon
    Similar improvements are achievable by Aitken’s Δ 2 -process, Wynn’s ϵ -algorithm, and other acceleration transformations. … For example, using double precision d 20 is found to agree with (2.11.31) to 13D. …
    9: 8 Incomplete Gamma and Related
    Functions
    10: 28 Mathieu Functions and Hill’s Equation