Je suis à la recherche d'un algorithme permettant le calcul des valeurs
propres et vecteurs propres d'une matrice symétrique (symetric eigen
problem). Cet algorithme devrait avoir les caractéristiques suivantes:
- être implémenté en C
- être parallèle (utilisation multithread/multiprocesseurs)
- flexible quand au nombre de vecteurs retournés (je m'intéresse soit à
tout les vecteurs propres, soit à ceux ayant les plus grandes valeurs
propres).
J'ai exploré PLAPACK et ScaLAPACK, mais ces produits sont difficiles à
installer et peu portables. La méthode powerEigen me semble prometteuse,
mais je n'ai pas trouvé d'implémentation en parallèle.
Avez vous des idées à ce sujet?
Merci,
Olivier