Discussion :

[Yin Zhen]

Je souhaite faire passer un code d'une parallélisation MPI à une parallélisation hybride MPI + openMP, et cela m'amène à me poser la question suivante :

Actuellement, je manipule des tableaux du style : r(1:n,1:d), où n est très grand et où d représente la dimension (2 ou 3). Les calculs sont exécutés identiquement sur les indices de 1 à n, par exemple : r(:,1)=a(:,1)+b(:,1).

Afin d'obtenir de bonnes performances avec l'openMP, ai-je intérêt à passer de ce schéma à un schéma où les tableaux vont être du style : r(1:d,1:n), tau(1:d,1:d,1:n), l'indice "de grande taille" étant rejeté à la fin ?

Merci d'avance.

[Sylvain Bergeron]

En Fortran, les tableaux sont organisés comme suit en mémoire :

(1,1), (2,1), (3,1), (n,1), (1,2), (2,2), ... (n,m)

Pour que le traitement soit efficace, tu dois donc maintenir les indices tels qu'ils sont.