Discussion :

[Ulath0]

Bonjour,

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

Code :

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

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 !

[Hydro999]

A mon avis ce qui bloque c'est le if

Code :

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

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 :

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

if (A[,j]>0)

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

Code :

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

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 :

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

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

Cdt