Results 1  10
of
127
MorseSmale Complexes for Piecewise Linear 3Manifolds
, 2003
"... We define the MorseSmale complex of a Morse function over a 3manifold as the overlay of the descending and ascending manifolds of all critical points. In the generic case, its 3dimensional cells are shaped like crystals and are separated by quadrangular faces. In this paper, we give a combinatori ..."
Abstract

Cited by 105 (28 self)
 Add to MetaCart
We define the MorseSmale complex of a Morse function over a 3manifold as the overlay of the descending and ascending manifolds of all critical points. In the generic case, its 3dimensional cells are shaped like crystals and are separated by quadrangular faces. In this paper, we give a combinatorial algorithm for constructing such complexes for piecewise linear data.
A Simple Mesh Generator in MATLAB
 SIAM Review
, 2004
"... Abstract. Creating a mesh is the first step in a wide range of applications, including scientific computing and computer graphics. An unstructured simplex mesh requires a choice of meshpoints (vertex nodes) and a triangulation. We want to offer a short and simple MATLAB code, described in more detai ..."
Abstract

Cited by 42 (2 self)
 Add to MetaCart
Abstract. Creating a mesh is the first step in a wide range of applications, including scientific computing and computer graphics. An unstructured simplex mesh requires a choice of meshpoints (vertex nodes) and a triangulation. We want to offer a short and simple MATLAB code, described in more detail than usual, so the reader can experiment (and add to the code) knowing the underlying principles. We find the node locations by solving for equilibrium in a truss structure (using piecewise linear forcedisplacement relations) and we reset the topology by the Delaunay algorithm. The geometry is described implicitly by its distance function. In addition to being much shorter and simpler than other meshing techniques, our algorithm typically produces meshes of very high quality. We discuss ways to improve the robustness and the performance, but our aim here is simplicity. Readers can download (and edit) the codes from
Hardwareassisted visibility sorting for unstructured volume rendering
 IEEE Transactions on Visualization and Computer Graphics
, 2005
"... Abstract—Harvesting the power of modern graphics hardware to solve the complex problem of realtime rendering of large unstructured meshes is a major research goal in the volume visualization community. While, for regular grids, texturebased techniques are wellsuited for current GPUs, the steps ne ..."
Abstract

Cited by 41 (12 self)
 Add to MetaCart
Abstract—Harvesting the power of modern graphics hardware to solve the complex problem of realtime rendering of large unstructured meshes is a major research goal in the volume visualization community. While, for regular grids, texturebased techniques are wellsuited for current GPUs, the steps necessary for rendering unstructured meshes are not so easily mapped to current hardware. We propose a novel volume rendering technique that simplifies the CPUbased processing and shifts much of the sorting burden to the GPU, where it can be performed more efficiently. Our hardwareassisted visibility sorting algorithm is a hybrid technique that operates in both objectspace and imagespace. In objectspace, the algorithm performs a partial sort of the 3D primitives in preparation for rasterization. The goal of the partial sort is to create a list of primitives that generate fragments in nearly sorted order. In imagespace, the fragment stream is incrementally sorted using a fixeddepth sorting network. In our algorithm, the objectspace work is performed by the CPU and the fragmentlevel sorting is done completely on the GPU. A prototype implementation of the algorithm demonstrates that the fragmentlevel sorting achieves rendering rates of between one and six million tetrahedral cells per second on an ATI Radeon 9800. Index Terms—Volume visualization, graphics processors, visibility sorting. 1
Quality Meshing with Weighted Delaunay Refinement
 SIAM J. Comput
, 2002
"... Delaunay meshes with bounded circumradius to shortest edge length ratio have been proposed in the past for quality meshing. The only poor quality tetrahedra called slivers that can occur in such a mesh can be eliminated by the sliver exudation method. This method has been shown to work for periodic ..."
Abstract

Cited by 39 (7 self)
 Add to MetaCart
Delaunay meshes with bounded circumradius to shortest edge length ratio have been proposed in the past for quality meshing. The only poor quality tetrahedra called slivers that can occur in such a mesh can be eliminated by the sliver exudation method. This method has been shown to work for periodic point sets, but not with boundaries. Recently a randomized pointplacement strategy has been proposed to remove slivers while conforming to a given boundary. In this paper we present a deterministic algorithm for generating a weighted Delaunay mesh which respects the input boundary and has no poor quality tetrahedron including slivers. This success is achieved by combining the weight pumping method for sliver exudation and the Delaunay refinement method for boundary conformation. We show that an incremental weight pumping can be mixed seamlessly with vertex insertions in our weighted Delaunay refinement paradigm. 1
Collision detection for deforming necklaces
 IN SYMP. ON COMPUTATIONAL GEOMETRY
