ANZIAM J. 46(E) pp.C690--C703, 2005.

Sloan iteration and Richardson extrapolation for Walsh series solutions of integral equations

W. F. Blyth

V. Uljanov
(Received 14 February 2005, revised 24 June 2005)

Abstract

Algorithms have been developed by the authors and co-workers for the solution of both Volterra and Fredholm integral equations by using the discontinuous wavelet packets known as the Walsh functions. These Walsh function methods are typically globally convergent of order one (and locally of order two).

The usual Walsh function method with m = 2n terms approximates the solution with a piecewise constant function (constant on each sub-interval of width 1/m). For Fredholm integral equations, Sloan iteration of the Walsh series solution is globally convergent of order two.
For linear Fredholm integral equations of the second kind, we show in this paper that the Sloan iterates for Walsh series solution with m and 2m terms can be extrapolated by the Richardson method to give a function which approximates the solution with global fourth order convergence. This is very easy and efficient to implement. A second Richardson iteration for smooth problems results in sixth order convergence.
With Kulkarni's example (where the solution is only twice differentiable), a second Richardson iteration results in fifth order convergence. Our new method outperforms Kulkarni's method: it is easier implement, more accurate with 64 or more terms and the order of convergence is one higher.
As a further example, we use the nonlinear Chandrasekar integral equation for which a new solution method was recently proposed by the authors. Here, our numerical experiment shows that the Richardson extrapolation of the Sloan iteration functions is third order convergent.
The Walsh function methods are intuitively simple and robust. The straightforward implementation of Walsh series leads to schemes with low order convergence. However the implementation of the acceleration of convergence techniques described here is easy and efficient and provides schemes of higher order convergence.

Download to your computer

Authors

W. F. Blyth
V. Uljanov
Dept. Mathematics & Statistics, RMIT University, Melbourne, Australia. mailto:bill.blyth@rmit.edu.au, vladimir.uljanov@rmit.edu.au

Published July 26, 2005, amended August 2, 2005. ISSN 1446-8735

References

  1. K. E. Atkinson. The Numerical Solution of Integral Equations of the Second Kind. Cambridge University Press, 1997.
  2. W. F. Blyth, R. L. May and P. Widyaningsih. The Solution of Integral Equations Using Walsh Functions and a Multigrid Approach. In J. Noye et al., editors, Computational Techniques and Applications: CTAC97, Proceedings, 99--106. World Scientific, Singapore, 1998.
  3. W. F. Blyth, R. L. May and P. Widyaningsih. Solution of Separable Volterra Integral Equations Using Walsh Functions and a Multigrid Approach. In M. Pemberton, I. Turner and P. Jacobs, editors, Engineering Mathematics and Applications Conference: EMAC 2002 Proceedings, 85--90. The Inst. of Engineers, Australia, 2002.
  4. W. F. Blyth, R. L. May and P. Widyaningsih. Volterra Integral Equations solved in Fredholm form using Walsh functions. ANZIAM J., 45(E), ppC269--C282, 2004. http://anziamj.austms.org.au/V45/CTAC2003/Blyt
  5. W. F. Blyth and V. Uljanov. Numerical Solution of Weakly Singular Fredholm Integral Equations using Walsh Functions. In R. L. May and A. K. Easton, editors, Computational Techniques and Applications: CTAC-95, Proceedings, 137--143. World Scientific, Singapore, 1996.
  6. W. F. Blyth and V. Uljanov. Numerical Solution of the Chandrasekar Nonlinear Fredholm Integral Equation using Walsh Functions. in preparation.
  7. N. J. Fine. On the Walsh Functions. Trans. Amer. Math. Soc., 65, 372--414, 1949.
  8. A. Haar. Zur Theorie der orthogonalen Funktionen systeme. Math. Ann., 69, 331--371, 1910.
  9. R. Kress. Linear Integral Equations. Springer--Verlag, 1989.
  10. R. Kulkarni. A superconvergence result for solutions of compact operator equations. Bull. Austral. Math. Soc., 68, 517--528, 2003.
  11. R. E. A. C. Paley. A remarkable series of orthogonal functions (I & II). London Math. Soc., 34, 241--264 & 265--279, 1932.
  12. I. H. Sloan. Improvement by Iteration for Compact Operator Equations. Math. Comp. 30, 758--764 , 1976.
  13. B. G. Sloss and W. F. Blyth. A Walsh function method for a nonlinear Volterra integral equation. J. Franklin Inst., 340, 25--41, 2003.
  14. P. Wojtaszczyk. A mathematical introduction to wavelets. Cambridge University Press, 1997.
Copyright Austral. Mathematical Soc.; prepared August 2, 2005, from Blyt.tex with anziamje.sty (v2.3f)