117,672 research outputs found
Inner and Outer Rounding of Boolean Operations on Lattice Polygonal Regions
Robustness problems due to the substitution of the exact computation on real
numbers by the rounded floating point arithmetic are often an obstacle to
obtain practical implementation of geometric algorithms. If the adoption of the
--exact computation paradigm--[Yap et Dube] gives a satisfactory solution to
this kind of problems for purely combinatorial algorithms, this solution does
not allow to solve in practice the case of algorithms that cascade the
construction of new geometric objects. In this report, we consider the problem
of rounding the intersection of two polygonal regions onto the integer lattice
with inclusion properties. Namely, given two polygonal regions A and B having
their vertices on the integer lattice, the inner and outer rounding modes
construct two polygonal regions with integer vertices which respectively is
included and contains the true intersection. We also prove interesting results
on the Hausdorff distance, the size and the convexity of these polygonal
regions
An exact general remeshing scheme applied to physically conservative voxelization
We present an exact general remeshing scheme to compute analytic integrals of
polynomial functions over the intersections between convex polyhedral cells of
old and new meshes. In physics applications this allows one to ensure global
mass, momentum, and energy conservation while applying higher-order polynomial
interpolation. We elaborate on applications of our algorithm arising in the
analysis of cosmological N-body data, computer graphics, and continuum
mechanics problems.
We focus on the particular case of remeshing tetrahedral cells onto a
Cartesian grid such that the volume integral of the polynomial density function
given on the input mesh is guaranteed to equal the corresponding integral over
the output mesh. We refer to this as "physically conservative voxelization".
At the core of our method is an algorithm for intersecting two convex
polyhedra by successively clipping one against the faces of the other. This
algorithm is an implementation of the ideas presented abstractly by Sugihara
(1994), who suggests using the planar graph representations of convex polyhedra
to ensure topological consistency of the output. This makes our implementation
robust to geometric degeneracy in the input. We employ a simplicial
decomposition to calculate moment integrals up to quadratic order over the
resulting intersection domain.
We also address practical issues arising in a software implementation,
including numerical stability in geometric calculations, management of
cancellation errors, and extension to two dimensions. In a comparison to recent
work, we show substantial performance gains. We provide a C implementation
intended to be a fast, accurate, and robust tool for geometric calculations on
polyhedral mesh elements.Comment: Code implementation available at https://github.com/devonmpowell/r3
Engineering Art Galleries
The Art Gallery Problem is one of the most well-known problems in
Computational Geometry, with a rich history in the study of algorithms,
complexity, and variants. Recently there has been a surge in experimental work
on the problem. In this survey, we describe this work, show the chronology of
developments, and compare current algorithms, including two unpublished
versions, in an exhaustive experiment. Furthermore, we show what core
algorithmic ingredients have led to recent successes
Flood propagation modelling with the Local Inertia Approximation: theoretical and numerical analysis of its physical limitations
Attention of the researchers has increased towards a simplification of the
complete Shallow water Equations called the Local Inertia Approximation (LInA),
which is obtained by neglecting the advection term in the momentum conservation
equation. In the present paper it is demonstrated that a shock is always
developed at moving wetting-drying frontiers, and this justifies the study of
the Riemann problem on even and uneven beds. In particular, the general exact
solution for the Riemann problem on horizontal frictionless bed is given,
together with the exact solution of the non-breaking wave propagating on
horizontal bed with friction, while some example solution is given for the
Riemann problem on discontinuous bed. From this analysis, it follows that
drying of the wet bed is forbidden in the LInA model, and that there are
initial conditions for which the Riemann problem has no solution on smoothly
varying bed. In addition, propagation of the flood on discontinuous sloping bed
is impossible if the bed drops height have the same order of magnitude of the
moving-frontier shock height. Finally, it is found that the conservation of the
mechanical energy is violated. It is evident that all these findings pose a
severe limit to the application of the model. The numerical analysis has proven
that LInA numerical models may produce numerical solutions, which are
unreliable because of mere algorithmic nature, also in the case that the LInA
mathematical solutions do not exist. The applicability limits of the LInA model
are discouragingly severe, even if the bed elevation varies continuously. More
important, the non-existence of the LInA solution in the case of discontinuous
topography and the non-existence of receding fronts radically question the
viability of the LInA model in realistic cases. It is evident that classic SWE
models should be preferred in the majority of the practical applications
- …