, 2002
"... In this paper, we propose to study deformable necklaces — flexible chains of balls, called beads, in which only adjacent balls may intersect. Such objects can be used to model macromolecules, muscles, rope, and other ‘linear ’ objects in the physical world. In this paper, we exploit this linearity ..."
Abstract

Cited by 34 (11 self)
 Add to MetaCart
In this paper, we propose to study deformable necklaces — flexible chains of balls, called beads, in which only adjacent balls may intersect. Such objects can be used to model macromolecules, muscles, rope, and other ‘linear ’ objects in the physical world. In this paper, we exploit this linearity to develop geometric structures associated with necklaces that are useful in physical simulations. We show how these structures can be implemented efficiently and maintained under necklace deformation. In particular, we study a bounding volume hierarchy based on spheres built on a necklace. Such a hierarchy is easy to compute and is suitable for maintenance when the necklace deforms, as our theoretical and experimental results show. This hierarchy can be used for collision and selfcollision detection. In particular, we achieve an upper bound of O(nlog n) in two dimensions and O(n 2−2/d) in ddimensions, d ≥ 3, for collision checking. To our knowledge, this is the first subquadratic bound proved for a collision detection algorithm using predefined hierarchies. In addition, we show that the power diagram, with the help of some additional mechanisms, can be also used to detect selfcollisions of a necklace in certain ways complementary to the sphere hierarchy.
Dense Point Sets Have Sparse Delaunay Triangulations
"... Delaunay triangulations and Voronoi diagrams are one of the most thoroughly studies objects in computational geometry, with numerous applications including nearestneighbor searching, clustering, finiteelement mesh generation, deformable surface modeling, and surface reconstruction. Many algorithms ..."
Abstract

Cited by 30 (2 self)
 Add to MetaCart
Delaunay triangulations and Voronoi diagrams are one of the most thoroughly studies objects in computational geometry, with numerous applications including nearestneighbor searching, clustering, finiteelement mesh generation, deformable surface modeling, and surface reconstruction. Many algorithms in these application domains begin by constructing the Delaunay triangulation or Voronoi diagram of a set of points in R³. Since threedimensional Delaunay triangulations can have complexity Ω(n²) in the worst case, these algorithms have worstcase running time \Omega (n2). However, this behavior is almost never observed in practice except for highlycontrived inputs. For all practical purposes, threedimensional Delaunay triangulations appear to have linear complexity. This frustrating
Quality Meshing for Polyhedra with Small Angles
, 2004
"... We present an algorithm to compute a Delaunay mesh conforming to a polyhedron possibly with small input angles. The radiusedge ratio of most output tetrahedra are bounded by a constant, except possibly those that are provably close to small angles. Furthermore, the mesh is not unnecessarily dense i ..."
Abstract

Cited by 29 (8 self)
 Add to MetaCart
We present an algorithm to compute a Delaunay mesh conforming to a polyhedron possibly with small input angles. The radiusedge ratio of most output tetrahedra are bounded by a constant, except possibly those that are provably close to small angles. Furthermore, the mesh is not unnecessarily dense in the sense that the edge lengths are at least a constant fraction of the local feature sizes at the edge endpoints. This algorithm is simple to implement as it eliminates most of the computation of local feature sizes and explicit protective zones. Our experimental results validate that few skinny tetrahedra remain and they lie close to small acute input angles. 1
Discrete Laplace operators: No free lunch
, 2007
"... Discrete Laplace operators are ubiquitous in applications spanning geometric modeling to simulation. For robustness and efficiency, many applications require discrete operators that retain key structural properties inherent to the continuous setting. Building on the smooth setting, we present a set ..."
Abstract

Cited by 28 (0 self)
 Add to MetaCart
Discrete Laplace operators are ubiquitous in applications spanning geometric modeling to simulation. For robustness and efficiency, many applications require discrete operators that retain key structural properties inherent to the continuous setting. Building on the smooth setting, we present a set of natural properties for discrete Laplace operators for triangular surface meshes. We prove an important theoretical limitation: discrete Laplacians cannot satisfy all natural properties; retroactively, this explains the diversity of existing discrete Laplace operators. Finally, we present a family of operators that includes and extends wellknown and widelyused operators.