TY - JOUR
T1 - Adaption allows efficient integration of functions with unknown singularities
AU - Plaskota, Leszek
AU - Wasilkowski, Grzegorz W.
PY - 2005/11
Y1 - 2005/11
N2 - 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.
AB - 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.
UR - http://www.scopus.com/inward/record.url?scp=27844505175&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=27844505175&partnerID=8YFLogxK
U2 - 10.1007/s00211-005-0640-3
DO - 10.1007/s00211-005-0640-3
M3 - Article
AN - SCOPUS:27844505175
SN - 0029-599X
VL - 102
SP - 123
EP - 144
JO - Numerische Mathematik
JF - Numerische Mathematik
IS - 1
ER -