1

openMP コードの何が問題になっていますか? 常に 1 つのスレッドしか必要とせず、非並列バージョンと同じ時間で動作します

template <typename T>
Matrix<T>* Matrix<T>::OMPMultiplication(Matrix<T>* A, Matrix<T>* B){ 

    if(A->ySize != B->xSize)
      throw;

    Matrix<T>* C = new Matrix<T>(A->xSize, B->ySize);

    sizeType i, j, k;
    T element;

    #pragma omp parallel for private(i, j)
    {
      #pragma omp for private(i, j)
      for( i = 0; i < A->xSize; i++ )
          cout<<"There are "<<omp_get_num_threads()<<" threads"<<endl;

          for(j = 0; j < B->ySize; j++){

              C->matrix[i][j] = 0;
              for(k = 0; k < A->ySize; k++){
                  C->matrix[i][j] += A->matrix[i][k] * B->matrix[k][j]; 
              }   

      }   
    }   
    return C;
}
4

2 に答える 2

2

まず第一に、ループの一部が欠落{}しており、iループの反復ごとに変数kをプライベートにする必要がありiます。parallelただし、プラグマとforプラグマの組み合わせ方も混同していると思います。parallelfor ループを正常に並列化するには、それをプラグマ内に配置してから、プラグマ内に配置する必要がありますfor。これを行うには、コードを次のように変更できます

#pragma omp parallel private(i, j, k)
{
    #pragma omp for
    for( i = 0; i < A->xSize; i++ ) {
        cout<<"There are "<<omp_get_num_threads()<<" threads"<<endl;

        for(j = 0; j < B->ySize; j++) {

            C->matrix[i][j] = 0;
            for(k = 0; k < A->ySize; k++){
                C->matrix[i][j] += A->matrix[i][k] * B->matrix[k][j]; 
            }   

        }
    }
}

または組み合わせparallel for表記を使用する

#pragma omp parallel for private(i, j, k)
for( i = 0; i < A->xSize; i++ ) {
    ...
}

また、ここで OpenMP に複数のスレッドを使用するように指示していることを確認してください。これは、omp_set_num_threads(<number of threads here>)のような環境変数を使用して設定することによって行うことができますOMP_NUM_THREADS。

並列化されることを願っています。:)

于 2013-05-23T12:09:39.327 に答える
1

このコードを使用すると、4 コアでわずかに高速な結果が得られます。

    omp_set_num_threads(4);
    #pragma omp parallel for
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            c[i] += b[j] * a[j][i];
        }
    }

フルプログラム

#include <stdio.h>
#include <time.h>
#include <omp.h>
#include <stdlib.h>


int main() {
    int i, j, n, a[719][719], b[719], c[719];

    clock_t start = clock();

    n = 100; //Max 719

    printf("Matrix A\n");

    for (i = 0; i < n; ++i) {
        for (j = 0; j < n; ++j) {
            a[i][j] = 10;
            printf("%d ", a[i][j]);
        }
        printf("\n");
    }

    printf("\nMatrix B\n");

#pragma omp parallel private(i) shared(b)
    {
#pragma omp for
        for (i = 0; i < n; ++i) {
            b[i] = 5;
            printf("%d\n", b[i]);
        }
    }

    printf("\nA * B\n");

#pragma omp parallel private(i) shared(c)
    {
#pragma omp for
        for (i = 0; i < n; ++i) {
            c[i] = 0;
        }
    }

#pragma omp parallel private(i,j) shared(n,a,b,c)
    {
#pragma omp for schedule(dynamic)
        for (i = 0; i < n; ++i) {
            for (j = 0; j < n; ++j) {
                c[i] += b[j] * a[j][i];
            }
        }
    }


#pragma omp parallel private(i) shared(c)
    {
#pragma omp for
        for (i = 0; i < n; ++i) {
            printf("%d\n", c[i]);
        }
    }

    clock_t stop = clock();
    double elapsed = (double) (stop - start) / CLOCKS_PER_SEC;
    printf("\nTime elapsed: %.5f\n", elapsed);
    start = clock();
    printf("Matrix A\n");

    for (i = 0; i < n; ++i) {
        for (j = 0; j < n; ++j) {
            a[i][j] = 10;
            printf("%d ", a[i][j]);
        }
        printf("\n");
    }

    printf("\nMatrix B\n");

#pragma omp parallel private(i) shared(b)
    {
#pragma omp for
        for (i = 0; i < n; ++i) {
            b[i] = 5;
            printf("%d\n", b[i]);
        }
    }
    printf("\nA * B\n");
    omp_set_num_threads(4);
#pragma omp parallel for
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            c[i] += b[j] * a[j][i];
        }
    }
    stop = clock();
    elapsed = (double) (stop - start) / CLOCKS_PER_SEC;
    printf("\nTime elapsed: %.5f\n", elapsed);
    return 0;
}

最初の方法は

経過時間: 0.03442

2番目の方法

経過時間: 0.02630

于 2016-09-11T16:09:38.533 に答える