Discussion :

[maxi_simo]

Bonjour,
je suis débutant sur Matlab et j'ai réussi à écrire la somme des matrice ci dessous en utilisant les boucle for. est ce qu'il est possible d'optimiser ce calcul ?
Merci d'avance pour votre aide :

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
np=10;N=1000;
cm=0;Nn=0;
Bchap=[0 0;0 0];
Cchap=[0 0;0 0];
for i = 1:N
    Bchap1(:,:)=[0 0;0 0] ;
    Cchap1(:,:)=[0 0;0 0] ;
    for s = 1:np
        for j = 1:ne(1,s)
            for k = 1:ne(1,s)
                if k==j
                    Bchap1(:,:)=Bchap1(:,:)+(w(1,s).^2)*([x(j,s,i) y(j,s,i)]' * [x(j,s,i) y(j,s,i)]);
                end
                if k~=j
                    Cchap1(:,:)=Cchap1(:,:)+(w(1,s).^2)*([x(j,s,i) y(j,s,i)]' * [x(k,s,i) y(k,s,i)]);
                end
            end
        end
    end
    Bchap=Bchap+Bchap1(:,:)/sum(ne);
    Cchap=Cchap+Cchap1(:,:)/sum(ne);
    i
end

[shaiHulud]

Bonjour,

Les boucles sont relativement optimisées dans les dernières release, et elles évitent des éventuelles out-of-memory error. Commence peut être par parraléliser la boucle extérieur (doc parfor).

Difficile en tout cas à la simple vue du code de proposer une version vectorisée

[Merel]

Bonjour,

j'ajouterai juste une précision au message de VV33D.

Tu peux utiliser la fonction profile de matlab pour savoir quelle ligne prend du temps à l'execution. Ca donne déjà une idée...

Code :

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

1
2
3
4
 
profile on
% Code
profile viewer

[Jerome Briot]

Une des deux conditions devrait pouvoir être sortie de la boucle sur k :

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
np=10;N=1000;
cm=0;Nn=0;
Bchap=zeros(2);
Cchap=zeros(2);
Bchap1=zeros(2);
Cchap1=zeros(2);
for i = 1:N
    Bchap1(:)=0;
    Cchap1(:)=0;
    for s = 1:np
        for j = 1:ne(1,s)            
            Bchap1(:,:)=Bchap1(:,:)+(w(1,s).^2)*([x(j,s,i) ; y(j,s,i)] * [x(j,s,i) y(j,s,i)]);
            for k = 1:ne(1,s)
                if k~=j
                    Cchap1(:,:)=Cchap1(:,:)+(w(1,s).^2)*([x(j,s,i) ; y(j,s,i)] * [x(k,s,i) y(k,s,i)]);
                end
            end
        end
    end
    Bchap=Bchap+Bchap1(:,:)/sum(ne);
    Cchap=Cchap+Cchap1(:,:)/sum(ne);
end

Sinon, il nous faudrait des valeurs pour les variables w, x, y, np et ne