Question

我正在尝试并行化一段代码，但是存在一些问题，在共享映射中插入的关键部分会在我的执行中造成瓶颈。有一种方法可以并行插入吗？

第二个问题涉及如何在地图上使用openMP迭代预定义次数。这是代码，我在Clion上使用g ++ - 5.

else if (PERF_ROWS <= MAX_ROWS && function_switch == false)
{
    int array_dist_perf[PERF_ROWS];
    int array_dist[MAX_ROWS];

    multimap<float, int> temp_map;

    #pragma omp parallel for schedule(dynamic) private(array_dist_perf) shared(temp_map)
    for (int i = 0; i < MAX_COLUMNS; i++)
    {
        if (i % PERF_CLMN == 1) continue;

        for (int j = 0; j < PERF_ROWS; j++) //truncation perforation
        {
            array_dist_perf[j] = abs(input[j] - input_matrix[j][i]);
        }

        float av = mean(PERF_ROWS, array_dist_perf);

        float temp_score = score_func(av);

        #pragma omp critical(map_update)
        temp_map.insert({temp_score, i});
    }

    map<float,int>::reverse_iterator rit;

    int iter = 0;

    #pragma omp parallel for schedule(dynamic) private(array_dist) shared(iter)
    for (rit=temp_map.rbegin(); rit!=temp_map.rend(); rit++)
    {
        int s = rit->second;

        for (int k = 0; k < MAX_ROWS; k++)
        {
            array_dist[k] = abs(input[k] - input_matrix[k][s]);
        }

        float av_real = mean(MAX_ROWS, array_dist);

        float score_real = score_func(av_real);

        rank_function(score_real, s);

        #pragma omp atomic
        iter++;

        if (iter == THRESHOLD_NUM)

            break;
    }
}

Answer 1

在Open MP 4（应该由g ++ 5支持）中，可以定义自己的缩减。我不想做一些不必要的重复，所以请看看这个答案：C++ OpenMP Parallel For Loop - Alternatives to std::vector

如果我理解你的问题是正确的，你可以使用order and schedule（static，chuncksize），其中chunksize是一个线程获得的迭代次数。

Maps和openMP并行调度

1 个答案: