English

A New Algorithm for the Higher-Order $G$-Transformation

Numerical Analysis 2017-06-07 v1

Abstract

Let the scalars An(j)A^{(j)}_n be defined via the linear equations Al=An(j)+k=1nαˉkuk+l1,  l=j,j+1,,j+n .A_l=A^{(j)}_n+\sum^n_{k=1}\bar{\alpha}_ku_{k+l-1},\ \ l=j,j+1,\ldots,j+n\ . Here the AiA_i and uiu_i are known and the αˉk\bar{\alpha}_k are additional unknowns, and the quantities of interest are the An(j)A^{(j)}_n. This problem arises, for example, when one computes infinite-range integrals by the higher-order GG-transformation of Gray, Atchison, and McWilliams. One efficient procedure for computing the An(j)A^{(j)}_n is the rs-algorithm of Pye and Atchison. In the present work, we develop yet another procedure that combines the FS-algorithm of Ford and Sidi and the qd-algorithm of Rutishauser, and we denote it the FS/qd-algorithm. We show that the FS/qd-algorithm has a smaller operation count than the rs-algorithm. We also show that the FS/qd algorithm can also be used to implement the transformation of Shanks, and compares very favorably with the ε\varepsilon-algorithm of Wynn that is normally used for this purpose.

Keywords

Cite

@article{arxiv.1706.01786,
  title  = {A New Algorithm for the Higher-Order $G$-Transformation},
  author = {Avram Sidi},
  journal= {arXiv preprint arXiv:1706.01786},
  year   = {2017}
}
R2 v1 2026-06-22T20:10:37.080Z