Parallélisation OpenMP en C

Version imprimable

Bonjour à tous,

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

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:

#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

10/12/2013, 16h49
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.
10/12/2013, 17h07
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 ?
11/12/2013, 08h23
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.

Citation:

Envoyé par pité60

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.