Discussion :

[uriotcea]

Bonjour,

J'ai un message d'erreur concernant cette directive dans l'example simplifi� ci-dessous.
Le message est "r declared 'threadprivate' after first use"
Ca vient du fait que ma classe a un constructeur mais ca je n'y peut rien.
Cet operator est appel� dans une boucle parall�lis�e et qui merde parce va variable stratic 'r' est commune aux threads, d'ou la tentaative d'utilisation de la directive "threadprivate"

Je ne sais pas comment faire autrement.

Merci de vos suggestions

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
 
 
class  my_class
{
 public :
  double x,y;
  my_class() { x=y=0; }
};
 
 
my_class&  my_class::operator*(const my_calss &p)
{
// déclaré en static parce que son constructeur prend trop de temps
  static my_class r; 
  #pragma omp threadprivate(r)
 
  r.x=r.y=0.;
  // differents calculs
 
   return(r);
}

[screetch]

avec un pointeur ?

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
 
class  my_class
{
 public :
  double x,y;
  my_class() { x=y=0; }
};
 
 
my_class&  my_class::operator*(const my_calss &p)
{
// déclaré en static parce que son constructeur prend trop de temps
  static my_class *r = new my_class(); 
  #pragma omp threadprivate(r)
 
  *r.x=*r.y=0.;
  // differents calculs
 
   return(*r);
}

il se peut que tu doives faire un if pour l'initialisation

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
 
class  my_class
{
 public :
  double x,y;
  my_class() { x=y=0; }
};
 
 
my_class&  my_class::operator*(const my_calss &p)
{
// déclaré en static parce que son constructeur prend trop de temps
  static my_class *r;
  #pragma omp threadprivate(r)
  if(!r) r = new my_class(); 
 
  *r.x=*r.y=0.;
  // differents calculs
 
   return(*r);
}

[uriotcea]

Oui �videmment,j 'aurais du y penser, merci beaucoup.

Une autre question au passage:
Quelle est la diff�rence entre omp_set_num_thread(0) & omp_set_num_thread(1), en d'autre termes si je veux iniber le mode //, quelle est la facon la plus �fficace ?

[screetch]

alors la aucune id�e, j'avais jamais utilis� openmp avant ta question hier ou avant-hier, donc je s�che :p