+ Répondre à la discussion
Affichage des résultats 1 à 9 sur 9
  1. #1
    Membre du Club
    Inscrit en
    septembre 2007
    Messages
    99
    Détails du profil
    Informations forums :
    Inscription : septembre 2007
    Messages : 99
    Points : 47
    Points
    47

    Par défaut Inverser une matrice en parallele quelles librairies ?

    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.

  2. #2
    Membre habitué
    Inscrit en
    mai 2010
    Messages
    122
    Détails du profil
    Informations forums :
    Inscription : mai 2010
    Messages : 122
    Points : 132
    Points
    132

    Par défaut

    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

  3. #3
    Membre du Club
    Inscrit en
    septembre 2007
    Messages
    99
    Détails du profil
    Informations forums :
    Inscription : septembre 2007
    Messages : 99
    Points : 47
    Points
    47

    Par défaut

    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.

  4. #4
    Dut
    Dut est déconnecté
    Rédacteur/Modérateur

    Avatar de Dut
    Homme Profil pro Jérôme Briot
    Ingénieur hospitalier
    Inscrit en
    novembre 2006
    Messages
    17 805
    Détails du profil
    Informations personnelles :
    Nom : Homme Jérôme Briot
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Ingénieur hospitalier
    Secteur : Santé

    Informations forums :
    Inscription : novembre 2006
    Messages : 17 805
    Points : 49 577
    Points
    49 577

    Par défaut

    Je ne code pas en Fortran mais tu peux jeter un œil à PLASMA
    Sur mon blog : Installer GNU Octave 3.8.1 sur Linux Xubuntu 14.10

    Ma page personnelle avec mes contributions MATLAB (R2011b/R2014b - Windows, Linux et Mac) et Hardware

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  5. #5
    Rédacteur/Modérateur

    Homme Profil pro Jean-Marc Blanc
    Comme retraité, des masses
    Inscrit en
    avril 2007
    Messages
    2 961
    Détails du profil
    Informations personnelles :
    Nom : Homme Jean-Marc Blanc
    Âge : 74
    Localisation : Suisse

    Informations professionnelles :
    Activité : Comme retraité, des masses
    Secteur : Industrie

    Informations forums :
    Inscription : avril 2007
    Messages : 2 961
    Points : 5 065
    Points
    5 065

    Par défaut

    Salut!
    je dois inverser de grosses matrices
    A quoi cela sert-il d'inverser ces matrices?
    Jean-Marc Blanc
    Calcul numérique de processus industriels
    Formation, conseil, développement

    Point n'est besoin d'espérer pour entreprendre, ni de réussir pour persévérer. (Guillaume le Taiseux)

  6. #6
    Membre éprouvé
    Inscrit en
    mars 2007
    Messages
    372
    Détails du profil
    Informations forums :
    Inscription : mars 2007
    Messages : 372
    Points : 491
    Points
    491

    Par défaut

    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.

  7. #7
    Membre du Club
    Inscrit en
    septembre 2007
    Messages
    99
    Détails du profil
    Informations forums :
    Inscription : septembre 2007
    Messages : 99
    Points : 47
    Points
    47

    Par défaut

    Bonjour,
    merci pour vos reponses. Pour repondre dans l'ordre :

    Dut, je viens de compiler plasma sur ma machine et de feuilleter rapidement la doc. Ca m'a l'air au poil ! Malgre tout, je ne vois rien concernant le link de ces libs lors de la compilation. Je vais creuser un peu.

    Pour repondre a Jean-Marc j'ai ecrit un code d'elements finis qui utilise un algorithme de type Newton-Raphson pour resoudre une equation differentielle non-lineaire. Je dois donc a resoudre le systeme [A]{X}={Y}, avec [A] une grosse matrice, et {.} des vecteurs. Inversion de matrice est donc un abus de language de ma part. Il s'agit de la resolution du systeme suivant [A]{X}={Y}

    Enfin Ehouarn j'ai en effet vu ce type de solution, notament avec ScaLapack. Cependant je ne suis vraiment pas un specialiste en la matiere, et j'avoue que ce type de chose depasse tres largement mes competences...

    Merci beaucoup,
    Yann.

  8. #8
    Membre régulier Avatar de moomba
    Inscrit en
    avril 2007
    Messages
    134
    Détails du profil
    Informations forums :
    Inscription : avril 2007
    Messages : 134
    Points : 77
    Points
    77

    Par défaut

    Simple petit ajout : si tu ne cherche qu'à résoudre le système, une méthode itérative type PCG ou BICGSTAB semble plus indiqué, surtout si ta matrice est grosse en mémoire. Tout dépend de l'allure de ta matrice, mais avec des éléments fini ça doit être bon...
    Ces méthodes sont facilement parallélisable sur des architectures parallèles avec MPI ou OpenMP.
    "Celui qui à le pouvoir de faire le mal, mais qui ne le fait pas, celui là est le prince de l'univers." (shakespeare)

  9. #9
    Invité de passage
    Inscrit en
    septembre 2012
    Messages
    6
    Détails du profil
    Informations forums :
    Inscription : septembre 2012
    Messages : 6
    Points : 2
    Points
    2

    Par défaut

    Citation Envoyé par yann_m Voir le message
    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

+ Répondre à la discussion
Cette discussion est résolue.

Liens sociaux

Règles de messages

  • Vous ne pouvez pas créer de nouvelles discussions
  • Vous ne pouvez pas envoyer des réponses
  • Vous ne pouvez pas envoyer des pièces jointes
  • Vous ne pouvez pas modifier vos messages
  •