Discussion :

[cpalperou]

Bonjour,
mon problème est résumé dans le titre de c epost.
J'ai trouvé une méthode mais y a t'il mieux?
Soit K ma matrice initiale de dimension (36*36).

Code :

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

1
2
3
4
5
 
Kdiag = diak(K);
Kdiag = diag(Kdiag);
Knotd = 2.*(K-Kdiag);
K = Kdiag + Knotd;

merci

[Jean Dumoncel]

Bonjour,

mieux en terme de quoi? rapidité? occupation mémoire?

Ta méthode peut se simplifier :

Code :

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

Kp = 2*K - diag(diag(K));

[Jerome Briot]

Une autre approche :

Code :

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

1
2
3
4
N = 36;
tmp = K(1:N+1:N*N);
K = 2*K;
K(1:N+1:N*N) = tmp;

[cpalperou]

merci,
je voulais juste une formulation plus élégante que la mienne (qui est assez bourrin!)
Merci bein à vous 2