III.C.4 Code Helmholtz bidimensionnel dans le vide.

Sur le cas de la section III.C.4, on donne les temps de calcul des parties essentielles du programme. Nous remarquerons que le code bidimensionnel n'est pas parfaitement optimisé. La performance globale du programme est de l'ordre de 110 MFlops^III.C1.

Tableau III.C.3: Performances du code Helmholtz 2D.

Tâche	Temps (s)	MFlops
10	4.15	143.7
2	1.15	0.1
9	0.331	0.2
5	0.43	125.1
6	0.186	15.6
7	0.0564	123.7
8	0.0316	98.9

Le temps de calcul du code dépend donc par ordre décroissant des procédures suivantes.

i)

L'algorithme itératif (I.2.31) qui consiste à effectuer des produits matrice-vecteur. En pratique, un nombre d'itérations de 600 s'est révélé suffisant.

ii)

Calculer les matrices D et C.

iii)

Inverser les matrices D_k (de petite taille, inférieure, en pratique, à 15 x 15).