Discussion :

[asmasbai]

Bonjour,
J'ai remarqué une grande lenteur lors de l'execution de mon code et ce au niveau de la construction des m equations:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
 
m=500;
%CONSTRUIRE LES EQUATIONS
eq=sym('eq',[m,1]);
    Y=sym('Y',[m,1]);
    for a=1:m
        eq(a)=0;
        for i=1:n
            for j=1:n
                eq(a) = (U(i,j,a))* (Y(a)^ P(i,j,a)) + eq(a);
            end
        end
        eq(a)=eq(a)-V(a);
      end
%RESOUDRE LES EQUATIONS
y=zeros(1,m);
    for a=1:m
        solution=solve(eq(a),Y(a),'PrincipalValue',true); 
       if isempty(solution)
           y(a)=0;
       else
           y(a)=solution;
       end
    end

Pourriez vous m'aider? car je dois encore augmenter le nombre d'equations de 500 vers 1000

Merci d'Avance

[Kratein]

Tu peux utiliser la methode parfor.

Elle va te permettre d'utiliser des threads pour diminuer le temps de traitement.

[Dranen]

L'utilisation de parfor peut être utile mais avant même d'utiliser cela il est préférable d'optimiser le code. Matlab préfère nettement quand tu présente le calcul matriciellement.
Ainsi en modifiant ton code de la façon suivant j’obtiens une multiplication de la vitesse d’exécution par environ 30 sur un test rapide.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
for a=1:m
    eq2(a) = sum(sum(U(:,:,a).*(Y(a).^P(:,:,a))));
    eq2(a)=eq2(a)-V(a);
end

Il y a certainement moyen de faire disparaître la dernière boucle mais je ne suis pas sûr du comportement de Y.^P quand Y est un vecteur et P un tableau à trois dimensions. Je te laisse faire les tests.