Discussion :

[yann_m]

Bonjour a vous,

je suis actuellement entrain de paralleliser un code dans lequel je dois inverser de grosses matrices.
Aussi, aurais-je aime savoir si quelqu'un peut me renseigner sur les differents types de librairies open sources concuent a cet effet, et utilisant de preference OpenMP.
Enfin ou se procurer des references (liens ou documents) contenant des exemples de mises en oeuvres de l'utilisation de tels librairies ?

Merci par avance.

Bien cordialement,
Yann.

[Marlan]

Salut,

Tu peux par exemple utiliser OPENMP à l'appel de ta subroutine et non en son sein. Ainsi la parallélisation sera efficace, sans que tu n'aies à te poser la question du fait que la subroutine soit parallèle ou non.

ie : Ne pas paralléliser dans la subroutine mais faire entrer en argument des portions de tableau associés à chaque THREADS

Ex : Si ma subroutine s'appelle INV_MAT(A,B) avec A la matrice de départ et B la matrice d'arriver, il suffit de faire rentrer une matrice partagée entre les différents threads pour B.

Note : cette astuce ne marche que dans le cas de résolutions explicites de matrices [les coefficients de B s'expriment directement et uniquement en fonction de ceux de A] (pas avec des méthodes de pivots par exemple)

En espérant avoir pu t'aider,

Marlan

[yann_m]

Bonjour,

merci pour ta réponse. Ca aurait été en effet absolument parfait, puisque ce que tu me suggères correspond exactement à ce que je souhaite faire :

faire entrer en argument des portions de tableau associés à chaque THREADS

Cependant outre les inversions de matrices, j'ai également des décompositions LU à réaliser, et par souci d'optimisation je souhaite pouvoir utiliser une librairie de type Lapack (parallèle) fonctionnant avec openMP.

Merci beaucoup,
Yann.

[Jerome Briot]

Je ne code pas en Fortran mais tu peux jeter un œil à PLASMA

[FR119492]

Salut!

je dois inverser de grosses matrices

A quoi cela sert-il d'inverser ces matrices?
Jean-Marc Blanc

[Ehouarn]

Bonjour,

Il y a PLAPACK en contexte mémoire distribuée (MPI) qui semble pas mal. Peut-être que l'utiliser en mode mixte OpenMP/MPI est possible (moyennant sans doute une certaine dose d'huile de coude...).

Bonne continuation.

[Bast24]

Bonjour,

merci pour ta réponse. Ca aurait été en effet absolument parfait, puisque ce que tu me suggères correspond exactement à ce que je souhaite faire :

faire entrer en argument des portions de tableau associés à chaque THREADS

Cependant outre les inversions de matrices, j'ai également des décompositions LU à réaliser, et par souci d'optimisation je souhaite pouvoir utiliser une librairie de type Lapack (parallèle) fonctionnant avec openMP.

Merci beaucoup,
Yann.

Lapack et scalapack fonctionnent bien avec MPI, je sais pas trop avec OpenMp