3,219 research outputs found
Zolotarev Quadrature Rules and Load Balancing for the FEAST Eigensolver
The FEAST method for solving large sparse eigenproblems is equivalent to
subspace iteration with an approximate spectral projector and implicit
orthogonalization. This relation allows to characterize the convergence of this
method in terms of the error of a certain rational approximant to an indicator
function. We propose improved rational approximants leading to FEAST variants
with faster convergence, in particular, when using rational approximants based
on the work of Zolotarev. Numerical experiments demonstrate the possible
computational savings especially for pencils whose eigenvalues are not well
separated and when the dimension of the search space is only slightly larger
than the number of wanted eigenvalues. The new approach improves both
convergence robustness and load balancing when FEAST runs on multiple search
intervals in parallel.Comment: 22 pages, 8 figure
The geometric mean of two matrices from a computational viewpoint
The geometric mean of two matrices is considered and analyzed from a
computational viewpoint. Some useful theoretical properties are derived and an
analysis of the conditioning is performed. Several numerical algorithms based
on different properties and representation of the geometric mean are discussed
and analyzed and it is shown that most of them can be classified in terms of
the rational approximations of the inverse square root functions. A review of
the relevant applications is given
GPU-accelerated discontinuous Galerkin methods on hybrid meshes
We present a time-explicit discontinuous Galerkin (DG) solver for the
time-domain acoustic wave equation on hybrid meshes containing vertex-mapped
hexahedral, wedge, pyramidal and tetrahedral elements. Discretely energy-stable
formulations are presented for both Gauss-Legendre and Gauss-Legendre-Lobatto
(Spectral Element) nodal bases for the hexahedron. Stable timestep restrictions
for hybrid meshes are derived by bounding the spectral radius of the DG
operator using order-dependent constants in trace and Markov inequalities.
Computational efficiency is achieved under a combination of element-specific
kernels (including new quadrature-free operators for the pyramid), multi-rate
timestepping, and acceleration using Graphics Processing Units.Comment: Submitted to CMAM
Rational approximation to the fractional Laplacian operator in reaction-diffusion problems
This paper provides a new numerical strategy to solve fractional in space
reaction-diffusion equations on bounded domains under homogeneous Dirichlet
boundary conditions. Using the matrix transform method the fractional Laplacian
operator is replaced by a matrix which, in general, is dense. The approach here
presented is based on the approximation of this matrix by the product of two
suitable banded matrices. This leads to a semi-linear initial value problem in
which the matrices involved are sparse. Numerical results are presented to
verify the effectiveness of the proposed solution strategy
- …