Discussion :

[locate]

Bonjour

Je souhaite paralléliser une résolution d'un système tridiagonale par MPI en utilisant 'Petsc' (qui peut appeler des librairies externes comme superLU,....) .

Je souhaite savoir quel est le meilleur solveur à utiliser et ses paramètres associés (au niveau vitesse )

La matrice A du système Ax=b est tridiagonale symétrique

Merci

[plumedesiles]

Salut,

Un bon début pour la résolution de systèmes linéaires http://jmblanc.developpez.com/algori...?page=sommaire

Toutes les méthodes n'y sont pas mais l'auteur devrait s'y ateler.
Pour répondre à ta question, la matrice étant sysmétrique, il faudrait voir si elle est définie positive. Si elle est, tu peux utiliser la méthode du gradient conjugué.

Il faut que tu saches quelle précision tu as besoin et quelle place mémoire ta future matrice occupera car cést cela qui derminera le choix de ta méthode numérique

[locate]

Salut

Merci pour la réponse.