Author
Leslie Greengard
Other affiliations: Yale University, National Institute of Standards and Technology, Courant Institute of Mathematical Sciences ...read more
Bio: Leslie Greengard is an academic researcher from New York University. The author has contributed to research in topics: Integral equation & Fast multipole method. The author has an hindex of 53, co-authored 205 publications receiving 17857 citations. Previous affiliations of Leslie Greengard include Yale University & National Institute of Standards and Technology.
Papers published on a yearly basis
Papers
More filters
••
TL;DR: An algorithm is presented for the rapid evaluation of the potential and force fields in systems involving large numbers of particles whose interactions are Coulombic or gravitational in nature, making it considerably more practical for large-scale problems encountered in plasma physics, fluid dynamics, molecular dynamics, and celestial mechanics.
5,003 citations
•
07 Apr 1988
TL;DR: In this paper, an algorithm for the rapid evaluation of the potential and force fields in large-scale ensembles of particles is presented, which requires an amount of work proportional to the number of particles.
Abstract: The evaluation of Coulombic or gravitational interactions in large-scale ensembles of particles is an integral part of the numerical simulation of a large number of physical processes. Examples include celestial mechanics, plasma physics, the vortex method in fluid dynamics, molecular dynamics, and the solution of the Laplace equation via potential theory. In a typical application, a numerical model follows the trajectories of a number of particles moving in accordance with Newton's second law of motion in a field generated by the whole ensemble. In many situations, in order to be of physical interest, the simulation has to involve thousands of particles (or more), and the fields have to be evaluated for a large number of configurations. Unfortunately, an amount of work of the order $O(N\sp 2)$ has traditionally been required to evaluate all pairwise interactions in a system of N particles, unless some approximation or truncation method is used. As a result, large-scale simulations have been extremely expensive in some cases, and prohibitive in others. We present an algorithm for the rapid evaluation of the potential and force fields in large-scale systems of particles. In order to evaluate all pairwise Coulombic interactions of N particles to within round-off error, the algorithm requires an amount of work proportional to N, and this estimate does not depend on the statistics of the distribution. Both two and three dimensional versions of the algorithm have been constructed, and we will discuss their applications to several problems in physics, chemistry, biology, and numerical complex analysis.
1,145 citations
••
TL;DR: A new version of the Fast Multipole Method for the evaluation of potential fields in three dimensions is introduced based on a new diagonal form for translation operators and yields high accuracy at a reasonable cost.
Abstract: We introduce a new version of the Fast Multipole Method for the evaluation of potential fields in three dimensions. It is based on a new diagonal form for translation operators and yields high accuracy at a reasonable cost.
896 citations
••
TL;DR: This paper observes that one of the standard interpolation or "gridding" schemes, based on Gaussians, can be accelerated by a significant factor without precomputation and storage of the interpolation weights, of particular value in two- and three- dimensional settings.
Abstract: The nonequispaced Fourier transform arises in a variety of application areas, from medical imaging to radio astronomy to the numerical solution of partial differential equations. In a typical problem, one is given an irregular sampling of N data in the frequency domain and one is interested in reconstructing the corresponding function in the physical domain. When the sampling is uniform, the fast Fourier transform (FFT) allows this calculation to be computed in O(N log N ) operations rather than O(N 2 ) operations. Unfortunately, when the sampling is nonuniform, the FFT does not apply. Over the last few years, a number of algorithms have been developed to overcome this limitation and are often referred to as nonuniform FFTs (NUFFTs). These rely on a mixture of interpolation and the judicious use of the FFT on an oversampled grid (A. Dutt and V. Rokhlin, SIAM J. Sci. Comput., 14 (1993), pp. 1368-1383). In this paper, we observe that one of the standard interpolation or "gridding" schemes, based on Gaussians, can be accelerated by a significant factor without precomputation and storage of the interpolation weights. This is of particular value in two- and three- dimensional settings, saving either 10 d N in storage in d dimensions or a factor of about 5-10 in CPUtime (independent of dimension).
714 citations
••
TL;DR: An adaptive fast multipole method for the Laplace equation in three dimensions that uses both new compression techniques and diagonal forms for translation operators to achieve high accuracy at a reasonable cost.
671 citations
Cited by
More filters
••
[...]
TL;DR: There is, I think, something ethereal about i —the square root of minus one, which seems an odd beast at that time—an intruder hovering on the edge of reality.
Abstract: There is, I think, something ethereal about i —the square root of minus one. I remember first hearing about it at school. It seemed an odd beast at that time—an intruder hovering on the edge of reality.
Usually familiarity dulls this sense of the bizarre, but in the case of i it was the reverse: over the years the sense of its surreal nature intensified. It seemed that it was impossible to write mathematics that described the real world in …
33,785 citations
••
TL;DR: In this article, three parallel algorithms for classical molecular dynamics are presented, which can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors.
32,670 citations
01 May 1993
TL;DR: Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems.
Abstract: Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.
29,323 citations
••
TL;DR: An N⋅log(N) method for evaluating electrostatic energies and forces of large periodic systems is presented based on interpolation of the reciprocal space Ewald sums and evaluation of the resulting convolutions using fast Fourier transforms.
Abstract: An N⋅log(N) method for evaluating electrostatic energies and forces of large periodic systems is presented. The method is based on interpolation of the reciprocal space Ewald sums and evaluation of the resulting convolutions using fast Fourier transforms. Timings and accuracies are presented for three large crystalline ionic systems.
24,332 citations
••
TL;DR: It is demonstrated that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N), which is comparable to that of a simple truncation method of 10 A or less.
Abstract: The previously developed particle mesh Ewald method is reformulated in terms of efficient B‐spline interpolation of the structure factors This reformulation allows a natural extension of the method to potentials of the form 1/rp with p≥1 Furthermore, efficient calculation of the virial tensor follows Use of B‐splines in place of Lagrange interpolation leads to analytic gradients as well as a significant improvement in the accuracy We demonstrate that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N) For biomolecular systems with many thousands of atoms this method permits the use of Ewald summation at a computational cost comparable to that of a simple truncation method of 10 A or less
17,897 citations