0

我想为 Intel Xeon Phi 协处理器（61 核）编写有效的并行应用程序，它执行五点模板计算。我写了两个版本的代码。

首先：我使用 OpenMP“#pragma omp parralel for”

void ParallelStencil(const double* macierzIn, double* macierzOut, const int m, const int n)
{
    int m_real = m + 2;
    int n_real = n + 2;

    TimeCPU t;
    t.start();
    #pragma omp parallel for schedule(static,1) shared(macierzIn, macierzOut)
    for(int i=1; i<m_real-1; ++i)
    {
        for(int j=1; j<n-1; ++j)
        {
            macierzOut[i * n_real + j] = Max(macierzIn[i * n_real + j], macierzIn[(i - 1) * n_real + j], macierzIn[(i + 1) * n_real + j],
                                             macierzIn[i * n_real + (j - 1)], macierzIn[i * n_real + (j + 1)]);
       }
    }
    t.stop();
    cout << "\nTime: " << t.time();
}

第二：我在 61 个核心之间划分矩阵。矩阵的每一部分由每个内核运行的 4 个硬件线程计算。在这个版本中，我尝试通过计算同一 L2 缓存周围的 4 个线程来减少缓存未命中。

void ParallelStencil(const double* macierzIn, double* macierzOut, int m, int n)
{
    int m_real = m + 2;
    int n_real = m + 2;
    int coreCount = threadsCount / 4;
    int tID, coreNum, start, stop, step;

    TimeCPU t;
    t.start();
    #pragma omp parallel shared(macierzIn, macierzOut, m, n, m_real, n_real, coreCount) private(tID, coreNum, start, stop, step)
    {
        tID = omp_get_thread_num();
        coreNum = tID / 4;
        start = tID % 4 + ((m / coreCount) * coreNum) + 1;
        stop = (m / coreCount) * (coreNum + 1) + 1;
        if(coreNum == coreCount - 1 && stop != m_real - 1)
        {
                stop = m_real -1;
        }
        step = 4;

        for(int i=start; i<stop; i+=step)
        {
            for(int j=1; j<n+1; ++j)
            {
                macierzOut[i * n_real + j] = Max(macierzIn[i * n_real + j], macierzIn[(i - 1) * n_real + j], macierzIn[(i + 1) * n_real + j],
                                                 macierzIn[i * n_real + (j - 1)], macierzIn[i * n_real + (j + 1)]);

            }
        }
    }
    t.stop();
    cout << "\nTime: " << t.time();
}

在这个 wersion 循环中，矩阵每个部分的迭代都以这种方式执行：
i=0 -> 线程 0
i=1 -> 线程 1
i=2 -> 线程 2
i=3 -> 线程 3
i=4 -> 线程0
...

运行此代码后。第二个版本比较慢。但为什么？

4

1 回答 1

0

这可能更多的是评论而不是答案。在深入研究有效利用缓存的问题之前，您应该修复这两个代码，以使它们等效。目前他们不是。

差异#1

第一个代码：

int m_real = m + 2;
int n_real = n + 2;

第二个代码：

int m_real = m + 2;
int n_real = m + 2; // <---- m instead of n

差异#2

第一个代码：

for(int j=1; j<n-1; ++j)

第二个代码：

for(int j=1; j<n+1; ++j) // <---- n+1 instead of n-1

如果您的矩阵碰巧不是平方和m > n，那么第二个代码肯定会更慢，因为它必须计算更多。

于 2014-09-04T15:49:55.653 回答