Meshless moment-free quadrature formulas arising from numerical differentiation
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_ | 1866915458394882048 |
|---|---|
| author | Davydov, Oleg Esposti, Bruno Degli |
| author_facet | Davydov, Oleg Esposti, Bruno Degli |
| contents | We suggest a method for simultaneously generating high order quadrature weights for integrals over Lipschitz domains and their boundaries that requires neither meshing nor moment computation. The weights are determined on pre-defined scattered nodes as a minimum norm solution of a sparse underdetermined linear system arising from a discretization of a suitable boundary value problem by either collocation or meshless finite differences. The method is easy to implement independently of the domain's representation, since it only requires as inputs the position of all quadrature nodes and the direction of outward-pointing normals at each node belonging to the boundary. Numerical experiments demonstrate the robustness and high accuracy of the method on a number of smooth and piecewise smooth domains in 2D and 3D, including some with reentrant corners and edges. Comparison with quadrature schemes provided by the state-of-the-art open source packages Gmsh and MFEM shows that the new method is competitive in terms of accuracy for a given number of nodes. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2409_03567 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Meshless moment-free quadrature formulas arising from numerical differentiation Davydov, Oleg Esposti, Bruno Degli Numerical Analysis 65D32, 41A55 We suggest a method for simultaneously generating high order quadrature weights for integrals over Lipschitz domains and their boundaries that requires neither meshing nor moment computation. The weights are determined on pre-defined scattered nodes as a minimum norm solution of a sparse underdetermined linear system arising from a discretization of a suitable boundary value problem by either collocation or meshless finite differences. The method is easy to implement independently of the domain's representation, since it only requires as inputs the position of all quadrature nodes and the direction of outward-pointing normals at each node belonging to the boundary. Numerical experiments demonstrate the robustness and high accuracy of the method on a number of smooth and piecewise smooth domains in 2D and 3D, including some with reentrant corners and edges. Comparison with quadrature schemes provided by the state-of-the-art open source packages Gmsh and MFEM shows that the new method is competitive in terms of accuracy for a given number of nodes. |
| title | Meshless moment-free quadrature formulas arising from numerical differentiation |
| topic | Numerical Analysis 65D32, 41A55 |
| url | https://arxiv.org/abs/2409.03567 |