Bonjour à tous, j'ai besoin pour un projet d'un algorithme permettant d'obtenir les valeurs/vecteurs propres d'une matrice tridiagonale.
Je me suis tourné vers numerical recipes in C, et plus particulièrement la fonction tqli (ici : http://www.library.cornell.edu/nr/bookcpdf/c11-3.pdf , pages 6 et 7 du document)
Or, certains éléments de ce programme me semblent bizarres :
(ne vous souciez pas de pythag et SIGN)
-> Les indices des tableaux commencent à 1 ????
-> On ne sort jamais de la boucle do /while, sauf dans des cas particuliers (matrice de dimension n*n, n=2, ou si certains éléments de la subdiagonale sont nuls, ce qui n'est pas mon cas....).
Bref, que pensez vous de ce programme, ai-je une chance de le faire fonctionner ?
Et que pensez-vous de Numerical recipes en général ? Est plutôt fiable ?
Merci à vous