Discussion :

[S4sha]

Bonjour, je suis confronté au problème suivant.

J'ai une matrice 3-D A de taille (n,n,n), et je voudrais calculer la matrice B telle que :

pour tout i,j,k dans [2,n-1],

Code :

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

1
2
 
B(i,j,k) = f(A(i,j,k) - A(i+1,j,k)) + f(A(i,j,k) - A(i-1,j,k))  +  f(A(i,j+1,k) - A(i-1,j,k)) +  f(A(i,j,k) - A(i,j-1,k)) + f(A(i,j,k) - A(i,j,k+1)) + f(A(i,j,k) - A(i,j,k-1))

où f est une fonction de R dans R donnée. Si un point (i,j,k) est au bord, alors B(i,j,k) est la somme des différences entre A(i,j,k) et les valeurs aux points voisins de A, avec application de la fonction f (la même définition que celle donnée au dessus avec des termes qui sautent). Je voudrais faire ce calcul sans faire de boucle !

Merci du coup de main

merci !

[Jerome Briot]

Je voudrais faire ce calcul sans faire de boucle !

Pourquoi ?

[mr_samurai]

Salut,

Il faut que la fonction f soit vectorisé ( accepte un vecteur/matrice en entrée).

Sinon, le code devrait ressembler à ceci (sans boucle) :

Code :

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

1
2
ix = 2:n-1;
B = f(A(ix,ix,ix) - A(ix+1,ix,ix)) + ...

++

[S4sha]

Dut je ne veux pas faire de boucle car A est de taille astronomique.

Samourai ton calcul me plait bien. Sinon il y a la fonction "diff" qui est pas mal non plus

[Jerome Briot]

La fonction CONVN pourrait peut être être également utile ici ?

[S4sha]

Oui ça peut servir, merci à vous deux !

[mr_samurai]

Je doute que l'opération peut être fait par convolution.

Code :

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

1
2
3
%Avec DIFF le code deviendra très sympathique :D 
 
  B = f(-diff(A,1,1)) + ....

++

[S4sha]

Je préfère ta solution, car ça fait les différences à gauche et à droite en chaque points, alors qu'avec diff il n'y a qu'il seul sens je pense

[Jerome Briot]

Je doute que l'opération peut être fait par convolution.

Tu as raison, je n'avais pas remarqué que f s'appliquait aux termes deux à deux...