OpenMP reduction en C++

Version imprimable

04/07/2011, 15h55
uriotcea

OpenMP reduction en C++
Bonjour,

Dans un code C++, j'essaye d'utiliser la direction d'openMP 3.0 REDUCTION sur une variable de type tableau ou de type pointeur et cela me génére un erreur à la compilation:

"'bb' has invalid type for 'reduction'"
Code:

1 2 3 4 5 6 7 N=100; double bb[100]; #pragma omp parallel for reduction(+:bb) for (int i=0;i<N;i++) { bb[i]++; }
Normallement depuis openMP 2.0 on peut faire des reduction sur tableau, ou alors il y a quelquechose qui m'echape.
04/07/2011, 16h13
Joel F

bb[i]++ n'est pas une reduction.

??

Si je fais ce qui suit sur 4 processeurs j'obtient bien 100

Code:

1
2
3
4
5
6
7
8
9
 
  N=100;
  double bb;
  #pragma omp parallel for reduction(+:bb) 
  for (int i=0;i<N;i++) {
     bb++;
  }
 
  printf("%d\n",bb)

04/07/2011, 18h29
Joel F

bb++ deja c'est pas bb[i]++

Une reduction comme son nom l'indique passe de N dimension a N-1 ce qui n'est pas le cas de ton code.
04/07/2011, 18h58
uriotcea

Passer de bb++ à bb[i]++ était l'objet de ma question!
Es-tu sure de bien saisir la directive "Reduction" de openMP. Son role est de paralléliser par reduction une operation comme l'addition par example. C'est bien ce que je fais ici.
04/07/2011, 19h36
Matthieu Brucher

Réduction = par exemple somme des bb, pas bb[i]++ !
04/07/2011, 20h27
Joel F

Citation:

Envoyé par uriotcea

Passer de bb++ à bb[i]++ était l'objet de ma question!
Es-tu sure de bien saisir la directive "Reduction" de openMP. Son role est de paralléliser par reduction une operation comme l'addition par example. C'est bien ce que je fais ici.

Non, la reduction serait de faire la somme de bb[i]. la tu incrementes chaque bb[i] separement.
04/07/2011, 20h35
Aleph69

Pour démarrer en openMP :
http://www.idris.fr/data/cours/parallel/openmp/
04/07/2011, 23h27
uriotcea
Extrait du cours ci-dessous au chapitre "reduction"
Code:

1 2 3 4 5 6 7 8 9 10 !$OMP PARALLEL !$OMP DO REDUCTION(+:s) REDUCTION(*:p,r) do i = 1, n s = s + 1 p = p * 2 r = r * 3 end do !$OMP END DO !$OMP END PARALLEL
A fin, c'est certain s=N; Je ne sais pas pourquoi ca semble choquer tout le monde, je fais ca depuis des années et ca fonctionne trés bien !!!!
04/07/2011, 23h55
Joel F

relis avec le doigt. Si tu fait

for(i=0;i<100;i++) bb[i]++;

tu fais un +1 sur chaque case du tableau, y a pas de reduction.
Lire c'est bien, comprendre c'est mieux.
05/07/2011, 00h28
uriotcea

Dit tu te la joue là!
Evidement que je veux faire +1 à chaque case du tableau.
Lit bien ma question initiale avec ton doigt aussi si c'est nécessaire.
05/07/2011, 03h46
Joel F

non.

Si tu veux faire a[i]++, ce n'est pas une reduction, c'ets une boucle parallele normale.

Bon mon example est peu etre trop minimaliste alors voici ce que je veux faire et qui ne fonctionne pas

Code:

1
2
3
4
5
6
7
8
 
#pragma omp parallel for reduction(+:tab,tabp) private(p,power)
  for (int i=0;i<100000;i++) {
    power=cord[i]*kpower;
    p=(unsigned int)(k*(cord[i]-min));
    tab[p]++;
    tabp[p]+=power;
  }

Et voici ce que je suis obligé de faire pour que ca fonctionne

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
   int pp;
#pragma omp parallel for private(p,pp,power)
  for (int i=0;i<100000;i++) {
    pp=omp_get_thread_num();
 
    power=cord[i]*kpower;
    p=(unsigned int)(k*(cord[i]-min));
    new_tab[pp][p]++;
    new_tabp[pp][p]+=power;
  }
 
  // J'effectue la reduction moi meme
  for (int j=0;j<omp_get_max_threads();j++) {
    tab[i]+=new_tab[j][i];
    tabp[i]+=new_tabp[j][i];
  }

Citation:

Envoyé par uriotcea

Bon mon example est peu etre trop minimaliste alors voici ce que je veux faire et qui ne fonctionne pas

Code:

1
2
3
4
5
6
7
8
 
#pragma omp parallel for reduction(+:tab,tabp) private(p,power)
  for (int i=0;i<100000;i++) {
    power=cord[i]*kpower;
    p=(unsigned int)(k*(cord[i]-min));
    tab[p]++;
    tabp[p]+=power;
  }

Et voici ce que je suis obligé de faire pour que ca fonctionne

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
   int pp;
#pragma omp parallel for private(p,pp,power)
  for (int i=0;i<100000;i++) {
    pp=omp_get_thread_num();
 
    power=cord[i]*kpower;
    p=(unsigned int)(k*(cord[i]-min));
    new_tab[pp][p]++;
    new_tabp[pp][p]+=power;
  }
 
  // J'effectue la reduction moi meme
  for (int j=0;j<omp_get_max_threads();j++) {
    for (int i=0;i<step;i++) {
      tab[i]+=new_tab[j][i];
      tabp[i]+=new_tabp[j][i];
    }
  }

05/07/2011, 09h02
Matthieu Brucher

Citation:

Envoyé par uriotcea

Bon mon example est peu etre trop minimaliste alors voici ce que je veux faire et qui ne fonctionne pas

Ca n'a strictement rien à voir avec ce que tu nous a montré...
Le problème, c'est que tu fais ta réduction sur tab et non sur tab[p], donc utilises une variable intermédiaire, ça ira mieux.
05/07/2011, 09h38
uriotcea

Ce que j'ai intialement montré, générait une erreur de compilation que je ne savais et ne sais toujours pas resoudre. c'était le but de l'example. La discution a un peu divergé sur le font plutot que sur la forme.

maitenant je ne vois pas trop ce que tu entends par passer par une variable intermediére.
05/07/2011, 11h00
Matthieu Brucher

Tu fais une réduction sur tab[p], mais je ne sais pas si c'est exprimable dans OpenMP. Donc crées une variable temporaire tabp sur laquelle tu fais ta réduction, et tu feras à la sortie de ta boucle OpenMP un tab[p] = tabp.
05/07/2011, 11h02
uriotcea

Oui, c'est ce que j'ai fini par faire.
J'ai aussi fini par comprendre que dans openMP 3.0 la reduction sur tableau etait possible en fortran mais pas en C. Cette possibilité semble arriver avec la derniére mouture, la 3.1 je crois.
Merci de vos remarques et commentaires à tous.

Citation:

Envoyé par uriotcea

J'ai aussi fini par comprendre que dans openMP 3.0 la reduction sur tableau etait possible en fortran mais pas en C.

l'exemple fortran que tu donne n'utilise pas de tableau :koi:
Normalement la réduction c'est uniquement pour des type primitif et une opération simple.
Aprés tu peut utiliser une section critique pour faire la reduction. Un truc genre
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 #pragma omp parallel { //definition des variables propre à un thread #pragma omp for for (int i=0;i<100000;i++) { //code de la boucle } pragma omp critical { //code de reduction } }

L'example que j'ai donnée en C sur un tableau qui ne compile pas fonctionne bien en fortran ou la reduction est valide sur le tableau. C'est ce que j'observe et c'est également indiqué dans les specif openMP 3.0 (voir ci-dessous

pour le C
Code:

1 2 3 Aggregate types (including arrays), pointer types and reference types may not appear in a reduction clause
Pour le fortran
Code:

1 2 3 4 An original list item with the ALLOCATABLE attribute must be in the allocated state at entry to the construct containing the reduction clause. Additionally, the list item must not be deallocated and/or allocated within the region.