7,391 research outputs found
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
Trefftz Difference Schemes on Irregular Stencils
The recently developed Flexible Local Approximation MEthod (FLAME) produces
accurate difference schemes by replacing the usual Taylor expansion with
Trefftz functions -- local solutions of the underlying differential equation.
This paper advances and casts in a general form a significant modification of
FLAME proposed recently by Pinheiro & Webb: a least-squares fit instead of the
exact match of the approximate solution at the stencil nodes. As a consequence
of that, FLAME schemes can now be generated on irregular stencils with the
number of nodes substantially greater than the number of approximating
functions. The accuracy of the method is preserved but its robustness is
improved. For demonstration, the paper presents a number of numerical examples
in 2D and 3D: electrostatic (magnetostatic) particle interactions, scattering
of electromagnetic (acoustic) waves, and wave propagation in a photonic
crystal. The examples explore the role of the grid and stencil size, of the
number of approximating functions, and of the irregularity of the stencils.Comment: 28 pages, 12 figures; to be published in J Comp Phy
- …