Discussion :

[Sylvain_F]

Bonjour,
J'ai un programme avec une grande quantité de boucles, son exécution est donc assez lente.
Il y a 3 boucles inbriquées qui traite les 3 dimesions x, y et z, sur des matrices 3D.
J'aimerais pouvoir transformer coder les 3 boucles sous forme de matrices mais il y a un élément qui gène.
Voici une partie du code en question :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
for p=1:np % boucle sur le nombre de points
   x(p)=M(p,2);y(p)=M(p,3);z(p)=M(p,4);f(p)=M(p,5);A(p)=M(p,6);
   w(p)=2*pi*f(p); 
   for ix=1:n
        for iy=1:n
            for iz=1:n
                r(p,ix,iy,iz)=sqrt((X(ix,iy,iz)-x(p))^2+(Y(ix,iy,iz)-y(p))^2+(Z(ix,iy,iz)-z(p))^2);
            end
        end
   end
end

La 2me ligne récupère juste les valeurs de x, y, z, f et A dans une matrice M, ces valeurs sont des vecteurs de dimension (np,1).
X, Y et Z sont de taille (n,n,n).
J'aimerais remplacer les 3 boucles ix, iy et iz par une commande de ce type :

Code :

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

r(p,:,:,:)=sqrt((X-x(p).*ones(n,n,n))^2+(Y-y(p)*ones(n,n,n))^2+(Z-z(p)*ones(n,n,n))^2);

Bien sur cela ne fonctionne pas car il y a un problème de dimension cer r est de dimension (np,n,n,n).
Merci d'avance pour votre aide.

[Jean Dumoncel]

Bonjour,

Tu dois avoir une erreur au niveau des puissances :

Code :

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

r(p,:,:,:)=sqrt((X-x(p).*ones(n,n,n)).^2+(Y-y(p)*ones(n,n,n)).^2+(Z-z(p)*ones(n,n,n)).^2);

Sinon x(p).*ones(n,n,n) peut aussi s'écrire :

Code :

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

repmat(x(p),n,n,n);

[Sylvain_F]

Ok, merci, le code que tu m'as donné fonctionne bien, je ne pensais pas que c'était possible étant donné que les dimensions des deux membres n'était pas identiques.
Par contre repmat ne fonctionne pas (Too many input arguments).
Ensuite j'ai le même type de problème que le premier.
J'aimerais replacer :

Code :

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

1
2
3
4
5
6
7
8
9
for p=1:np
    for ix=1:n
        for iy=1:n
            for iz=1:n
                pres(ix,iy,iz)=A(p)*exp(1i*(w(p)*t-k(p)*r(p,ix,iy,iz)))/(k(p)*r(p,ix,iy,iz));
            end
        end
    end
end

par :

Code :

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

1
2
3
for p=1:np
    pres(:,:,:)=A(p).*exp(1i.*(w(p)*t.*ones(n,n,n)-k(p).*r(p,:,:,:)))./(k(p).*r(p,:,:,:));
end

[Invité]

Bonjour,

x(p), y(p) et z(p) étant de simples scalaires, il n'y a pas besoin de les ajuster aux dimensions de X,Y et Z. Ainsi le code de magelan peut se simplifier en:

Code :

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

r1(p,:,:,:)=sqrt((X-x(p)).^2+(Y-y(p)).^2+(Z-z(p)).^2);

En ce qui concerne REPMAT, au delà de 2 dimensions, il faut fournir un vecteur:

Code :

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

repmat(x(p),[n n n]);

En ce qui concerne cette nouvelle simplification, ne serai-ce pas plutôt pres(p,ix,iy,iz) ?
Quoiqu'il en soit, le principe est le même: dès que tu veux multiplier/diviser membre-à-membre des tableaux, il te faut mettre un point devant:

Code :

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

1
2
A./B
A.*B

[Sylvain_F]

Ok, merci pour ces précisions.
Non il n'y a pas d'indice p sur pres car à la fin de la boucle il y a

Code :

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

sp=sp+real(pres);

Je fais donc la somme de tous les indices, je n'ai donc pas besoin de conserver la valeur pour chaque indice p.

L'erreur que me retourne matlab est :

??? Error using ==> minus
Number of array dimensions must match for binary array op.

Ce qui me semble normal car il y a l'indice p en trop dans r.

[Jean Dumoncel]

Bonjour,

x(p), y(p) et z(p) étant de simples scalaires, il n'y a pas besoin de les ajuster aux dimensions de X,Y et Z. Ainsi le code de magelan peut se simplifier en:

Code :

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

r1(p,:,:,:)=sqrt((X-x(p)).^2+(Y-y(p)).^2+(Z-z(p)).^2);

En ce qui concerne REPMAT, au delà de 2 dimensions, il faut fournir un vecteur:

Code :

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

repmat(x(p),[n n n]);

Oula! j'était pas bien réveillé moi ce matin!

Peux-tu nous montrer le message d'erreur complet ainsi que le code qui génère cette erreur?

[Sylvain_F]

Ok, comme ceci ça fonctionne bien :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
sp=zeros(n,n,n);pres=zeros(n,n,n);
for p=1:np   
    for ix=1:n
        for iy=1:n
            for iz=1:n
                pres(ix,iy,iz)=A(p)*exp(1i*(w(p)*t-k(p)*r(p,ix,iy,iz)))/(k(p)*r(p,ix,iy,iz));
            end
        end
    end
    sp=sp+real(pres);
end

Par contre comme celà ça ne marche pas :

Code :

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

1
2
3
4
sp=zeros(n,n,n)
for p=1:np
    sp=sp+real(A(p).*exp(1i.*(w(p)*t.*ones(n,n,n)-k(p).*r(p,:,:,:)))./(k(p).*r(p,:,:,:)));   
end

Je n'ai pas d'autre message d'erreur que ceci :

Code :

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

1
2
3
4
5
??? Error using ==> minus
Number of array dimensions must match for binary array op.
 
Error in ==> Source_3D at 36
           sp=sp+A(p).*exp(1i.*(w(p)*t.*ones(n,n,n)-k(p).*r(p,:,:,:)))./(k(p).*r(p,:,:,:));

[Invité]

Même remarque que précédemment pour le ones(n,n,n) restant. Maintenant tu utilises beaucoup de points inutilement: seul 1 reste indispensable ici.

Le problème est que ici size(r(p,:,:,:)) = [1 n n n] différent de [n n n]. Tu peux pour contrer cela utiliser la fonction SQUEEZE:

Code :

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

1
2
3
for p=1:np
    sp = sp + squeeze( A(p)*exp(1i*(w(p)*t-k(p)*r(p,:,:,:))) ./ (k(p)*r(p,:,:,:)) );   
end

[Sylvain_F]

Merci beaucoup !
Cela me fait gagner un temps précieux.