Проблема распараллеливания этого кода c с openmp - PullRequest
Новая
       8

Проблема распараллеливания этого кода c с openmp

1 голос
/ 12 марта 2012

Как я могу распараллелить этот код с 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;
}

1 Ответ

1 голос
/ 12 марта 2012

вместо того, чтобы причудливо увеличивать вашу вспомогательную переменную, вы должны иметь закрытую форму для индекса каждой из них в терминах k, j, i и max. Тогда вам следует убедиться, что между ними нет псевдонимов.

...