3. Solveurs et algèbre linéaire#

Concernant les méthodes de numérotation des inconnues, de stockage des matrices assemblées et de résolution des systèmes linéaires sur lesquels s’appuient les différents algorithmes, on dispose aujourd’hui de plusieurs méthodes: solveur direct ou itératif, séquentiel ou parallèle.

En «natif» dans code_aster:

  • Solveurs directs: méthode multi-frontale (parallélisé OpenMP) et méthode \({\mathit{LDL}}^{T}\) .

  • Solveur itératif: gradient conjugué préconditionné.

Ces méthodes sont associées à des algorithmes de renumérotation des degrés de liberté permettant d’optimiser la taille mémoire nécessaire pour stocker les matrices.

En externe, il est possible d’utiliser les solveurs suivants:

  • MUMPS: solveurs direct parallélisé MPI out-of-core;

  • PeTSC: bibliothèque de solveurs itératifs, parallélisées MPI;

Pour le calcul modal: méthodes de Lanczos, Sorensen, les méthodes par puissance inverse, pour les problèmes réels et complexes, symétrique ou non.