[Starpu-devel] Dépendances dynamiques

[Adrien Roussel <adrien.roussel@itwm.fraunhofer.de>]

Bonjour,

Je me permets de vous adresser ce mail car je suis actuellement en train de développer un petit programme d'algèbre linéaire où j'ai rencontré un problème d'_expression_ de dépendances via l'interface C de StarPU.
Mon problème est le suivant: j'aimerais programmer un produit scalaire entre deux vecteurs, qui prend la forme suivante.

double res = 0.0 ;
// Recuperation du nombre d'elements a calculer par partition
int local_size = x.size() / nb_partitions;
// boucle de calcul, où l'on cherche a assigner un indice de partition "i" a une tache
// On stocke le resultat dans le vecteur "tmp"
for(int i = 0; i < nb_partitions; ++i)
{
    int begin = i * local_size ;
    int end = local_size + begin ;

#pragma omp task depend(x[begin:end], y[begin:end]) depend(out: tmp[i])
    for(int k = begin; k < end; ++k)
      tmp[i] += x[k] * y[k]
}
#pragma omp taskwait
// Accumulation du vecteur tmp dans res
#pragma omp task depend(in: tmp[0], ..., tmp[nb_partitions - 1]) depend(out: res)
for(int i = 0; i < nb_partitions; ++i)
{
    res += tmp[i];
}

  Actuellement, pour palier à ce problème je suis obligé de mettre un point de synchronisation, modélisé par un "starpu_task_wait_for_all()" entre le calcul du produit scalaire, et l'accumulation des résultats "locaux" dans la variable résultat.
Malgré que cette solution marche et assure la cohérence du résultat, elle casse le parallélisme et peut provoquer des dommages sur la performance de la méthode.
Je ne peux pas créer de tâche de réduction prenant en dépendance d'entrée chacun de mes élèments stockés dans mon vecteur "tmp" car je ne connais sa taille qu'une fois l'exécution lancée. Il faudrait que je sois capable de rajouter des dépendances dynamiquement à ma tâche lors de l'exécution de mon application.

Avec l'interface C de StarPU, comment pourrais-je casser cette synchronisation et ainsi ne pas altérer le parallélisme de l'application ?

Merci beaucoup d'avance.

Cordialement,

-- 
|---------------------------------------------------------|
| Dr. Adrien Roussel
| Competence Center High Performance Computing
| Fraunhofer Institut ITWM
| Tel.   : +49 631 31600-4984
| E-Mail : adrien.roussel@itwm.fraunhofer.de
| Web    : https://adrien-roussel.eu
| Address: Fraunhofer-Platz 1
|	   D-67663 Kaiserslautern
|---------------------------------------------------------|