Université Bordeaux I, ENSEIRB-MATMECA Année 2012-2013
U.F.R de Mathématiques et Informatique
PROJET de Programmation CHP
Calcul Parallèle
Partie I: Maillages non structurés
• Faire un tutorial de l’utilisation des logiciels scotch et metis pour le partitionnement
de maillages non structurés (illustrer le document avec des figures de maillages);
• Décrire et expliquer en quoi consiste le parallélisme dans le code éléments finis (EF)
fourni en précisant l’équation résolue par le code;
• Discuter les courbes de speed-up et d’efficacité obtenues pour ce code EF avec scotch
et metis.
Partie II: maillage cartésien structuré
Résolution de l’équation de la chaleur instationnaire
On reprend le problème de l’an dernier. On se place dans le domaine [0, Lx ] × [0, Ly ] de IR2
dans lequel on résoud l’équation de la chaleur:
Γ0
Γ1 Γ1
∂t u(x, y, t) − D∆u(x, y, t) = f (x, y, t)
(1) u|
Γ0 = g(x, y, t) Γ0
u| = h(x, y, t)
Γ1
Conditions de bord pour les cas tests de validation
On utilisera les cas test suivants pour valider le travail avec Lx = Ly = 1, D = 1:
• Les solutions stationnaires résultant des conditions suivantes
(2) f = 2 ∗ (x − x2 + y − y 2 ) g=0 h=0
Puis
(3) f = sin(x) + cos(y) g = sin(x) + cos(y) h = sin(x) + cos(y)
1
• La solution instationnaire périodique résultant des conditions suivantes
2
Ly
x− L2x
2 − y− π
(4) f = e−( ) e 2
cos( t) g=0 h=1
2
Implémentation informatique - Coder 2 méthodes de Décomposition de Domaine
basées sur l’algorithme de Schwarz :
1. On ne stockera pas la matrice au complet. On développera une routine faisant le produit
matrice vecteur particulièrement adaptée à la structure de la matrice;
2. On proposera 3 méthodes itératives de résolution du système linéaire (Jacobi, Gauss-Seidel,
Gradient Conjugué);
3. Coder la méthode de Schwarz additive.
4. Coder une méthode de Schwarz multiplicative.
Il faudra fournir les documents suivants:
1. Un rapport contenant
• La comparaison des 3 méthodes: le parallélisme d’opérateur (année dernière), la
décomposition de domaine avec Schwarz additif, et la décomposition de domaine
avec Schwarz multiplicatif.
• Les courbes de convergence des différentes méthodes.
• Une description détaillée des communications réalisées.
• Les courbes du temps de calcul en fonction du nombre de processeurs avec le ratio
entre le temps de communication et le temps de calcul.
• L’analyse des courbes de Speed-up et d’efficacité de vos programmes.
• Pour la décomposition de domaine (méthodes de Schwarz)
– Faire la comparaison des différents solveurs choisis.
– Tracer les courbes du nombre d’itérations de Schwarz en fonction du nombre
d’inconnues, du solveur, et autres paramètres que vous choisirez.
– Etudier l’influence du recouvrement, et si possible des conditions de bord trans-
mises.
• Vos commentaires et conclusions.
2. Les codes // documentés et commentés, avec les commandes de compilation.
