1,701 research outputs found
Constructing Intrinsic Delaunay Triangulations of Submanifolds
We describe an algorithm to construct an intrinsic Delaunay triangulation of
a smooth closed submanifold of Euclidean space. Using results established in a
companion paper on the stability of Delaunay triangulations on -generic
point sets, we establish sampling criteria which ensure that the intrinsic
Delaunay complex coincides with the restricted Delaunay complex and also with
the recently introduced tangential Delaunay complex. The algorithm generates a
point set that meets the required criteria while the tangential complex is
being constructed. In this way the computation of geodesic distances is
avoided, the runtime is only linearly dependent on the ambient dimension, and
the Delaunay complexes are guaranteed to be triangulations of the manifold
The Harris-Luck criterion for random lattices
The Harris-Luck criterion judges the relevance of (potentially) spatially
correlated, quenched disorder induced by, e.g., random bonds, randomly diluted
sites or a quasi-periodicity of the lattice, for altering the critical behavior
of a coupled matter system. We investigate the applicability of this type of
criterion to the case of spin variables coupled to random lattices. Their
aptitude to alter critical behavior depends on the degree of spatial
correlations present, which is quantified by a wandering exponent. We consider
the cases of Poissonian random graphs resulting from the Voronoi-Delaunay
construction and of planar, ``fat'' Feynman diagrams and precisely
determine their wandering exponents. The resulting predictions are compared to
various exact and numerical results for the Potts model coupled to these
quenched ensembles of random graphs.Comment: 13 pages, 9 figures, 2 tables, REVTeX 4. Version as published, one
figure added for clarification, minor re-wordings and typo cleanu
Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations
One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. In the past, limited computer resources have motivated the development of efficient numerical schemes in computational fluid dynamics (CFD) utilizing structured meshes. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Unstructured grids on the other hand offer an alternative to modeling complex geometries. Unstructured meshes have irregular connectivity and usually contain combinations of triangles, quadrilaterals, tetrahedra, and hexahedra. The generation and use of unstructured grids poses new challenges in CFD. The purpose of this note is to present recent developments in the unstructured grid generation and flow solution technology
A probabilistic approach to reducing the algebraic complexity of computing Delaunay triangulations
Computing Delaunay triangulations in involves evaluating the
so-called in\_sphere predicate that determines if a point lies inside, on
or outside the sphere circumscribing points . This
predicate reduces to evaluating the sign of a multivariate polynomial of degree
in the coordinates of the points . Despite
much progress on exact geometric computing, the fact that the degree of the
polynomial increases with makes the evaluation of the sign of such a
polynomial problematic except in very low dimensions. In this paper, we propose
a new approach that is based on the witness complex, a weak form of the
Delaunay complex introduced by Carlsson and de Silva. The witness complex
is defined from two sets and in some metric space
: a finite set of points on which the complex is built, and a set of
witnesses that serves as an approximation of . A fundamental result of de
Silva states that if .
In this paper, we give conditions on that ensure that the witness complex
and the Delaunay triangulation coincide when is a finite set, and we
introduce a new perturbation scheme to compute a perturbed set close to
such that . Our perturbation
algorithm is a geometric application of the Moser-Tardos constructive proof of
the Lov\'asz local lemma. The only numerical operations we use are (squared)
distance comparisons (i.e., predicates of degree 2). The time-complexity of the
algorithm is sublinear in . Interestingly, although the algorithm does not
compute any measure of simplex quality, a lower bound on the thickness of the
output simplices can be guaranteed.Comment: 24 page
Delaunay Stability via Perturbations
We present an algorithm that takes as input a finite point set in Euclidean space, and performs a perturbation that guarantees that the Delaunay triangulation of the resulting perturbed point set has quantifiable stability with respect to the metric and the point positions. There is also a guarantee on the quality of the simplices: they cannot be too flat. The algorithm provides an alternative tool to the weighting or refinement methods to remove poorly shaped simplices in Delaunay triangulations of arbitrary dimension, but in addition it provides a guarantee of stability for the resulting triangulation
Adaptive, Anisotropic and Hierarchical cones of Discrete Convex functions
We address the discretization of optimization problems posed on the cone of
convex functions, motivated in particular by the principal agent problem in
economics, which models the impact of monopoly on product quality. Consider a
two dimensional domain, sampled on a grid of N points. We show that the cone of
restrictions to the grid of convex functions is in general characterized by N^2
linear inequalities; a direct computational use of this description therefore
has a prohibitive complexity. We thus introduce a hierarchy of sub-cones of
discrete convex functions, associated to stencils which can be adaptively,
locally, and anisotropically refined. Numerical experiments optimize the
accuracy/complexity tradeoff through the use of a-posteriori stencil refinement
strategies.Comment: 35 pages, 11 figures. (Second version fixes a small bug in Lemma 3.2.
Modifications are anecdotic.
- …