107 research outputs found
Some observations on weighted GMRES
We investigate the convergence of the weighted GMRES method for solving linear systems. Two different weighting variants are compared with unweighted GMRES for three model problems, giving a phenomenological explanation of cases where weighting improves convergence, and a case where weighting has no effect on the convergence. We also present new alternative implementations of the weighted Arnoldi algorithm which may be favorable in terms of computational complexity, and examine stability issues connected with these implementations. Two implementations of weighted GMRES are compared for a large number of examples. We find that weighted GMRES may outperform unweighted GMRES for some problems, but more often this method is not competitive with other Krylov subspace methods like GMRES with deflated restarting or BICGSTAB, in particular when a preconditioner is used
A new level-dependent coarsegrid correction scheme for indefinite Helmholtz problems
In this paper we construct and analyse a level-dependent coarsegrid
correction scheme for indefinite Helmholtz problems. This adapted multigrid
method is capable of solving the Helmholtz equation on the finest grid using a
series of multigrid cycles with a grid-dependent complex shift, leading to a
stable correction scheme on all levels. It is rigourously shown that the
adaptation of the complex shift throughout the multigrid cycle maintains the
functionality of the two-grid correction scheme, as no smooth modes are
amplified in or added to the error. In addition, a sufficiently smoothing
relaxation scheme should be applied to ensure damping of the oscillatory error
components. Numerical experiments on various benchmark problems show the method
to be competitive with or even outperform the current state-of-the-art
multigrid-preconditioned Krylov methods, like e.g. CSL-preconditioned GMRES or
BiCGStab.Comment: 21 page
Recycling Krylov Subspaces for Efficient Partitioned Solution of Aerostructural Adjoint Systems
Robust and efficient solvers for coupled-adjoint linear systems are crucial
to successful aerostructural optimization. Monolithic and partitioned
strategies can be applied. The monolithic approach is expected to offer better
robustness and efficiency for strong fluid-structure interactions. However, it
requires a high implementation cost and convergence may depend on appropriate
scaling and initialization strategies. On the other hand, the modularity of the
partitioned method enables a straightforward implementation while its
convergence may require relaxation. In addition, a partitioned solver leads to
a higher number of iterations to get the same level of convergence as the
monolithic one.
The objective of this paper is to accelerate the fluid-structure
coupled-adjoint partitioned solver by considering techniques borrowed from
approximate invariant subspace recycling strategies adapted to sequences of
linear systems with varying right-hand sides. Indeed, in a partitioned
framework, the structural source term attached to the fluid block of equations
affects the right-hand side with the nice property of quickly converging to a
constant value. We also consider deflation of approximate eigenvectors in
conjunction with advanced inner-outer Krylov solvers for the fluid block
equations. We demonstrate the benefit of these techniques by computing the
coupled derivatives of an aeroelastic configuration of the ONERA-M6 fixed wing
in transonic flow. For this exercise the fluid grid was coupled to a structural
model specifically designed to exhibit a high flexibility. All computations are
performed using RANS flow modeling and a fully linearized one-equation
Spalart-Allmaras turbulence model. Numerical simulations show up to 39%
reduction in matrix-vector products for GCRO-DR and up to 19% for the nested
FGCRO-DR solver.Comment: 42 pages, 21 figure
Absolute value preconditioning for symmetric indefinite linear systems
We introduce a novel strategy for constructing symmetric positive definite
(SPD) preconditioners for linear systems with symmetric indefinite matrices.
The strategy, called absolute value preconditioning, is motivated by the
observation that the preconditioned minimal residual method with the inverse of
the absolute value of the matrix as a preconditioner converges to the exact
solution of the system in at most two steps. Neither the exact absolute value
of the matrix nor its exact inverse are computationally feasible to construct
in general. However, we provide a practical example of an SPD preconditioner
that is based on the suggested approach. In this example we consider a model
problem with a shifted discrete negative Laplacian, and suggest a geometric
multigrid (MG) preconditioner, where the inverse of the matrix absolute value
appears only on the coarse grid, while operations on finer grids are based on
the Laplacian. Our numerical tests demonstrate practical effectiveness of the
new MG preconditioner, which leads to a robust iterative scheme with minimalist
memory requirements
On large-scale diagonalization techniques for the Anderson model of localization
We propose efficient preconditioning algorithms for an eigenvalue problem arising in quantum physics, namely the computation of a few interior eigenvalues and their associated eigenvectors for large-scale sparse real and symmetric indefinite matrices of the Anderson model
of localization. We compare the Lanczos algorithm in the 1987 implementation by Cullum and Willoughby with the shift-and-invert techniques in the implicitly restarted Lanczos method and in the JacobiāDavidson method. Our preconditioning approaches for the shift-and-invert symmetric indefinite linear system are based on maximum weighted matchings and algebraic multilevel incomplete
LDLT factorizations. These techniques can be seen as a complement to the alternative idea of using more complete pivoting techniques for the highly ill-conditioned symmetric indefinite Anderson matrices. We demonstrate the effectiveness and the numerical accuracy of these algorithms. Our numerical examples reveal that recent algebraic multilevel preconditioning solvers can accelerate the computation of a large-scale eigenvalue problem corresponding to the Anderson model of localization
by several orders of magnitude
GMRES implementations and residual smoothing techniques for solving ill-posed linear systems
AbstractThere are verities of useful Krylov subspace methods to solve nonsymmetric linear system of equations. GMRES is one of the best Krylov solvers with several different variants to solve large sparse linear systems. Any GMRES implementation has some advantages. As the solution of ill-posed problems are important. In this paper, some GMRES variants are discussed and applied to solve these kinds of problems. Residual smoothing techniques are efficient ways to accelerate the convergence speed of some iterative methods like CG variants. At the end of this paper, some residual smoothing techniques are applied for different GMRES methods to test the influence of these techniques on GMRES implementations
- ā¦