3,667 research outputs found
Hierarchical Schur complement preconditioner for the stochastic Galerkin finite element methods
Use of the stochastic Galerkin finite element methods leads to large systems
of linear equations obtained by the discretization of tensor product solution
spaces along their spatial and stochastic dimensions. These systems are
typically solved iteratively by a Krylov subspace method. We propose a
preconditioner which takes an advantage of the recursive hierarchy in the
structure of the global matrices. In particular, the matrices posses a
recursive hierarchical two-by-two structure, with one of the submatrices block
diagonal. Each one of the diagonal blocks in this submatrix is closely related
to the deterministic mean-value problem, and the action of its inverse is in
the implementation approximated by inner loops of Krylov iterations. Thus our
hierarchical Schur complement preconditioner combines, on each level in the
approximation of the hierarchical structure of the global matrix, the idea of
Schur complement with loops for a number of mutually independent inner Krylov
iterations, and several matrix-vector multiplications for the off-diagonal
blocks. Neither the global matrix, nor the matrix of the preconditioner need to
be formed explicitly. The ingredients include only the number of stiffness
matrices from the truncated Karhunen-Lo\`{e}ve expansion and a good
preconditioned for the mean-value deterministic problem. We provide a condition
number bound for a model elliptic problem and the performance of the method is
illustrated by numerical experiments.Comment: 15 pages, 2 figures, 9 tables, (updated numerical experiments
Sparse approximate inverse preconditioners on high performance GPU platforms
Simulation with models based on partial differential equations often requires the solution of (sequences of) large and sparse algebraic linear systems. In multidimensional domains, preconditioned Krylov iterative solvers are often appropriate for these duties. Therefore, the search for efficient preconditioners for Krylov subspace methods is a crucial theme. Recent developments, especially in computing hardware, have renewed the interest in approximate inverse preconditioners in factorized form, because their application during the solution process can be more efficient. We present here some experiences focused on the approximate inverse preconditioners proposed by Benzi and Tůma from 1996 and the sparsification and inversion proposed by van Duin in 1999. Computational costs, reorderings and implementation issues are considered both on conventional and innovative computing architectures like Graphics Programming Units (GPUs)
Hybrid preconditioning for iterative diagonalization of ill-conditioned generalized eigenvalue problems in electronic structure calculations
The iterative diagonalization of a sequence of large ill-conditioned
generalized eigenvalue problems is a computational bottleneck in quantum
mechanical methods employing a nonorthogonal basis for {\em ab initio}
electronic structure calculations. We propose a hybrid preconditioning scheme
to effectively combine global and locally accelerated preconditioners for rapid
iterative diagonalization of such eigenvalue problems. In partition-of-unity
finite-element (PUFE) pseudopotential density-functional calculations,
employing a nonorthogonal basis, we show that the hybrid preconditioned block
steepest descent method is a cost-effective eigensolver, outperforming current
state-of-the-art global preconditioning schemes, and comparably efficient for
the ill-conditioned generalized eigenvalue problems produced by PUFE as the
locally optimal block preconditioned conjugate-gradient method for the
well-conditioned standard eigenvalue problems produced by planewave methods
A domain decomposing parallel sparse linear system solver
The solution of large sparse linear systems is often the most time-consuming
part of many science and engineering applications. Computational fluid
dynamics, circuit simulation, power network analysis, and material science are
just a few examples of the application areas in which large sparse linear
systems need to be solved effectively. In this paper we introduce a new
parallel hybrid sparse linear system solver for distributed memory
architectures that contains both direct and iterative components. We show that
by using our solver one can alleviate the drawbacks of direct and iterative
solvers, achieving better scalability than with direct solvers and more
robustness than with classical preconditioned iterative solvers. Comparisons to
well-known direct and iterative solvers on a parallel architecture are
provided.Comment: To appear in Journal of Computational and Applied Mathematic
Parallel accelerated cyclic reduction preconditioner for three-dimensional elliptic PDEs with variable coefficients
We present a robust and scalable preconditioner for the solution of
large-scale linear systems that arise from the discretization of elliptic PDEs
amenable to rank compression. The preconditioner is based on hierarchical
low-rank approximations and the cyclic reduction method. The setup and
application phases of the preconditioner achieve log-linear complexity in
memory footprint and number of operations, and numerical experiments exhibit
good weak and strong scalability at large processor counts in a distributed
memory environment. Numerical experiments with linear systems that feature
symmetry and nonsymmetry, definiteness and indefiniteness, constant and
variable coefficients demonstrate the preconditioner applicability and
robustness. Furthermore, it is possible to control the number of iterations via
the accuracy threshold of the hierarchical matrix approximations and their
arithmetic operations, and the tuning of the admissibility condition parameter.
Together, these parameters allow for optimization of the memory requirements
and performance of the preconditioner.Comment: 24 pages, Elsevier Journal of Computational and Applied Mathematics,
Dec 201
A multidomain spectral method for solving elliptic equations
We present a new solver for coupled nonlinear elliptic partial differential
equations (PDEs). The solver is based on pseudo-spectral collocation with
domain decomposition and can handle one- to three-dimensional problems. It has
three distinct features. First, the combined problem of solving the PDE,
satisfying the boundary conditions, and matching between different subdomains
is cast into one set of equations readily accessible to standard linear and
nonlinear solvers. Second, touching as well as overlapping subdomains are
supported; both rectangular blocks with Chebyshev basis functions as well as
spherical shells with an expansion in spherical harmonics are implemented.
Third, the code is very flexible: The domain decomposition as well as the
distribution of collocation points in each domain can be chosen at run time,
and the solver is easily adaptable to new PDEs. The code has been used to solve
the equations of the initial value problem of general relativity and should be
useful in many other problems. We compare the new method to finite difference
codes and find it superior in both runtime and accuracy, at least for the
smooth problems considered here.Comment: 31 pages, 8 figure
- …