1

如何使用 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;
}
4

1 回答 1

1

而不是奇怪地增加您的辅助变量，您应该根据k,j和. 那么你应该确保它们之间没有混叠。imax

于 2012-03-12T17:11:56.373 回答