OpenMP 中的陣列縮減是否可行？

OpenMP 中的陣列縮減

問題：

問題：

問題：

並縮減程式需要數組數但它在OpenMP 中被認為是不可能的。有替代方案嗎？

答案：

• 是的，OpenMP 中可以進行陣列縮減。這裡有兩種替代方法：
• 方法 1：使用關鍵部分進行私有求和

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];
        }
    }
}

在每個執行緒中建立陣列的私有版本。 並行執行縮減。

合併關鍵部分中的值以保留正確的值
• 方法2：不使用臨界區的私有求和

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;

建立一個維度為[10 * nthreads]的私有陣列。 並行執行歸約並將結果儲存在私有中數組。 將數值合併到原始陣列中，不包含臨界區。

以上是OpenMP 中的陣列縮減是否可行？的詳細內容。更多資訊請關注PHP中文網其他相關文章！