Discussion :

[pité60]

Bonjour à tous,

J'ai une routine en C à paralléliser en utilisant OpenMP :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
void spmv(int n, int *rowptr, int *colind, double *val, double alpha, double *x, double beta, double *y){
 
  int i, j;
 
  for(i=0; i<n; i++){
    y[i] = beta*y[i];
    for(j=rowptr[i]; j<rowptr[i+1]; j++){
      y[i] += alpha*val[j]*x[colind[j]];
    }
  }
 
  return;
}

J'ai déjà essayé de placer dans la seconde boucle un

Code :

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

#pragma omp parallel for reduction(+:y[i])

mais le langage n'accepte pas le type de y[i].
Au final, il n'y a pas d'interdépendances entre les y[i] dans la boucle principale, serait-il possible d'utiliser une parallélisation classique sur celle ci ?

Merci d'avance

[Ehonn]

Bonjour

Il suffit d'utiliser une variable intermédiaire s et d'écrire s dans y[i] à la fin de la boucle.

Tu peux, en effet, paralléliser la boucle externe.

[pité60]

Bonjour,

J'avais utilisé cette solution, seulement il me semble que le temps de calcul avait beaucoup augmenté.
Il n'y a pas moyen de paralléliser sans modifier le code, juste en ajoutant du code OpenMP ?

[Ehonn]

Ça peut dépendre du nombre de données à traiter; du nombre d'accès mémoire par rapport au nombre de calcul.

Il faudrait que tu nous donnes un exemple qui tourne et les temps que tu obtiens.

Il n'y a pas moyen de paralléliser sans modifier le code, juste en ajoutant du code OpenMP ?

?
En parallélisant la boucle externe.
Pour la boucle interne, passer par une variable intermédiaire n'est pas une grosse modification.