Computing Weak Distance between the 2-Sphere and Its Nonsmooth Approximations

From MaRDI portal
Publication:6154969

DOI10.1137/23M1571216arXiv2304.06340OpenAlexW4391353441MaRDI QIDQ6154969FDOQ6154969


Authors: Kazuki Koga Edit this on Wikidata


Publication date: 16 February 2024

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

Abstract: A novel algorithm is proposed for quantitative comparisons between compact surfaces embedded in the three-dimensional Euclidian space. The key idea is to identify those objects with the associated surface measures and compute a weak distance between them using the Fourier transform on the ambient space. In particular, the inhomogeneous Sobolev norm of negative order for a difference between two surface measures is evaluated via the Plancherel theorem, which amounts to approximating an weighted integral norm of smooth data on the frequency space. This approach allows several advantages including high accuracy due to fast-converging numerical quadrature rules, acceleration by the nonuniform fast Fourier transform, and parallelization on many-core processors. In numerical experiments, the 2-sphere, which is an example whose Fourier transform is explicitly known, is compared with its icosahedral discretization, and it is observed that the piecewise linear approximations converge to the smooth object at the quadratic rate up to small truncation.


Full work available at URL: https://arxiv.org/abs/2304.06340







Cites Work






This page was built for publication: Computing Weak Distance between the 2-Sphere and Its Nonsmooth Approximations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6154969)