logo标签列表

OpenMP嵌套并行化

c++cmultiprocessingopenmp
5

我有一个带有两个内部并行区域的外部并行区域。是否可以将2个线程放入外部并行区域，将每个内部并行区域中的4个线程放入其中？我尝试了类似于这样的东西，但它似乎不能按照我想要的方式工作。有什么建议吗？

start_r = omp_get_wtime();
omp_set_nested(1);
omp_set_num_threads(2);
#pragma omp parallel 
{
    printf("Thread %d executes the outer parallel region\n",omp_get_thread_num());
    omp_set_num_threads(4);
    #pragma omp parellel for private(i,j,color)schedule(guided, chunk) default(shared) 
    {

// Blur
    for (int i = 1; i < x-1; i++)
        for (int j = 1; j < y-1; j++)
            for (int k = 0; k < 3; k++)
            {   
                wynik = 0;
                wynik = ((color[(i-1)][((j - 1))][k] +
                        color[(i-1)][j][k] +
                        color[(i-1)][(j + 1)][k] +
                        color[i][(j - 1)][k] +
                        color[i][j][k] +
                        color[i][(j + 1)][k] +
                        color[(i+1)][(j - 1)][k] +
                        color[(i+1)][j][k] +
                        color[(i+1)][(j + 1)][k])/9);
                if (wynik>255)wynik = 255;
                if (wynik<0)wynik = 0;
                color2[i][j][k] =  wynik;
            }
            stop_r = omp_get_wtime();
            cout << "Wyostrzenie zejelo : " << (stop_r-start_r) <<" sekund"<< endl;
            cout<<omp_get_nested( )<<endl;
            cout<<"Ilość wątków dla rozmycia : "<<omp_get_num_threads( )<<endl;
            printf("Thread %d executes the inner parallel region\n",omp_get_thread_num());
        }
        omp_set_num_threads(4);
#pragma omp parellel for schedule(guided, chunk) privat(i,j,color) default(shared)
{
// Sharp
    for (int i = 1; i < x - 1; i++)
        for (int j = 1; j < y - 1; j++)
            for (int k = 0; k < 3; k++)
            {
                wynik = 0;
                wynik = (color[(i-1)][(j - 1)][k] * (0) +
                        color[(i-1)][j][k] * (-1) +
                        color[(i-1)][(j + 1)][k] * (0) +
                        color[i][(j - 1)][k] * (-1) +
                        color[i][j][k] * 20 +
                        color[i][(j + 1)][k] * (-1) +
                        color[(i+1)][(j - 1)][k] * (0) +
                        color[(i+1)][j][k] * (-1) +
                        color[(i+1)][(j + 1)][k] * (0))/16;
                wynik = wynik % 255;
                color3[i][j][k] = wynik;



            }
            cout<<omp_get_nested( )<<endl;
            cout<<"Ilość wątków dla wyostrzenia : "<<omp_get_num_threads( )<<endl;
            printf("Thread %d executes the inner parallel region\n",omp_get_thread_num());
            }
        }
    for (int j = 0; j < y; j++)
        for (int i = 0; i < x; i++)     
            {
                fwrite(color2[i][j], 1, 3, fp2);
                fwrite(color3[i][j], 1, 3, fp3);

            }



    fclose(fp);
    fclose(fp2);
    fclose(fp3);

    system("PAUSE");
    return 0;
}
}
2
我建议您也使用OpenMP标记。 - terence hill
2个回答
3

这适用于VS2012。

示例：

#include <iostream>
#include <omp.h>

int main()
{
    omp_set_nested(2);

    #pragma omp parallel num_threads( 2 )
    {
        int threadID1 = omp_get_thread_num();

        #pragma omp parallel num_threads( 4 )
        {
            int threadID2 = omp_get_thread_num();

            #pragma omp critical
            {
                std::cout << "tID1: " << threadID1 << std::endl;
                std::cout << "tID2: " << threadID2 << std::endl;
                std::cout << std::endl;
            }
        }
    }

    return EXIT_SUCCESS;
}

输出：

tID1: 0
tID2: 0

tID1: 0
tID2: 2

tID1: 0
tID2: 1

tID1: 0
tID2: 3

tID1: 1
tID2: 0

tID1: 1
tID2: 1

tID1: 1
tID2: 2

tID1: 1
tID2: 3
1
是的，我尝试将线程数放在内部区域，但它们仍然只有两个线程作为外部区域。 - Lobo
这应该可以设置在for(...)上工作的线程数量，而不是在一个区域上工作。 - terence hill
即使我将4个线程放入并行循环中，omp_get_num_threads（）-仍然给出相同的答案（该循环中只有两个线程在工作）。看起来这取决于我在外部区域之前设置了多少个线程，但即使我设置了8个线程，它们都在内部区域中工作（我无法将它们分成同时在两个区域中工作）。 - Lobo
很难回答，更新你的问题并发布可工作的代码，这样我才能测试它。 - terence hill
网页内容由stack overflow 提供, 点击上面的
可以查看英文原文，
原文链接