Source author record

Robert S. Womersley

Robert S. Womersley appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

5works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2022arXiv

Needlets Liberated

Spherical needlets were introduced by Narcowich, Petrushev, and Ward to provide a multiresolution sequence of polynomial approximations to functions on the sphere. The needlet construction makes use of integration rules that are exact for polynomials up to a given degree. The aim of the present paper is to relax the exactness of the integration rules by replacing them with QMC designs as introduced by Brauchart, Saff, Sloan, and Womersley (2014). Such integration rules (generalised here by allowing non-equal cubature weights) provide the same asymptotic order of convergence as exact rules for Sobolev spaces $\mathbb{H}^s$, but are easier to obtain numerically. With such rules we construct ``generalised needlets''. The paper provides an error analysis that allows the replacement of the original needlets by generalised needlets, and more generally, analyses a hybrid scheme in which the needlets for the lower levels are of the traditional kind, whereas the new generalised needlets are used for some number of higher levels. Numerical experiments complete the paper.

preprint2020arXiv

Numerical computation of triangular complex spherical designs with small mesh ratio

This paper provides triangular spherical designs for the complex unit sphere $Ω^d$ by exploiting the natural correspondence between the complex unit sphere in $d$ dimensions and the real unit sphere in $2d-1$. The existence of triangular and square complex spherical $t$-designs with the optimal order number of points is established. A variational characterization of triangular complex designs is provided, with particular emphasis on numerical computation of efficient triangular complex designs with good geometric properties as measured by their mesh ratio. We give numerical examples of triangular spherical $t$-designs on complex unit spheres of dimension $d=2$ to $6$.

preprint2016arXiv

Fully discrete needlet approximation on the sphere

Spherical needlets are highly localized radial polynomials on the sphere $\mathbb{S}^{d}\subset \mathbb{R}^{d+1}$, $d\ge 2$, with centers at the nodes of a suitable cubature rule. The original semidiscrete spherical needlet approximation of Narcowich, Petrushev and Ward is not computable, in that the needlet coefficients depend on inner product integrals. In this work we approximate these integrals by a second quadrature rule with an appropriate degree of precision, to construct a fully discrete needlet approximation. We prove that the resulting approximation is equivalent to filtered hyperinterpolation, that is to a filtered Fourier-Laplace series partial sum with inner products replaced by appropriate cubature sums. It follows that the $\mathbb{L}_{p}$-error of discrete needlet approximation of order $J$ for $1 \le p \le \infty$ and $s > d/p$ has for a function $f$ in the Sobolev space $\mathbb{W}_{p}^{s}(\mathbb{S}^{d})$ the optimal rate of convergence in the sense of optimal recovery, namely $\mathcal{O}\bigl(2^{-J s}\bigr)$. Moreover, this is achieved with a filter function that is of smoothness class $C^{\lfloor \frac{d+3}{2}\rfloor}$, in contrast to the usually assumed $C^{\infty}$. A numerical experiment for a class of functions in known Sobolev smoothness classes gives $\mathbb{L}_2$ errors for the fully discrete needlet approximation that are almost identical to those for the original semidiscrete needlet approximation. Another experiment uses needlets over the whole sphere for the lower levels together with high-level needlets with centers restricted to a local region. The resulting errors are reduced in the local region away from the boundary, indicating that local refinement in special regions is a promising strategy.

preprint2016arXiv

Random Point Sets on the Sphere --- Hole Radii, Covering, and Separation

Geometric properties of $N$ random points distributed independently and uniformly on the unit sphere $\mathbb{S}^{d}\subset\mathbb{R}^{d+1}$ with respect to surface area measure are obtained and several related conjectures are posed. In particular, we derive asymptotics (as $N \to \infty$) for the expected moments of the radii of spherical caps associated with the facets of the convex hull of $N$ random points on $\mathbb{S}^{d}$. We provide conjectures for the asymptotic distribution of the scaled radii of these spherical caps and the expected value of the largest of these radii (the covering radius). Numerical evidence is included to support these conjectures. Furthermore, utilizing the extreme law for pairwise angles of Cai et al., we derive precise asymptotics for the expected separation of random points on $\mathbb{S}^{d}$.

preprint2016arXiv

Riemann localisation on the sphere

This paper first shows that the Riemann localisation property holds for the Fourier-Laplace series partial sum for sufficiently smooth functions on the two-dimensional sphere, but does not hold for spheres of higher dimension. By Riemann localisation on the sphere $\mathbb{S}^{d}\subset\mathbb{R}^{d+1}$, $d\ge2$, we mean that for a suitable subset $X$ of $\mathbb{L}_{p}(\mathbb{S}^{d})$, $1\le p\le \infty$, the $\mathbb{L}_{p}$-norm of the Fourier local convolution of $f\in X$ converges to zero as the degree goes to infinity. The Fourier local convolution of $f$ at $\boldsymbol{x}\in\mathbb{S}^{d}$ is the Fourier convolution with a modified version of $f$ obtained by replacing values of $f$ by zero on a neighbourhood of $\boldsymbol{x}$. The failure of Riemann localisation for $d>2$ can be overcome by considering a filtered version: we prove that for a sphere of any dimension and sufficiently smooth filter the corresponding local convolution always has the Riemann localisation property. Key tools are asymptotic estimates of the Fourier and filtered kernels.