Comment effectuer une réduction de tableau dans OpenMP ?

Réduction des tableaux dans OpenMP

Question :

Comment l'effectuer dans OpenMP pour les tableaux Opération de réduction ?

Réponse :

OpenMP prend en charge la réduction de tableau. Il y a deux manières de procéder :

Méthode 1 : Utiliser des variables "privées"

pour créer une copie privée de S pour chaque thread, les remplir en parallèle, puis Fusionnez-les en S dans la zone :

int A[] = {84, 30, 95, 94, 36, 73, 52, 23, 2, 13};
int S[10] = {0};
#pragma omp parallel
{
    int S_private[10] = {0};
    #pragma omp for
    for (int n = 0; n < 10; ++n)
    {
        for (int m = 0; m <= n; ++m)
        {
            S_private[n] += A[m];
        }
    }
    #pragma omp critical
    {
        for (int n = 0; n < 10; ++n)
        {
            S[n] += S_private[n];
        }
    }
}

méthode 2 : Créez un tableau de dimensions 10*nthreads à l'aide d'un tableau multidimensionnel

, remplissez-le en parallèle, puis fusionnez-le dans S sans section critique :

int A[] = {84, 30, 95, 94, 36, 73, 52, 23, 2, 13};
int S[10] = {0};
int *S_private;
#pragma omp parallel
{
    const int nthreads = omp_get_num_threads();
    const int ithread = omp_get_thread_num();

    #pragma omp single
    {
        S_private = new int[10 * nthreads];
        for (int i = 0; i < (10 * nthreads); i++)
            S_private[i] = 0;
    }
    #pragma omp for
    for (int n = 0; n < 10; ++n)
    {
        for (int m = 0; m <= n; ++m)
        {
            S_private[ithread * 10 + n] += A[m];
        }
    }
    #pragma omp for
    for (int i = 0; i < 10; i++)
    {
        for (int t = 0; t < nthreads; t++)
        {
            S[i] += S_private[10 * t + i];
        }
    }
}
delete[] S_private;

Ce qui précède est le contenu détaillé de. pour plus d'informations, suivez d'autres articles connexes sur le site Web de PHP en chinois!