Bonjour,

Je suis bloqué dans un code que j'ai téléchargé sur le net (http://www.cis.upenn.edu/~jshi/software/) à cause d'une fonction nommé arpackc. Si j'ai bien compris, à partir de maintenant, je peux faire la même chose que cette fonction (calcul des valeurs propres d'une matrice je pense), en utilisant la fonction eigs. Malheureusement cela ne semble pas fonctionner et ne me donne pas la même chose.

Voilà les lignes originales :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
 
arpackc( 'dseupd', rvec, 'A', select, 'A', select, ...
                 d, v, ldv, sigma, ...
                 bmat, int32(n), whch, nev, tol, resid, ncv, ...
                 v, ldv, iparam, ipntr, workd, workl, lworkl, info );
que j'aurai espérer remplacer par une simple fonction :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
d = eigs(matrice)
Sauf que la fonction arpackc apparemment, me permet de calculer des vectors propres sur une matrice rectangulaire en passant par des Lanczos vectors (entre autres...). Ca a l'air bien compliqué et je me demande si je peux y faire quelque chose... et surtout si on peut faire la même chose avec la fonction eigs...

D'avance merci