| 1234567891011121314151617181920212223242526272829 |
- % résolution de Ax = b avec du multithreading
- n = 20000; % taille de la matrice n x n
- P = randn(n);
- A = P' + P; % A matrice symétrique
- %A = triu(rand(n)); % A triangulaire
- x = ones(n,1); % la solution est le vecteur unité
- b = A*x; % calcul du rhs
- fprintf('\n\n*** resolution systeme lineaire ***\n');
- fprintf('\n taille du systeme, n = %d \n\n',n);
- tic % top horloge
- y1 = A\b; % calcul de la solution (en principe vec unité!)
- time= toc; % top horloge
- fprintf('\n temps execution mldivide: %f \n',time);
- opts= struct('SYM', true);
- %opts= struct('UT', true);
- tic
- y2 = linsolve(A,b,opts);
- %y2 = linsolve(A,b);
- time= toc; % top horloge
- fprintf('\n temps execution linsolve: %f \n',time);
- % merci de quitter matlab pour libérer les jetons de licence.
- quit
|
