Quel est l'état actuel des préconditionneurs polynomiaux?


15

Je me demande ce qui est arrivé aux préconditionneurs polynomiaux. Je m'intéresse à eux, car ils semblent relativement élégants d'un point de vue mathématique, mais d'après ce que j'ai lu dans les enquêtes sur les méthodes de krylov, ils fonctionnent généralement très mal en tant que préconditionneurs. Pour reprendre les mots de Saad et van der Host, «L'intérêt actuel pour ces techniques a pratiquement disparu» (ici) . Néanmoins, des calculs multicœurs et GPU ont été utilisés récemment.

Quelqu'un peut-il me dire ou plutôt m'expliquer dans quels contextes ces méthodes sont toujours vivantes et où trouver un bon aperçu de l'état actuel de la technique?


Un article récent sur l'arxiv ( arxiv.org/pdf/1806.08020.pdf ) étudie les préconditionneurs polynomiaux d'Arnoldi. En particulier, ils le testent sur une variété de problèmes et obtiennent une bonne accélération. Ils concluent que la réduction des opérations vectorielles due au préconditionnement polynomial "est très prometteuse pour le calcul de valeurs propres évitant la communication sur des ordinateurs à hautes performances". Je ne suis pas l'un des auteurs.
amarney

Réponses:


12

Pour fonctionner raisonnablement, les préconditionneurs polynomiaux ont besoin d'estimations spectrales assez précises. Pour les problèmes elliptiques mal conditionnés, les plus petites valeurs propres sont généralement séparées de sorte que des méthodes comme Tchebychev sont loin d'être optimales. La propriété la plus intéressante des méthodes polynomiales est qu'elles ne nécessitent aucun produit interne.

Il est en fait assez populaire d'utiliser des lisseurs polynomiaux en multigrille. La principale différence avec un préconditionneur est que le lisseur n'est censé cibler qu'une partie du spectre. Un lisseur polynomial est actuellement la valeur par défaut dans le multigrille de PETSc, par exemple. Voir aussi Adams et al, Parallel multigrid smoother: polynomial versus Gauss-Seidel (2003) pour une comparaison.

Les préconditionneurs polynomiaux peuvent être utilisés uniquement pour réduire la fréquence des réductions. Bien qu'ils doivent être recalibrés pour chaque matrice, les économies peuvent être importantes sur le matériel dans lequel les réductions sont coûteuses (commun sur les grands supercalculateurs). Voir McInnes, Smith, Zhang et Mills, Hierarchical and Nested Krylov Methods for Extreme-Scale Computing (2012) pour plus d'informations.

En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.