Parsimonious convolution quadrature
Fuente:
arXiv
Salvato in:
| Autori principali: | , |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2024
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866914980501127168 |
|---|---|
| author | Melenk, Jens M. Nick, Jörg |
| author_facet | Melenk, Jens M. Nick, Jörg |
| contents | We present a method to rapidly approximate convolution quadrature (CQ) approximations, based on a piecewise polynomial interpolation of the Laplace domain operator, which we call the \emph{parsimonious} convolution quadrature method. For implicit Euler and second order backward difference formula based discretizations, we require $O(\sqrt{N}\log N)$ evaluations in the Laplace domain to approximate $N$ time steps of the convolution quadrature method to satisfactory accuracy. The methodology proposed here differentiates from the well-understood fast and oblivious convolution quadrature \cite{SLL06}, since it is applicable to Laplace domain operator families that are only defined and polynomially bounded on a positive half space, which includes acoustic and electromagnetic wave scattering problems. The methods is applicable to linear and nonlinear integral equations. To elucidate the core idea, we give a complete and extensive analysis of the simplest case and derive worst-case estimates for the performance of parsimonious CQ based on the implicit Euler method. For sectorial Laplace transforms, we obtain methods that require $O(\log^2 N)$ Laplace domain evaluations on the complex right-half space. We present different implementation strategies, which only differ slightly from the classical realization of CQ methods. Numerical experiments demonstrate the use of the method with a time-dependent acoustic scattering problem, which was discretized by the boundary element method in space. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2410_15079 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Parsimonious convolution quadrature Melenk, Jens M. Nick, Jörg Numerical Analysis 65M38, 65M80, 65R20 We present a method to rapidly approximate convolution quadrature (CQ) approximations, based on a piecewise polynomial interpolation of the Laplace domain operator, which we call the \emph{parsimonious} convolution quadrature method. For implicit Euler and second order backward difference formula based discretizations, we require $O(\sqrt{N}\log N)$ evaluations in the Laplace domain to approximate $N$ time steps of the convolution quadrature method to satisfactory accuracy. The methodology proposed here differentiates from the well-understood fast and oblivious convolution quadrature \cite{SLL06}, since it is applicable to Laplace domain operator families that are only defined and polynomially bounded on a positive half space, which includes acoustic and electromagnetic wave scattering problems. The methods is applicable to linear and nonlinear integral equations. To elucidate the core idea, we give a complete and extensive analysis of the simplest case and derive worst-case estimates for the performance of parsimonious CQ based on the implicit Euler method. For sectorial Laplace transforms, we obtain methods that require $O(\log^2 N)$ Laplace domain evaluations on the complex right-half space. We present different implementation strategies, which only differ slightly from the classical realization of CQ methods. Numerical experiments demonstrate the use of the method with a time-dependent acoustic scattering problem, which was discretized by the boundary element method in space. |
| title | Parsimonious convolution quadrature |
| topic | Numerical Analysis 65M38, 65M80, 65R20 |
| url | https://arxiv.org/abs/2410.15079 |