BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

J. D. Ward

J. D. Ward contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.CAmath.NA

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

preprint2012arXiv

Kernel Based Quadrature on Spheres and Other Homogeneous Spaces

Quadrature formulas for spheres, the rotation group, and other compact, homogeneous manifolds are important in a number of applications and have been the subject of recent research. The main purpose of this paper is to study coordinate independent quadrature (or cubature) formulas associated with certain classes of positive definite and conditionally positive definite kernels that are invariant under the group action of the homogeneous manifold. In particular, we show that these formulas are accurate -- optimally so in many cases -- and stable under an increasing number of nodes and in the presence of noise, provided the set X of quadrature nodes is quasi-uniform. The stability results are new in all cases. In addition, we may use these quadrature formulas to obtain similar formulas for manifolds diffeomorphic to $\mathbb{S}^n$, oblate spheroids for instance. The weights are obtained by solving a single linear system. For $\mathbb{S}^2$, and the restricted thin plate spline kernel $r^2 \log r$, these weights can be computed for two-thirds of a million nodes, using a preconditioned iterative technique introduced by us.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Localized bases for kernel spaces on the unit sphere

Approximation/interpolation from spaces of positive definite or conditionally positive definite kernels is an increasingly popular tool for the analysis and synthesis of scattered data, and is central to many meshless methods. For a set of $N$ scattered sites, the standard basis for such a space utilizes $N$ \emph{globally} supported kernels; computing with it is prohibitively expensive for large $N$. Easily computable, well-localized bases, with "small-footprint" basis elements - i.e., elements using only a small number of kernels -- have been unavailable. Working on $\sphere$, with focus on the restricted surface spline kernels (e.g. the thin-plate splines restricted to the sphere), we construct easily computable, spatially well-localized, small-footprint, robust bases for the associated kernel spaces. Our theory predicts that each element of the local basis is constructed by using a combination of only $\mathcal{O}((\log N)^2)$ kernels, which makes the construction computationally cheap. We prove that the new basis is $L_p$ stable and satisfies polynomial decay estimates that are stationary with respect to the density of the data sites, and we present a quasi-interpolation scheme that provides optimal $L_p$ approximation orders. Although our focus is on $\mathbb{S}^2$, much of the theory applies to other manifolds - $\mathbb{S}^d$, the rotation group, and so on. Finally, we construct algorithms to implement these schemes and use them to conduct numerical experiments, which validate our theory for interpolation problems on $\mathbb{S}^2$ involving over one hundred fifty thousand data sites.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Better bases for kernel spaces

In this article we investigate the feasibility of constructing stable, local bases for computing with kernels. In particular, we are interested in constructing families $(b_ξ)_{ξ\inΞ}$ that function as bases for kernel spaces $S(k,Ξ)$ so that each basis function is constructed using very few kernels. In other words, each function $b_ζ(x) = \sum_{ξ\inΞ} A_{ζ,ξ} k(x,ξ)$ is a linear combination of samples of the kernel with few nonzero coefficients $A_{ζ,ξ}$. This is reminiscent of the construction of the B-spline basis from the family of truncated power functions. We demonstrate that for a large class of kernels (the Sobolev kernels as well as many kernels of polyharmonic and related type) such bases exist. In fact, the basis elements can be constructed using a combination of roughly $O(\log N)^d$ kernels, where $d$ is the local dimension of the manifold and $N$ is the dimension of the kernel space (i.e. $N=#Ξ$). Viewing this as a preprocessing step -- the construction of the basis has computational cost $O(N(\log N)^d)$. Furthermore, we prove that the new basis is $L_p$ stable and satisfies polynomial decay estimates that are stationary with respect to the density of $Ξ$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Polyharmonic and Related Kernels on Manifolds: Interpolation and Approximation

This article is devoted to developing a theory for effective kernel interpolation and approximation in a general setting. For a wide class of compact, connected $C^\infty$ Riemannian manifolds, including the important cases of spheres and SO(3), we establish, using techniques involving differential geometry and Lie groups, that the kernels obtained as fundamental solutions of certain partial differential operators generate Lagrange functions that are uniformly bounded and decay away from their center at an algebraic rate, and in certain cases, an exponential rate. An immediate corollary is that the corresponding Lebesgue constants for interpolation as well as for $L_2$ minimization are uniformly bounded with a constant whose only dependence on the set of data sites is reflected in the mesh ratio, which measures the uniformity of the data. The kernels considered here include the restricted surface splines on spheres, as well as surface splines for SO(3), both of which have elementary closed-form representations that are computationally implementable. In addition to obtaining bounded Lebesgue constants in this setting, we also establish a "zeros lemma" for domains on compact Riemannian manifolds -- one that holds in as much generality as the corresponding Euclidean zeros lemma (on Lipschitz domains satisfying interior cone conditions) with constants that clearly demonstrate the influence of the geometry of the boundary (via cone parameters) as well as that of the Riemannian metric.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cardinal Interpolation with Gaussian Kernels

In this paper, interpolation by scaled multi-integer translates of Gaussian kernels is studied. The main result establishes $L_p$ Sobolev error estimates and shows that the error is controlled by the $L_p$ multiplier norm of a Fourier multiplier closely related to the cardinal interpolant, and comparable to the Hilbert transform. Consequently, its multiplier norm is bounded independent of the grid spacing when $1<p<\infty$, and involves a logarithmic term when $p=1$ or $\infty$.

0 citations0 reviewsNo code linkedOpen access