Adaption allows efficient integration of functions with unknown singularities

Leszek Plaskota, Grzegorz W. Wasilkowski

Research output: Contribution to journalArticlepeer-review

23 Scopus citations

Abstract

We study numerical integration [InlineMediaObject not available: see fulltext.] for functions f with singularities. Nonadaptive methods are inefficient in this case, and we show that the problem can be efficiently solved by adaptive quadratures at cost similar to that for functions with no singularities. Consider first a class [InlineMediaObject not available: see fulltext.] of functions whose derivatives of order up to r are continuous and uniformly bounded for any but one singular point. We propose adaptive quadratures Q* n , each using at most n function values, whose worst case errors [InlineMediaObject not available: see fulltext.] are proportional to n - r . On the other hand, the worst case error of nonadaptive methods does not converge faster than n -1. These worst case results do not extend to the case of functions with two or more singularities; however, adaption shows its power even for such functions in the asymptotic setting. That is, let F r be the class of r-smooth functions with arbitrary (but finite) number of singularities. Then a generalization of Q* n yields adaptive quadratures Q** n such that |I(f)-Q** n (f)|=O(n - r ) for any f F r . In addition, we show that for any sequence of nonadaptive methods there are 'many' functions in F r for which the errors converge no faster than n -1. Results of numerical experiments are also presented.

Original languageEnglish
Pages (from-to)123-144
Number of pages22
JournalNumerische Mathematik
Volume102
Issue number1
DOIs
StatePublished - Nov 2005

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Adaption allows efficient integration of functions with unknown singularities'. Together they form a unique fingerprint.

Cite this