Discussion :

[le_fou1]

Bonjour,

j'essaie de trouver une solution avec la fonction intrinseque "where"

Ce que j'aimerais obtenir, c'est le nombre des elements d'un vecteur (array) qui satisfait une certaine conditon logique et puis les indices des elements qui satisfont une expression logique :

Code :

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

1
2
3
4
5
6
7
8
9
do jpm=1, ny_c
    count_j_f=count( abs(y_c(jpm) - y_f) <= y_res_c / 1.99)  ! y_c et y_f sont deux arrays de 140 elements
    where( abs(y_c(jpm) - y_f) <= y_res_c / 1.99)
        B=1.0
    elsewhere
        B=0.0
    END WHERE
end do 
result_count=maxloc(ARRAY=B)

est ce quelqu'un peut m'expliquer pourquoi ca ne fonctionne pas et me suggerer une solution?

merci

[xflr6]

Bonjour,

les indices des elements qui satisfont une expression logique

Je ne sais pas si c'est votre problème, mais maxloc ne va retourner que l'indice du premier maximum rencontré.

Il y a peut-être aussi une erreur d'implémentation

Code :

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

1
2
3
do jpm=1, ny_c
    count_j_f=count( abs(y_c(jpm) - y_f) <= y_res_c / 1.99)  ! y_c et y_f sont deux arrays de 140 elements
    where( abs(y_c(jpm) - y_f) <= y_res_c / 1.99)

Les count, where etc, sont en général utilisés pour éviter les boucles, alors le parcours jpm et le masque à base de y_c(jpm) - y_f(:) peut être une piste d'erreur d'implémentation (peut aussi bien être nécessaire à l'algo certes).

Cordialement,
xflr6

[le_fou1]

Bonjour xflr6,

Je te donne un exemple de ce qui se passe avec la fonction where.

Si on suppose que 'abs(y_c(jpm) - y_f)' me donne les valeurs suivantes : (/0.75, 0.23, 0.1, 1, 2, 3, 4...), ces dernieres valeurs doivent satisfaire la condition logique tels que :


(/0.75, 0.23, 0.1, 1, 2, 3, 4...) <=0.5 ! sachant que 0.5 n'est autre que y_res_c / 1.99

donc si j'utilise where, la reponse doit etre la suivante :

le nombre d'elements est egal deux (car les elements qui satisfont la condition sont 0.23 et 0.1) et donc les indices doivent etre 2 et 3 respectivement. Mais en appliquant le code que j'ai utilise, j'arrive a detecter deux elements par contre je n'arrive pas a obtenir les bons indices, car ce que la fonction where fait dans mon code, elle me met les 1 a la fin du vecteur B et les valeurs zeros au debut du vecteur B de la facon suivante :

B=(/0, 0, 0, 0, 0,........1,1)

comme tu peux remarquer, j'ai obtenu deux valeurs egal 1 ds B, donc je vois que la fonction 'where' a reussi a detecter deux elements qui satisfont ma condition, par contre je n'obtiens pas le bon indice, et c'est ca ma question, comment trouver les bons indices ?

merci d'avance,

[xflr6]

Re-,

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 program test
   real, dimension(7) :: vmask, b
   vmask=(/0.75, 0.23, 0.1, 1., 2., 3., 4./)
   where ( vmask <=0.5 )
     b=1.
   elsewhere
     b=0.
   endwhere
   do i=1,7
     print *, i, vmask(i), b(i)
   enddo
 end program test

donne

Code :

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

1
2
3
4
5
6
7
           1  0.750000000       0.00000000    
           2  0.230000004       1.00000000    
           3  0.100000001       1.00000000    
           4   1.00000000       0.00000000    
           5   2.00000000       0.00000000    
           6   3.00000000       0.00000000    
           7   4.00000000       0.00000000

Le problème pour moi semble donc ailleurs que dans le fonctionnement du where, non ?

Edit: en d'autres termes, ça vient probablement de l'implémentation. Dans l'algo que vous avez écrit, vous redéfinissez b ny_c fois, la dernière fois avec y_c(ny_c) dans le masque.

Edit 2: j'espère que si les portées des tableaux utilisés dans le test du where et dans les assignements du where sont différentes, le compilateur crache. C'est quasi sûr dans les tableaux statiques, mais dans les dynamiques, ça passe probablement au travers de la vérif. Peut-être chercher par là pourrait s'avérer payant aussi.

Cordialement,
xflr6

[le_fou1]

Bonjour xflr6,
c'est resolu

merci pour la solution