Discussion :

[iznogoud_23]

Bonjour à tous

Voila ma question du jour
Comment, en fortran 90, récupérer la position (les indices) des valeurs d'un tableau qui vérifient un test logique.

Avec "WHERE", il est possible d'effectuer une action sur ces valeurs, mais pas de récupérer leur position.
Bref, comment faire?

Merci d'avance

[Ladgalen]

Bonjour

Tu ne peux pas les enregistrer dans un tableau au passage ?

Tu fais ta boucle where qui fait un test sur les valeurs de ton tableau et tu dis que si le test logique est vrai tu enregistres l'indice courant dans un tableau et tu incrémente l'indice du tableau en question.

[iznogoud_23]

Bonjour

Je veux bien, mais justement, comment avoir accès à l'indice courant?

C'est possible avec une structure en "DO ENDDO" avec un "IF" imbriqué, mais je voudrais éviter cette solution car et est beaucoup trop lente. Je cherche justement a optimiser un programme qui est très lent car comportant de nombreuses boucles et test imbriqué.
Le but était justement de récupérer ces indices afin de limiter l'étendue des boucle "DO ENDDO"

[Ladgalen]

Exact, je comprend mieux ton problème. Je pense malheureusement que le bloc WHERE n'est pas fait pour. Il est indiqué que les instructions dans un bloc where doivent concerner uniquement des affectations dans un tableau.

En fait l'instruction WHERE permet simplement de rajouter des filtres sur une instructions du type :

Code :

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

1
2
3
4
 
sqrta(:) = sqrt( a(:) )
 
where ( a > 0 ) sqrta(:) = sqrt( a(:) )

D'autre part (sauf si tu fais de la parallélisation) un bloc do avec un if ne doit pas avoir une efficacité très différentes de celle d'un bloc where. Peut être que where est bien adapté à l'utilisation d'instruction OpenMP mais je ne m'y connais pas trop. Pour te convaincre voici un petit test :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 
# user@dvp  ~
> ifort -O0 -o exe test.f90
 
# user@dvp  ~
> exe
bloc where
temps cpu      :    0.0330
temps systeme  :    0.0327
 
bloc do if then else end if end do
temps cpu      :    0.0120
temps systeme  :    0.0119
 
# user@dvp  ~
> gfortran -o exe test.f90
 
# user@dvp  ~
> exe
bloc where
temps cpu      :    0.0150
temps systeme  :    0.0150
 
bloc do if then else end if end do
temps cpu      :    0.0130
temps systeme  :    0.0130

Et voilà le code en question :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
 
program test
 
  implicit none
 
  integer            :: i
  integer, parameter :: n = 1000000
  integer            :: t2, t1, freq, maxt, nb_periodes
  real, dimension(n) :: a, sqrta
  real               :: x, cpu1, cpu2, tps_systeme
 
  call system_clock( count_rate=freq, count_max=maxt )
 
  ! a entre -1 et 1
  do i = 1, n
    call random_number(x)
    a(i) = 2.d0 * x - 1.d0
  end do
 
  !* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
  ! bloc where
  !* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
 
  call system_clock( t1 )
  call cpu_time( cpu1 )
 
  where ( a(:) > 0.d0 ) 
    sqrta(:) = sqrt(a(:))
  else where
    sqrta(:) = 0.d0
  end where
 
  call cpu_time( cpu2 )
  call system_clock( t2 )
 
  nb_periodes = t2 - t1
  do while ( nb_periodes < 0 ) 
    nb_periodes = nb_periodes + maxt
  end do
  tps_systeme = real(nb_periodes) / real(freq)
 
  write(*,"('bloc where')")
  write(*,"('temps cpu      :',F10.4)") cpu2 - cpu1
  write(*,"('temps systeme  :',F10.4)") tps_systeme
 
  !* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
  ! bloc do if else end do
  !* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
 
  call system_clock( t1 )
  call cpu_time( cpu1 )
 
  do i = 1, n
    if( a(i) > 0.d0 ) then
      sqrta(i) = sqrt(a(i))
    else
      sqrta(i) = 0.d0
    end if
  end do
 
  call cpu_time( cpu2 )
  call system_clock( t2 )
 
  write(*,*)
  write(*,"('bloc do if then else end if end do')")
 
  nb_periodes = t2 - t1
  do while ( nb_periodes < 0 ) 
    nb_periodes = nb_periodes + maxt
  end do
  tps_systeme = real(nb_periodes) / real(freq)
 
  write(*,"('temps cpu      :',F10.4)") cpu2 - cpu1
  write(*,"('temps systeme  :',F10.4)") tps_systeme
 
end program test

[iznogoud_23]

Oui c'est bien ce que je pensais.
Merci pour le test, finalement, il semble que where soit plus lent que if, je suis étonné.

Quoi qu'il en soit, j'ai trouvé une solution, a base de where/minloc et count. C'est sioux, mais ca marche. Maintenant reste a voir si j'y gagne en temps

[Ladgalen]

Oui au final c'est le gain de temps le plus important !

Il faut aussi voir que si c'est un script à usage unique le facteur limitant peut être l'écriture du script et il vaut mieux éviter les trucs trop sophistiqué !