75,927 research outputs found
Nonoverlapping domain decomposition preconditioners for discontinuous Galerkin finite element methods in -type norms
We analyse the spectral bounds of nonoverlapping domain decomposition preconditioners for -version discontinuous Galerkin finite element methods in -type norms, which arise in applications to fully nonlinear Hamilton--Jacobi--Bellman partial differential equations. We show that for a symmetric model problem, the condition number of the preconditioned system is at most of order , where and are respectively the coarse and fine mesh sizes, and and are respectively the coarse and fine mesh polynomial degrees. This represents the first result for this class of methods that explicitly accounts for the dependence of the condition number on , and its sharpness is shown numerically. The key analytical tool is an original optimal order approximation result between fine and coarse discontinuous finite element spaces.\ud
\ud
We then go beyond the model problem and show computationally that these methods lead to efficient and competitive solvers in practical applications to nonsymmetric, fully nonlinear Hamilton--Jacobi--Bellman equations
Nonoverlapping domain decomposition preconditioners for discontinuous Galerkin approximations of Hamilton--Jacobi--Bellman equations
We analyse a class of nonoverlapping domain decomposition preconditioners for
nonsymmetric linear systems arising from discontinuous Galerkin finite element
approximation of fully nonlinear Hamilton--Jacobi--Bellman (HJB) partial
differential equations. These nonsymmetric linear systems are uniformly bounded
and coercive with respect to a related symmetric bilinear form, that is
associated to a matrix . In this work, we construct a
nonoverlapping domain decomposition preconditioner , that is based
on , and we then show that the effectiveness of the preconditioner
for solving the} nonsymmetric problems can be studied in terms of the condition
number . In particular, we establish the
bound , where
and are respectively the coarse and fine mesh sizes, and and
are respectively the coarse and fine mesh polynomial degrees. This represents
the first such result for this class of methods that explicitly accounts for
the dependence of the condition number on ; our analysis is founded upon an
original optimal order approximation result between fine and coarse
discontinuous finite element spaces. Numerical experiments demonstrate the
sharpness of this bound. Although the preconditioners are not robust with
respect to the polynomial degree, our bounds quantify the effect of the coarse
and fine space polynomial degrees. Furthermore, we show computationally that
these methods are effective in practical applications to nonsymmetric, fully
nonlinear HJB equations under -refinement for moderate polynomial degrees
Non-Rigid Puzzles
Shape correspondence is a fundamental problem in computer graphics and vision, with applications in various problems including animation, texture mapping, robotic vision, medical imaging, archaeology and many more. In settings where the shapes are allowed to undergo non-rigid deformations and only partial views are available, the problem becomes very challenging. To this end, we present a non-rigid multi-part shape matching algorithm. We assume to be given a reference shape and its multiple parts undergoing a non-rigid deformation. Each of these query parts can be additionally contaminated by clutter, may overlap with other parts, and there might be missing parts or redundant ones. Our method simultaneously solves for the segmentation of the reference model, and for a dense correspondence to (subsets of) the parts. Experimental results on synthetic as well as real scans demonstrate the effectiveness of our method in dealing with this challenging matching scenario
Linear response separation of a solid into atomic constituents: Li, Al, and their evolution under pressure
We present the first realization of the generalized pseudoatom concept
introduced by Ball, and adopt the name enatom to minimize confusion. This
enatom, which consists of a unique decomposition of the total charge density
(or potential) of any solid into a sum of overlapping atomiclike contributions
that move rigidly with the nuclei to first order, is calculated using
(numerical) linear response methods, and is analyzed for both fcc Li and Al at
pressures of 0, 35, and 50 GPa. These two simple fcc metals (Li is fcc and a
good superconductor in the 20-40 GPa range) show different physical behaviors
under pressure, which reflects the increasing covalency in Li and the lack of
it in Al. The nonrigid (deformation) parts of the enatom charge and potential
have opposite signs in Li and Al; they become larger under pressure only in Li.
These results establish a method of construction of the enatom, whose potential
can be used to obtain a real-space understanding of the vibrational properties
and electron-phonon interaction in solids.Comment: 13 pages, 9 figures, 1 table, V2: fixed problem with Fig. 7, V3:
minor correction
- …