Möbius-Transformed Trapezoidal Rule

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Suzuki, Yuya, Hyvönen, Nuutti, Karvonen, Toni
Formato: Preprint
Publicado: 2024
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866917326145716224
author Suzuki, Yuya
Hyvönen, Nuutti
Karvonen, Toni
author_facet Suzuki, Yuya
Hyvönen, Nuutti
Karvonen, Toni
contents We study numerical integration by combining the trapezoidal rule with a Möbius transformation that maps the unit circle onto the real line. We prove that the resulting transformed trapezoidal rule attains the optimal rate of convergence if the integrand function lives in a weighted Sobolev space with a weight that is only assumed to be a positive Schwartz function decaying monotonically to zero close to infinity. Our algorithm only requires the ability to evaluate the weight at the selected nodes, and it does not require sampling from a probability measure defined by the weight nor information on its derivatives. In particular, we show that the Möbius transformation, as a change of variables between the real line and the unit circle, sends a function in the weighted Sobolev space to a periodic Sobolev space with the same smoothness. Since there are various results available for integrating and approximating periodic functions, we also describe several extensions of the Möbius-transformed trapezoidal rule, including function approximation via trigonometric interpolation, integration with randomized algorithms, and multivariate integration.
format Preprint
id arxiv_https___arxiv_org_abs_2407_13650
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Möbius-Transformed Trapezoidal Rule
Suzuki, Yuya
Hyvönen, Nuutti
Karvonen, Toni
Numerical Analysis
65D15, 65D30, 65D32, 68W20
We study numerical integration by combining the trapezoidal rule with a Möbius transformation that maps the unit circle onto the real line. We prove that the resulting transformed trapezoidal rule attains the optimal rate of convergence if the integrand function lives in a weighted Sobolev space with a weight that is only assumed to be a positive Schwartz function decaying monotonically to zero close to infinity. Our algorithm only requires the ability to evaluate the weight at the selected nodes, and it does not require sampling from a probability measure defined by the weight nor information on its derivatives. In particular, we show that the Möbius transformation, as a change of variables between the real line and the unit circle, sends a function in the weighted Sobolev space to a periodic Sobolev space with the same smoothness. Since there are various results available for integrating and approximating periodic functions, we also describe several extensions of the Möbius-transformed trapezoidal rule, including function approximation via trigonometric interpolation, integration with randomized algorithms, and multivariate integration.
title Möbius-Transformed Trapezoidal Rule
topic Numerical Analysis
65D15, 65D30, 65D32, 68W20
url https://arxiv.org/abs/2407.13650