C : Diagonalisation d'une matrice complexe

[invite62a84de7]

Bonjour à tous,

Je cherche à diagonaliser une matrice complexe qui n'a pas de propriété particulière (il se pourrait bien qu'elle soit unitaire, mais ce n'est pas sûr). J'ai pas mal cherché, mais j'ai l'impression qu'il n'y a pas en C de librairie standard qui réalise une telle opération. Par exemple, GSL ne semble pas savoir le faire. J'ai l'impression qu'il existe une librairie standard qui le fait en Fortran, LAPACK, mais il me semble improbable qu'il n'y ait pas d'analogue en C.

Merci beaucoup de vos suggestions,
oliv (qui n'a pas envie de se mettre au Fortran )
-----

[inviteb9f49292]

CLAPACK est une version C de LAPACK (de mémoire), mais tu as aussi ATLAS qui fait peut-être ça, ou bien suivre "numerical recipes" (google devrait t'y ammener), même si ce n'est paut-être pas ce qu'il y a de plus efficace

[invite62a84de7]

Bonjour,

Merci beaucoup de ta réponse. Malheureusement, ATLAS ne fait pas, mon édition de Numerical Recipes in C ne donne pas d'algorithme général de diagonalisation d'une matrice complexe, et je trouve CLAPACK difficile à installer - en plus, je ne suis pas sûr que pour des manipulations compliquées, une traduction de Fortran en C soit très optimale en temps, et encore moins en mémoire (et sur des matrices carrées de taille 2^13, ça commence à compter!). Je devrais pouvoir me débrouiller avec une décomposition en valeurs singulières, mais si quelqu'un a des lumières, je prends!

oliv