1

私の行列乗算コードは

int matMul(int ld, double** matrix)
{ 

  //local variables initialize

  omp_set_num_threads(nthreads);


  \#pragma omp parallel private(tid,diag,ld) shared(i,j,k,matrix)

  {
    /* Obtain and print thread id */

    tid = omp_get_thread_num();

    for ( k=0; k<ld; k++)  {
    if (matrix[k][k] == 0.0) {
      error = 1;
      return error;
    }
    diag = 1.0 / matrix[k][k];
\#pragma omp for 

    for ( i=k+1; i < ld; i++) {

      matrix[i][k] = diag * matrix[i][k];

    }
    for ( j=k+1; j<ld; j++) {

      for ( i=k+1; i<ld; i++) {

        matrix[i][j] = matrix[i][j] - matrix[i][k] * matrix[k][j];

      }

    }
  } 

  }  
  return error;

}

マトリックスオブジェクトのみが原因だと思いますが、パラメーターとして渡されているのになぜnullになるのでしょうか..

4

2 に答える 2

1

I ran into the same problem while compiling my code under Linux using GCC 4.2. The line that is causing the problem is: 

omp_set_num_threads(nthreads);

You should try to set the number of threads by specifying it under pragma omp for: 

#pragma omp for num_threads(nthreads)

Hope it helps!

于 2010-07-12T19:49:57.677 に答える
0

あなたの「マトリックス」は、実際のマトリックスの行へのポインターの配列であると想定しています。

double *matrix[NROWS];
for (i = 0; i < NROWS; ++i) {
    matrix[i] = malloc(sizeof(double)*NCOL);
}

しかし、「行列」がこのように定義されている場合

double matrix[NROWS][NCOL];

あなたのプログラムは動作しません。マトリックス ID が適切に定義されている場合、クラッシュの原因として考えられるのは、不適切なサイズ ("ld") である可能性があります。

よろしく

于 2010-03-27T10:13:34.067 に答える