Discussion :

[domdom974]

Bonjour,

Je suis a la recherche d'un solveur pour le systeme lineaire AX = B

avec
A matrice carree complexe
B et X matrices rectangulaires complexes (pas un vecteur)

sans pour autant passer par des grosses librairies du style Lapack...

Si quelqu'un en a un je suis tout oui!!!

Merci d'avance,

Domdom974

[domdom974]

Bon, j'ai trouve une subroutine qui fait ca, a condition que A soit triangulaire.

Donc je reformule ma requete:

je cherche soit un solveur pour systeme de matrices complexes comme precedemment soit un petit code qui trigonalise une matrice complexe.

Si quelqu'un en a un, je veux bien.


Merci d'avance,


Domdom974

[genteur slayer]

et ce n'est pas plus simple de convertir X et B en vecteur (en transformant A évidement pour garder le même problème) juste avant ta résolution puis de les remettre en rectangle ensuite? remarque cela ne change pas trop le schmilblick, mais tu trouvera peut-être plus de routines de résolution....

[FR119492]

Salut!

des grosses librairies du style Lapack

La bibliothèque Lapack est grosse parce qu'elle permet de résoudre une masse de problèmes différents, mais chaque sous-programme est assez simple à utiliser. Si tu veux encore plus simple, va voir dans Linpack.
Jean-Marc Blanc

[domdom974]

Merci pour les remarques.

Je vais essayer de faire quelque chose a partir de l'idee de genteur slayer.

Domdom974