Discussion :

[pdufranc]

Bonjour
Je dois résoudre Ax =b
avec A et b en nombre complexe

A1 A2 A3 0 0 0
A4 A5 A6 0 0 0
A= A7 A8 A9 0 0 0
0 0 0 A10 A11 A12
0 0 0 A13 A14 A15
0 0 0 A16 A17 A18

B1
B2
b= B3
B4
B5
B6
Ma résolution doit se faire en fortran et je peux utiliser la librairie MKL de Intel fortran, et/ou faire de l'OPENMP

L'objectif est d'être le plus rapide possible.
Est ce que quelqu'un aurait qqc chose sur étagère ?
Merci par avance

[FR119492]

Deux remarques:

1. Pour résoudre un système linéaire, l'utilisation de la matrice inverse n'est pas optimale.
2. Je pense que les bibliothèques LinPack et LAPack sont ce qu'il y a de mieux.

Jean-Marc Blanc

[Ehouarn]

Bonjour,

Vu sa structure, ton système 6x6 est plutôt deux systèmes 3x3 indépendants, non?
Pour des systèmes si petits on peut a priori juste appliquer la solution telle qu'exprimée par résolution analytique.
Sinon (ou de plus), Intel inclue la libraririe Lapack, y'a plus qu'à utiliser la routine concernée (et pour être optimal en l'appelant deux fois, une sur chaque système 3x3, pas sur le système 6x6 entier).

Ehouarn

[genteur slayer]

Bonjour,

Je suis d'accord avec Ehouarn, ton système est en fait deux système 3x3 indépendants.

cela dit, si tu veux être rapide, dans la librairie MKL d'intel, il y a déja tout ce qu'il faut, en particulier les librairie blas et lapack, ces solveurs sont d'ailleurs déjà optimisées et parallélisées....