Multipreconditioning for nonsymmetric problems: the case of orthomin and biCG
Résumé
Preconditioned Krylov subspace methods are powerful tools for solving linear systems but sometimes they converge very slowly, and often after a long stagnation. A natural way to fix this is by enlarging the space in which the solution is computed at each iteration. Following this idea, we propose in this note two multipreconditioned algorithms: multipreconditioned orthomin and multipreconditioned biCG which aim at solving general nonsingular linear systems in a small number of iterations. After describing the algorithms, we illustrate their behaviour on systems arising from the FETI domain decomposition method, where in order to enlarge the search space, each local component in the usual preconditioner is kept as a separate preconditioner.
Les solveurs de Krylov préconditionnés sont des outils performants pour la résolution de systèmes linéaires. Il arrive cependant qu'ils convergent lentement, souvent après une phase de stagnation. Une manière naturelle de remédier à cette situation est d'agrandir l'espace dans lequel on cherche la solution à chaque itération. En suivant cette idée, nous proposons dans cette note deux algorithmes multipréconditionnés : orthomin multipréconditionné (MPorthomin) et biCG multipréconditionné (MPbiCG) avec l'objectif de résoudre des systèmes linéaires généraux en un petit nombre d'itérations. Après avoir décrit les nouveaux algorithmes, nous illustrons leur comportement sur des systèmes linéaires issus de la méthode de décomposition de domaine FETI, où au lieu d'appliquer le préconditionneur habituel, on fait agir chaque contribution locale au préconditionneur séparément.
Origine : Fichiers produits par l'(les) auteur(s)
Loading...