91,157 research outputs found
A least-squares implicit RBF-FD closest point method and applications to PDEs on moving surfaces
The closest point method (Ruuth and Merriman, J. Comput. Phys.
227(3):1943-1961, [2008]) is an embedding method developed to solve a variety
of partial differential equations (PDEs) on smooth surfaces, using a closest
point representation of the surface and standard Cartesian grid methods in the
embedding space. Recently, a closest point method with explicit time-stepping
was proposed that uses finite differences derived from radial basis functions
(RBF-FD). Here, we propose a least-squares implicit formulation of the closest
point method to impose the constant-along-normal extension of the solution on
the surface into the embedding space. Our proposed method is particularly
flexible with respect to the choice of the computational grid in the embedding
space. In particular, we may compute over a computational tube that contains
problematic nodes. This fact enables us to combine the proposed method with the
grid based particle method (Leung and Zhao, J. Comput. Phys. 228(8):2993-3024,
[2009]) to obtain a numerical method for approximating PDEs on moving surfaces.
We present a number of examples to illustrate the numerical convergence
properties of our proposed method. Experiments for advection-diffusion
equations and Cahn-Hilliard equations that are strongly coupled to the velocity
of the surface are also presented
Integration over curves and surfaces defined by the closest point mapping
We propose a new formulation for integrating over smooth curves and surfaces
that are described by their closest point mappings. Our method is designed for
curves and surfaces that are not defined by any explicit parameterization and
is intended to be used in combination with level set techniques. However,
contrary to the common practice with level set methods, the volume integrals
derived from our formulation coincide exactly with the surface or line
integrals that one wish to compute. We study various aspects of this
formulation and provide a geometric interpretation of this formulation in terms
of the singular values of the Jacobian matrix of the closest point mapping.
Additionally, we extend the formulation - initially derived to integrate over
manifolds of codimension one - to include integration along curves in three
dimensions. Some numerical examples using very simple discretizations are
presented to demonstrate the efficacy of the formulation.Comment: Revised the pape
Computing largest circles separating two sets of segments
A circle separates two planar sets if it encloses one of the sets and its
open interior disk does not meet the other set. A separating circle is a
largest one if it cannot be locally increased while still separating the two
given sets. An Theta(n log n) optimal algorithm is proposed to find all largest
circles separating two given sets of line segments when line segments are
allowed to meet only at their endpoints. In the general case, when line
segments may intersect times, our algorithm can be adapted to
work in O(n alpha(n) log n) time and O(n \alpha(n)) space, where alpha(n)
represents the extremely slowly growing inverse of the Ackermann function.Comment: 14 pages, 3 figures, abstract presented at 8th Canadian Conference on
Computational Geometry, 199
Algorithms for Stable Matching and Clustering in a Grid
We study a discrete version of a geometric stable marriage problem originally
proposed in a continuous setting by Hoffman, Holroyd, and Peres, in which
points in the plane are stably matched to cluster centers, as prioritized by
their distances, so that each cluster center is apportioned a set of points of
equal area. We show that, for a discretization of the problem to an
grid of pixels with centers, the problem can be solved in time , and we experiment with two slower but more practical algorithms and
a hybrid method that switches from one of these algorithms to the other to gain
greater efficiency than either algorithm alone. We also show how to combine
geometric stable matchings with a -means clustering algorithm, so as to
provide a geometric political-districting algorithm that views distance in
economic terms, and we experiment with weighted versions of stable -means in
order to improve the connectivity of the resulting clusters.Comment: 23 pages, 12 figures. To appear (without the appendices) at the 18th
International Workshop on Combinatorial Image Analysis, June 19-21, 2017,
Plovdiv, Bulgari
Solving eigenvalue problems on curved surfaces using the Closest Point Method
Eigenvalue problems are fundamental to mathematics and science. We present a
simple algorithm for determining eigenvalues and eigenfunctions of the
Laplace--Beltrami operator on rather general curved surfaces. Our algorithm,
which is based on the Closest Point Method, relies on an embedding of the
surface in a higher-dimensional space, where standard Cartesian finite
difference and interpolation schemes can be easily applied. We show that there
is a one-to-one correspondence between a problem defined in the embedding space
and the original surface problem. For open surfaces, we present a simple way to
impose Dirichlet and Neumann boundary conditions while maintaining second-order
accuracy. Convergence studies and a series of examples demonstrate the
effectiveness and generality of our approach
Computing a Compact Spline Representation of the Medial Axis Transform of a 2D Shape
We present a full pipeline for computing the medial axis transform of an
arbitrary 2D shape. The instability of the medial axis transform is overcome by
a pruning algorithm guided by a user-defined Hausdorff distance threshold. The
stable medial axis transform is then approximated by spline curves in 3D to
produce a smooth and compact representation. These spline curves are computed
by minimizing the approximation error between the input shape and the shape
represented by the medial axis transform. Our results on various 2D shapes
suggest that our method is practical and effective, and yields faithful and
compact representations of medial axis transforms of 2D shapes.Comment: GMP14 (Geometric Modeling and Processing
- …