numerical recipes - erreurs ?

Question

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

turnerom · Answer

Salut

je ne connais pas les numerical recipes, mais je te conseil d'utiliser nt2 disponible sur sourceforge

TuRn3r

benlustr · Answer

ok, je vais regarder :)

Merci

chouchou79 · Answer

j'ai besoin d'un code source qui calcule un intégrale à deux dimensions

Numerical recipes -> erreurs ?

3 réponses

Votre réponse

Discussions similaires