BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ian H. Sloan

Ian H. Sloan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath.CAmath.PRmath.AP

Catalog footprint

What is connected

14works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

14 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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fast approximation by periodic kernel-based lattice-point interpolation with application in uncertainty quantification

This paper deals with the kernel-based approximation of a multivariate periodic function by interpolation at the points of an integration lattice -- a setting that, as pointed out by Zeng, Leung, Hickernell (MCQMC2004, 2006) and Zeng, Kritzer, Hickernell (Constr. Approx., 2009), allows fast evaluation by fast Fourier transform, so avoiding the need for a linear solver. The main contribution of the paper is the application to the approximation problem for uncertainty quantification of elliptic partial differential equations, with the diffusion coefficient given by a random field that is periodic in the stochastic variables, in the model proposed recently by Kaarnioja, Kuo, Sloan (SIAM J. Numer. Anal., 2020). The paper gives a full error analysis, and full details of the construction of lattices needed to ensure a good (but inevitably not optimal) rate of convergence and an error bound independent of dimension. Numerical experiments support the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quasi-Monte Carlo finite element analysis for wave propagation in heterogeneous random media

We propose and analyze a quasi-Monte Carlo (QMC) algorithm for efficient simulation of wave propagation modeled by the Helmholtz equation in a bounded region in which the refractive index is random and spatially heterogenous. Our focus is on the case in which the region can contain multiple wavelengths. We bypass the usual sign-indefiniteness of the Helmholtz problem by switching to an alternative sign-definite formulation recently developed by Ganesh and Morgenstern (Numerical Algorithms, 83, 1441-1487, 2020). The price to pay is that the regularity analysis required for QMC methods becomes much more technical. Nevertheless we obtain a complete analysis with error comprising stochastic dimension truncation error, finite element error and cubature error, with results comparable to those obtained for the diffusion problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Uncertainty quantification using periodic random variables

Many studies in uncertainty quantification have been carried out under the assumption of an input random field in which a countable number of independent random variables are each uniformly distributed on an interval, with these random variables entering linearly in the input random field (the so-called affine model). In this paper we consider an alternative model of the random field, in which the random variables have the same uniform distribution on an interval, but the random variables enter the input field as periodic functions. The field is constructed in such a way as to have the same mean and covariance function as the affine random field. Higher moments differ from the affine case, but in general the periodic model seems no less desirable. The new model of the random field is used to compute expected values of a quantity of interest arising from an elliptic PDE with random coefficients. The periodicity is shown to yield a higher order cubature convergence rate of $\mathcal{O}(n^{-1/p})$ independently of the dimension when used in conjunction with rank-1 lattice cubature rules constructed using suitably chosen smoothness-driven product and order dependent weights, where $n$ is the number of lattice points and $p$ is the summability exponent of the fluctuations in the series expansion of the random coefficient. We present numerical examples that assess the performance of our method.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Infinite-dimensional integration and the multivariate decomposition method

