Results 1  10
of
86
Analysis of iterative methods for the forcebased quasicontinuum method. manuscript
, 2009
"... Abstract. Forcebased atomisticcontinuum hybrid methods are the only known pointwise consistent methods for coupling a general atomistic model to a finite element continuum model. For this reason, and due to their algorithmic simplicity, forcebased coupling methods have become a popular class of a ..."
Abstract

Cited by 31 (18 self)
 Add to MetaCart
Abstract. Forcebased atomisticcontinuum hybrid methods are the only known pointwise consistent methods for coupling a general atomistic model to a finite element continuum model. For this reason, and due to their algorithmic simplicity, forcebased coupling methods have become a popular class of atomisticcontinuum hybrid models as well as other types of multiphysics models. However, the recently discovered unusual stability properties of the linearized forcebased quasicontinuum (QCF) approximation, especially its indefiniteness, present a challenge to the development of efficient and reliable iterative methods. We present analytic and computational results for the generalized minimal residual (GMRES) solution of the linearized QCF equilibrium equations. We show that the GMRES method accurately reproduces the stability of the forcebased approximation and conclude that an appropriately preconditioned GMRES method results in a reliable and efficient solution method. 1.
Computing and deflating eigenvalues while solving multiple right hand side linear systems with an application to quantum chromodynamics
, 2008
"... Abstract. We present a new algorithm that computes eigenvalues and eigenvectors of a Hermitian positive definite matrix while solving a linear system of equations with Conjugate Gradient (CG). Traditionally, all the CG iteration vectors could be saved and recombined through the eigenvectors of the t ..."
Abstract

Cited by 21 (2 self)
 Add to MetaCart
(Show Context)
Abstract. We present a new algorithm that computes eigenvalues and eigenvectors of a Hermitian positive definite matrix while solving a linear system of equations with Conjugate Gradient (CG). Traditionally, all the CG iteration vectors could be saved and recombined through the eigenvectors of the tridiagonal projection matrix, which is equivalent theoretically to unrestarted Lanczos. Our algorithm capitalizes on the iteration vectors produced by CG to update only a small window of about ten vectors that approximate the eigenvectors. While this window is restarted in a locally optimal way, the CG algorithm for the linear system is unaffected. Yet, in all our experiments, this small window converges to the required eigenvectors at a rate identical to unrestarted Lanczos. After the solution of the linear system, eigenvectors that have not accurately converged can be improved in an incremental fashion by solving additional linear systems. In this case, eigenvectors identified in earlier systems can be used to deflate, and thus accelerate, the convergence of subsequent systems. We have used this algorithm with excellent results in lattice QCD applications, where hundreds of right hand sides may be needed. Specifically, about 70 eigenvectors are obtained to full accuracy after solving 24 right hand sides. Deflating these from the large number of subsequent right hand sides removes the dreaded critical slowdown, where the conditioning of the matrix increases as the quark mass reaches a critical value. Our experiments show almost a constant number of iterations for our method, regardless of quark mass, and speedups of 8 over original CG for light quark masses.
THE EXPONENTIALLY CONVERGENT TRAPEZOIDAL RULE
"... Abstract. It is well known that the trapezoidal rule converges geometrically when applied to analytic functions on periodic intervals or the real line. The mathematics and history of this phenomenon are reviewed and it is shown that far from being a curiosity, it is linked with computational methods ..."
Abstract

Cited by 17 (3 self)
 Add to MetaCart
(Show Context)
Abstract. It is well known that the trapezoidal rule converges geometrically when applied to analytic functions on periodic intervals or the real line. The mathematics and history of this phenomenon are reviewed and it is shown that far from being a curiosity, it is linked with computational methods all across scientific computing, including algorithms related to inverse Laplace transforms, special functions, complex analysis, rational approximation, integral equations, and the computation of functions and eigenvalues of matrices and operators.
INTERPRETING IDR AS A PETROVGALERKIN METHOD
, 2009
"... The IDR method of Sonneveld and van Gijzen [SIAM J. Sci. Comput., 31:1035–1062, 2008] is shown to be a PetrovGalerkin (projection) method with a particular choice of left Krylov subspaces; these left subspaces are rational Krylov spaces. Consequently, other methods, such as BiCGStab and ML(s)BiCG ..."
Abstract

Cited by 13 (2 self)
 Add to MetaCart
The IDR method of Sonneveld and van Gijzen [SIAM J. Sci. Comput., 31:1035–1062, 2008] is shown to be a PetrovGalerkin (projection) method with a particular choice of left Krylov subspaces; these left subspaces are rational Krylov spaces. Consequently, other methods, such as BiCGStab and ML(s)BiCG, which are mathematically equivalent to some versions of IDR, can also be interpreted as PetrovGalerkin methods. The connection with rational Krylov spaces inspired a new version of IDR, called RitzIDR, where the poles of the rational function are chosen as certain Ritz values. Experiments are presented illustrating the effectiveness of this new version.
Deflated and augmented Krylov subspace methods: A framework for deflated . . .
, 2013
"... We present an extension of the framework of Gaul et al. (SIAM J. Matrix Anal. Appl. 34, 495–518 (2013)) for deflated and augmented Krylov subspace methods satisfying a Galerkin condition to more general Petrov–Galerkin conditions. The main goal is to apply the framework also to the biconjugate gra ..."
Abstract

Cited by 13 (2 self)
 Add to MetaCart
(Show Context)
We present an extension of the framework of Gaul et al. (SIAM J. Matrix Anal. Appl. 34, 495–518 (2013)) for deflated and augmented Krylov subspace methods satisfying a Galerkin condition to more general Petrov–Galerkin conditions. The main goal is to apply the framework also to the biconjugate gradient method (BiCG) and some of its generalizations, including BiCGStab approach does not depend on particular recurrences and thus simplifies the derivation of theoretical results. It easily leads to a variety of realizations by specific algorithms. We do not go into algorithmic details, but we show that for every method there are two different approaches for extending it by augmentation and deflation: one that explicitly takes care of the augmentation space in every step, and one that applies the unchanged basic algorithm to a projected problem but requires a correction step at the end. Both typically generate a Krylov space for a singular operator that is associated with the projected problem. The deflated biconjugate gradient requires two such Krylov spaces, but it also allows us to solve two dual linear systems at once. Deflated Lanczostype product methods fit in our new framework too. The question of how to extract the augmentation and deflation subspace is not addressed here.
Block preconditioning of realvalued iterative algorithms for complex linear systems
, 2008
"... ..."
ON MUTUAL IMPACT OF NUMERICAL LINEAR ALGEBRA AND LARGESCALE OPTIMIZATION WITH FOCUS ON INTERIOR POINT METHODS
, 2008
"... ..."