Accélérer une boucle imbriquée

Bonjour,

Je me retrouve à devoir écrire cette boucle très inefficace :

Code:

---

1 2 3 4 5 6

for (i in 1:n) { for (j in 1:n) { if (A[i,j]>0) W_nc[nums[i],nums[j]] = W[i,j] } }

---

où "nums" est un vecteur de taille n, contenant des indices dans [1,nc] correspondant à i et j dans [1,n]. Les autres éléments de W_nc valent 0.

Y-a-t'il un moyen d'écrire quelque chose comme "for ( (i,j) in all_positive_pairs(A) ) do W_nc[nums[i],nums[j]] = W[i,j]" ? Ou faut-il le coder en C ? (Auquel cas il faudra que j'apprenne à le faire, je code en C/C++ et un peu en R, mais n'ai pas encore mélangé les deux :-) )

Merci d'avance !

A mon avis ce qui bloque c'est le if

Code:

---

if (A[i,j]>0)

---

Tu parcours A élément après élément et il y en a n^2.
Il doit être possible de gagner du temps en testant

Code:

---

if (A[,j]>0)

---

Cela revient à tester d'un coups tout les éléments du vecteur.

Code:

---

1 2 3 4 5 6 7 8

M=matrix(1:9,nrow=3) [,1] [,2] [,3] [1,] 1 4 7 [2,] 2 5 8 [3,] 3 6 9   M[,1]>4 [1] FALSE FALSE FALSE

---

et pour les plus pressés

Code:

---

1 2 3 4 5

> M>4 [,1] [,2] [,3] [1,] FALSE FALSE TRUE [2,] FALSE TRUE TRUE [3,] FALSE TRUE TRUE

---

Cdt