We further develop the \emph{Multivariate Decomposition Method} (MDM) for the Lebesgue integration of functions of infinitely many variables $x_1,x_2,x_3,\ldots$ with respect to a corresponding product of a one dimensional probability measure. Although a number of concepts of infinite-dimensional integrals have been used in the literature, questions of uniqueness and compatibility have mostly not been studied. We show that, under appropriate convergence conditions, the Lebesgue integral equals the `anchored' integral, independently of the anchor. The MDM assumes that point values of $f_{\mathfrak{u}}$ are available for important subsets ${\mathfrak{u}}$, at some known cost. In this paper we introduce a new setting, in which it is assumed that each $f_{\mathfrak{u}}$ belongs to a normed space $F_{\mathfrak{u}}$, and that bounds $B_{\mathfrak{u}}$ on $\|f_{\mathfrak{u}}\|_{F_{\mathfrak{u}}}$ are known. This contrasts with the assumption in many papers that weights $γ_{\mathfrak{u}}$, appearing in the norm of the infinite-dimensional function space, are somehow known. Often such weights $γ_{\mathfrak{u}}$ were determined by minimizing an error bound depending on the $B_{\mathfrak{u}}$, the $γ_{\mathfrak{u}}$ \emph{and} the chosen algorithm, resulting in weights that depend on the algorithm. In contrast, in this paper only the bounds $B_{\mathfrak{u}}$ are assumed known. We give two examples in which we specialize the MDM: in the first case $F_{\mathfrak{u}}$ is the $|{\mathfrak{u}}|$-fold tensor product of an anchored reproducing kernel Hilbert space, and in the second case it is a particular non-Hilbert space for integration over an unbounded domain.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multilevel Quasi-Monte Carlo Methods for Lognormal Diffusion Problems

In this paper we present a rigorous cost and error analysis of a multilevel estimator based on randomly shifted Quasi-Monte Carlo (QMC) lattice rules for lognormal diffusion problems. These problems are motivated by uncertainty quantification problems in subsurface flow. We extend the convergence analysis in [Graham et al., Numer. Math. 2014] to multilevel Quasi-Monte Carlo finite element discretizations and give a constructive proof of the dimension-independent convergence of the QMC rules. More precisely, we provide suitable parameters for the construction of such rules that yield the required variance reduction for the multilevel scheme to achieve an $\varepsilon$-error with a cost of $\mathcal{O}(\varepsilon^{-θ})$ with $θ< 2$, and in practice even $θ\approx 1$, for sufficiently fast decaying covariance kernels of the underlying Gaussian random field inputs. This confirms that the computational gains due to the application of multilevel sampling methods and the gains due to the application of QMC methods, both demonstrated in earlier works for the same model problem, are complementary. A series of numerical experiments confirms these gains. The results show that in practice the multilevel QMC method consistently outperforms both the multilevel MC method and the single-level variants even for non-smooth problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Needlet approximation for isotropic random fields on the sphere

In this paper we establish a multiscale approximation for random fields on the sphere using spherical needlets --- a class of spherical wavelets. We prove that the semidiscrete needlet decomposition converges in mean and pointwise senses for weakly isotropic random fields on $\mathbb{S}^{d}$, $d\ge2$. For numerical implementation, we construct a fully discrete needlet approximation of a smooth $2$-weakly isotropic random field on $\mathbb{S}^{d}$ and prove that the approximation error for fully discrete needlets has the same convergence order as that for semidiscrete needlets. Numerical examples are carried out for fully discrete needlet approximations of Gaussian random fields and compared to a discrete version of the truncated Fourier expansion.

0 citations0 reviewsNo code linkedOpen access
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}$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On filtered polynomial approximation on the sphere

This paper considers filtered polynomial approximations on the unit sphere $\mathbb{S}^d\subset \mathbb{R}^{d+1}$, obtained by truncating smoothly the Fourier series of an integrable function $f$ with the help of a "filter" $h$, which is a real-valued continuous function on $[0,\infty)$ such that $h(t)=1$ for $t\in[0,1]$ and $h(t)=0$ for $t\ge2$. The resulting "filtered polynomial approximation" (a spherical polynomial of degree $2L-1$) is then made fully discrete by approximating the inner product integrals by an $N$-point cubature rule of suitably high polynomial degree of precision, giving an approximation called "filtered hyperinterpolation". In this paper we require that the filter $h$ and all its derivatives up to $\lfloor\tfrac{d-1}{2}\rfloor$ are absolutely continuous, while its right and left derivatives of order $\lfloor \tfrac{d+1}{2}\rfloor$ exist everywhere and are of bounded variation. Under this assumption we show that for a function $f$ in the Sobolev space $W^s_p(\mathbb{S}^d),\ 1\le p\le \infty$, both approximations are of the optimal order $ L^{-s}$, in the first case for $s>0$ and in the second fully discrete case for $s>d/p$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two-parameter regularization of ill-posed spherical pseudo-differential equations in the space of continuous functions

In this paper, a two-step regularization method is used to solve an ill-posed spherical pseudo-differential equation in the presence of noisy data. For the first step of regularization we approximate the data by means of a spherical polynomial that minimizes a functional with a penalty term consisting of the squared norm in a Sobolev space. The second step is a regularized collocation method. An error bound is obtained in the uniform norm, which is potentially smaller than that for either the noise reduction alone or the regularized collocation alone. We discuss an a posteriori parameter choice, and present some numerical experiments, which support the claimed superiority of the two-step method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multi-level quasi-Monte Carlo finite element methods for a class of elliptic partial differential equations with random coefficients

Quasi-Monte Carlo (QMC) methods are applied to multi-level Finite Element (FE) discretizations of elliptic partial differential equations (PDEs) with a random coefficient, to estimate expected values of linear functionals of the solution. The expected value is considered as an infinite-dimensional integral in the parameter space corresponding to the randomness induced by the random coefficient. We use a multi-level algorithm, with the number of QMC points depending on the discretization level, and with a level-dependent dimension truncation strategy. In some scenarios, we show that the overall error is $\mathcal{O}(h^2)$, where $h$ is the finest FE mesh width, or $\mathcal{O}(N^{-1+δ})$ for arbitrary $δ>0$, where $N$ is the maximal number of QMC sampling points. For these scenarios, the total work is essentially of the order of one single PDE solve at the finest FE discretization level. The analysis exploits regularity of the parametric solution with respect to both the physical variables (the variables in the physical domain) and the parametric variables (the parameters corresponding to randomness). Families of QMC rules with "POD weights" ("product and order dependent weights") which quantify the relative importance of subsets of the variables are found to be natural for proving convergence rates of QMC errors that are independent of the number of parametric variables.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stability and preconditioning for a hybrid approximation on the sphere

This paper proposes a new preconditioning scheme for a linear system with a saddle-point structure arising from a hybrid approximation scheme on the sphere, an approximation scheme that combines (local) spherical radial basis functions and (global) spherical polynomials. Making use of a recently derived inf-sup condition [13] and the Brezzi stability and convergence theorem for this approximation scheme, we show that the linear system can be optimally preconditioned with a suitable block-diagonal preconditioner. Numerical experiments with a non-uniform distribution of data points support the theoretical conclusions.

0 citations0 reviewsNo code linkedOpen access