[F90] fonction where

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:

---

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

Bonjour,

Citation:

---

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:

---

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

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,

Re-,

Code:

---

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:

---

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

Bonjour xflr6,
c'est resolu

merci pour la solution