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

Leslie Greengard

Leslie Greengard contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAmath-phmath.MPmath.APNumerical Analysiseess.IVmath.DGphysics.comp-phphysics.plasm-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
9topics
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

10 published item(s)

preprint2023arXiv

FMM-LU: A fast direct solver for multiscale boundary integral equations in three dimensions

We present a fast direct solver for boundary integral equations on complex surfaces in three dimensions using an extension of the recently introduced recursive strong skeletonization scheme. For problems that are not highly oscillatory, our algorithm computes an ${LU}$-like hierarchical factorization of the dense system matrix, permitting application of the inverse in $\mathcal O(n)$ time, where $n$ is the number of unknowns on the surface. The factorization itself also scales linearly with the system size, albeit with a somewhat larger constant. The scheme is built on a level-restricted adaptive octree data structure, and therefore it is compatible with highly nonuniform discretizations. Furthermore, the scheme is coupled with high-order accurate locally-corrected Nyström quadrature methods to integrate the singular and weakly-singular Green's functions used in the integral representations. Our method has immediate applications to a variety of problems in computational physics. We concentrate here on studying its performance in acoustic scattering (governed by the Helmholtz equation) at low to moderate frequencies, and provide rigorous justification for compression of submatrices via proxy surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Geometry of the Phase Retrieval Problem

One of the most powerful approaches to imaging at the nanometer or subnanometer length scale is coherent diffraction imaging using X-ray sources. For amorphous (non-crystalline) samples, the raw data can be interpreted as the modulus of the continuous Fourier transform of the unknown object. Making use of prior information about the sample (such as its support), a natural goal is to recover the phase through computational means, after which the unknown object can be visualized at high resolution. While many algorithms have been proposed for this phase retrieval problem, careful analysis of its well-posedness has received relatively little attention. In this paper, we show that the problem is, in general, not well-posed and describe some of the underlying issues that are responsible for the ill-posedness. We then show how this analysis can be used to develop experimental protocols that lead to better conditioned inverse problems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

High-order discretization of a stable time-domain integral equation for 3D acoustic scattering

We develop a high-order, explicit method for acoustic scattering in three space dimensions based on a combined-field time-domain integral equation. The spatial discretization, of Nyström type, uses Gaussian quadrature on panels combined with a special treatment of the weakly singular kernels arising in near-neighbor interactions. In time, a new class of convolution splines is used in a predictor-corrector algorithm. Experiments on a torus and a perturbed torus are used to explore the stability and accuracy of the proposed scheme. This involved around one thousand solver runs, at up to 8th order and up to around 20,000 spatial unknowns, demonstrating 5-9 digits of accuracy. In addition we show that parameters in the combined field formulation, chosen on the basis of analysis for the sphere and other convex scatterers, work well in these cases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Debye Sources, Beltrami Fields, and a Complex Structure on Maxwell Fields

The Debye source representation for solutions to the time harmonic Maxwell equations is extended to bounded domains with finitely many smooth boundary components. A strong uniqueness result is proved for this representation. Natural complex structures are identified on the vector spaces of time-harmonic Maxwell fields. It is shown that in terms of Debye source data, these complex structures are uniformized, that is, represented by a fixed linear map on a fixed vector space, independent of the frequency. This complex structure relates time-harmonic Maxwell fields to constant-k Beltrami fields, i.e. solutions of the equation curl(E) = kE. A family of self-adjoint boundary conditions are defined for the Beltrami operator. This leads to a proof of the existence of zero-flux, constant-k, force-free Beltrami fields for any bounded region in R^3, as well as a constructive method to find them. The family of self-adjoint boundary value problems defines a new spectral invariant for bounded domains in R^3.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Efficient sum-of-exponentials approximations for the heat kernel and their applications

In this paper, we show that efficient separated sum-of-exponentials approximations can be constructed for the heat kernel in any dimension. In one space dimension, the heat kernel admits an approximation involving a number of terms that is of the order $O(\log(\frac{T}δ) (\log(\frac{1}ε)+\log\log(\frac{T}δ)))$ for any $x\in\bbR$ and $δ\leq t \leq T$, where $ε$ is the desired precision. In all higher dimensions, the corresponding heat kernel admits an approximation involving only $O(\log^2(\frac{T}δ))$ terms for fixed accuracy $ε$. These approximations can be used to accelerate integral equation-based methods for boundary value problems governed by the heat equation in complex geometry. The resulting algorithms are nearly optimal. For $N_S$ points in the spatial discretization and $N_T$ time steps, the cost is $O(N_S N_T \log^2 \frac{T}δ)$ in terms of both memory and CPU time for fixed accuracy $ε$. The algorithms can be parallelized in a straightforward manner. Several numerical examples are presented to illustrate the accuracy and stability of these approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Norm-preserving discretization of integral equations for elliptic PDEs with internal layers I: the one-dimensional case

We investigate the behavior of integral formulations of variable coefficient elliptic partial differential equations (PDEs) in the presence of steep internal layers. In one dimension, the equations that arise can be solved analytically and the condition numbers estimated in various L_p norms. We show that high-order accurate Nyström discretization leads to well-conditioned finite dimensional linear systems if and only if the discretization is both norm-preserving in a correctly chosen L_p space and adaptively refined in the internal layer.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the convergence of local expansions of layer potentials

In a recently developed quadrature method (quadrature by expansion or QBX), it was demonstrated that weakly singular or singular layer potentials can be evaluated rapidly and accurately on surface by making use of local expansions about carefully chosen off-surface points. In this paper, we derive estimates for the rate of convergence of these local expansions, providing the analytic foundation for the QBX method. The estimates may also be of mathematical interest, particularly for microlocal or asymptotic analysis in potential theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the efficient representation of the half-space impedance Green's function for the Helmholtz equation

A classical problem in acoustic (and electromagnetic) scattering concerns the evaluation of the Green's function for the Helmholtz equation subject to impedance boundary conditions on a half-space. The two principal approaches used for representing this Green's function are the Sommerfeld integral and the (closely related) method of complex images. The former is extremely efficient when the source is at some distance from the half-space boundary, but involves an unwieldy range of integration as the source gets closer and closer. Complex image-based methods, on the other hand, can be quite efficient when the source is close to the boundary, but they do not easily permit the use of the superposition principle since the selection of complex image locations depends on both the source and the target. We have developed a new, hybrid representation which uses a finite number of real images (dependent only on the source location) coupled with a rapidly converging Sommerfeld-like integral. While our method applies in both two and three dimensions, we restrict the detailed analysis and numerical experiments here to the two-dimensional case.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Debye Sources and the Numerical Solution of the Time Harmonic Maxwell Equations, II

In this paper, we develop a new integral representation for the solution of the time harmonic Maxwell equations in media with piecewise constant dielectric permittivity and magnetic permeability in R^3. This representation leads to a coupled system of Fredholm integral equations of the second kind for four scalar densities supported on the material interface. Like the classical Muller equation, it has no spurious resonances. Unlike the classical approach, however, the representation does not suffer from low frequency breakdown. We illustrate the performance of the method with numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fast multi-particle scattering: a hybrid solver for the Maxwell equations in microstructured materials

A variety of problems in device and materials design require the rapid forward modeling of Maxwell's equations in complex micro-structured materials. By combining high-order accurate integral equation methods with classical multiple scattering theory, we have created an effective simulation tool for materials consisting of an isotropic background in which are dispersed a large number of micro- or nano-scale metallic or dielectric inclusions.

0 citations0 reviewsNo code linkedOpen access