previous up next top index
Précédent : Calcul formel Remonter : Algorithmique et applications Suivant : Dynamique moléculaire


Algèbre linéaire sur des structures creuses

Ce travail se fait en collaboration avec le CEA-Grenoble (B. Brun, Direction des Réacteurs Nucléaires).

En calcul scientifique, beaucoup de problèmes se ramènent à la résolution numérique de très grands systèmes linéaires creux (optimisation par méthode de Newton, éléments finis...) qui est l'archétype des problèmes irréguliers. Deux approches complémentaires sont suivies. Nous étudions la parallélisation d'un code itératif (basé sur les méthodes de type Gradient Conjugué, CS-Stab et GMRES). Par ailleurs, les méthodes directes de résolution sont une alternative intéressante pour la parallélisation à cause de l'espace mémoire fourni par les plates-formes à mémoire distribuée. Dans ce cadre, au problème intrinsèque de remplissage de la matrice lors de la factorisation, s'ajoute celui de la régulation de charge en parallèle. Pour ce problème, il est attendu une analyse comparative fine entre régulation dynamique (ordonnancement en ligne avec prise en compte de la localité) et placement statique (tirant parti des connaissances sur le graphe d'élimination).