Question

如何将此代码与OpenMP并行化？我得到的结果是不正确的。

我尝试使用临时变量p1aux，p2aux和psumaux，因为Reduction子句不能与指针或内部函数一起使用。但正如我所说，结果不正确。

当我说“结果不正确”时，我会说：array1 + array2的计算存储在矩阵sumsse中。在组件sumsse[50] [50] [50]之前，计算是正确的，或多或少，但矩阵的其他组件始终为0。

你知道会发生什么吗？

如果有人可以帮助我，非常感谢你。

#include <stdio.h>
#ifdef __SSE2__
#include <emmintrin.h>
#endif

#include <omp.h>
#include <stdlib.h>
#include <math.h>

double __attribute__((aligned(16))) array1[256][256][256],array2[256][256][256];
double __attribute__((aligned(16))) sumsse[256][256][256];



int main ()
{

    int k,j,i,dim;
    dim=256;


    __m128d *p1= (__m128d*) &array1[0][0][0];
    __m128d *p2= (__m128d*) &array2[0][0][0];
    __m128d *psum= (__m128d*) &sumsse[0][0][0];

    __m128d *p1aux= (__m128d*) &array1[0][0][0];
    __m128d *p2aux= (__m128d*) &array2[0][0][0];
    __m128d *psumaux= (__m128d*) &sumsse[0][0][0];

    int nthreads =(8);
    omp_set_num_threads(nthreads);


    //initializa array2

    for(k = 0; k < dim; ++k){
        for(j = 0; j < dim; ++j){
            for(i = 0; i < dim; ++i){
                array2[k][j][i] = 1.0;
            }
        }
    }

    //initialize array1

    for(k = 0; k < dim; ++k){
        for(j = 0; j < dim; ++j){
            for(i = 0; i < dim; ++i){
                array1[k][j][i] = i + j + k;
            }
        }
    }

    //initialize array sumsse

    for(k = 0; k < dim; ++k){
        for(j = 0; j < dim; ++j){
            for(i = 0; i < dim; ++i){
                sumsse[k][j][i] = 0.0;
            }
        }
    }

    //add array1 + array2 with sse 

#pragma omp parallel firstprivate(p1,p2,p1aux,p2aux)
    {

    for(k = 0; k < dim; ++k){

        for(j = 0; j < dim; ++j){

#pragma omp for private(i)

                for(i = 0; i < dim; i += 2){

                    *psum = _mm_add_pd(*p1,*p2);

                    psum = psumaux + 1;
                    p1 = p1aux+1;
                    p2 = p2aux+1;
                    psumaux = psum;
                    p1aux = p1;
                    p2aux = p2;
                }
            }
        }

    }// end parallel


    //show some datas


    printf("Elementosse=10 sumsse=%f",sumsse[10][10][10]);
    printf("\n");
    printf("Elementosse=50 sumsse=%f",sumsse[50][50][50]);   
    printf("Elementosse=100 sumsse=%f",sumsse[100][100][100]);
    printf("\n");
    printf("Elementosse=120 sumsse=%f",sumsse[120][120][120]);
    printf("\n");
    printf("Elementosse=200 sumsse=%f",sumsse[200][200][200]);
    printf("\n");

    return 0;
}

Answer 1

而不是增加你的辅助变量，奇怪的是你应该以{{1}}，k，j和{{1}的方式为每个变量的索引设置一个封闭形式的东西。 }。然后你应该确保它们之间没有别名。

使用openmp将此c代码并行化的问题

1 个答案: