Numerical study of a double preconditioning strategy

• Cédric Doucet ^[1] ; Isabelle Charpentier ^[3] ; Jean-Louis Coulomb ^[2] ; Christophe Guérin ^[4] ; Yann Le Floch ^[5] ; Gerard Meunier ^[2]
  1. [1] Laboratoire Jean Kuntzmann

     Laboratoire Jean Kuntzmann

     Arrondissement de Grenoble, Francia

     
  2. [2] Laboratoire de Génie Electrique de Grenoble

     Laboratoire de Génie Electrique de Grenoble

     Arrondissement de Grenoble, Francia

     
  3. [3] Université de Metz
  4. [4] Cedrat SA
  5. [5] CEDRAT

  Mostrar afiliaciones +
• Localización: Compel: International journal for computation and mathematics in electrical and electronic engineering, ISSN 0332-1649, Vol. 27, Nº 4 (Special Issue: ISEF 2007), 2008, págs. 897-903
• Idioma: inglés
• Enlaces
  • Texto completo
• Resumen

  • Purpose – The aim of this paper is to accelerate the convergence of iterative methods on ill‐conditioned linear systems of equations.

    Design/methodology/approach – First a brief numerical analysis is given of left preconditioners on ill‐conditioned linear systems of equations. From this result, it is deduced that a double preconditioning approach may be better. Then, a double preconditioner based on an iterative diagonal scaling method and an incomplete factorization method is proposed. The efficiency of this approach is illustrated on two finite element models produced by computational electromagnetism.

    Findings – The double preconditioning approach is efficient for 2D and 3D finite element problems. The bi‐conjugate gradient algorithm always converges when it is double preconditioned. This is not the case when a simple incomplete factorization method is applied. Furthermore, when the two preconditioning techniques lead to the convergence of the iterative solving method, the double preconditioner significantly reduces the number of iterations in comparison with the simple preconditioner. On the proposed 2D problem, the speed‐up is between 6 and 32. On the proposed 3D problem, the speed‐up is between 13 and 20. Finally, the approach seems to reduce the growth of the condition number when higher‐order finite elements are used.

    Research limitations/implications – The paper proposes a particular double preconditioning approach which can be applied to any invertible linear system of equations. A numerical evaluation on a singular linear system is also provided but no proof or analysis of stability is given for this case.

    Originality/value – The paper presents a new preconditioning technique based on the combination of two very simple and elementary methods: a diagonal scaling method and an incomplete factorization process. Acceleration obtained from this approach is quite impressive.