C openMP 程序的完全空白输出

Question

我编写了以下矩阵乘法程序:

#include <stdio.h>
#include <omp.h>
#include <time.h>
#define NRA 500
#define NCA 500
#define NCB 500

int mat_mul() ;
int i , j , k ;
int main(){
    double start , end ;
    start  =  omp_get_wtime() ;
    mat_mul() ;
    end = omp_get_wtime() ; 
    printf("Time taken : \n %lf " , (end - start) );
    return 0;
}

int mat_mul(){
    int mat1[NRA][NCA] , mat2[NCA][NCB] , mat3[NRA][NCB] ;
    //double start , end ;
    //start  =  omp_get_wtime() ;
    #pragma omp parallel private( i , j , k ) shared(mat1 , mat2 , mat3)
    #pragma omp for

    for( i = 0 ; i < NRA ; i++){
        for( j = 0 ; j < NCB ; j++){
            mat1[i][j] = mat2[i][j] = rand() ;
            for ( k = 0 ; k <  NCB ; k++){
                mat3[i][k] += mat1[i][k] * mat2[k][j] ;
            }
        }
    }
    //end = omp_get_wtime() ; 
    printf("REsult : \n");
    for( i = 0 ; i < NRA ; i ++ ) {
        for( j =0 ; j < k ; j ++)
            printf("%lf" , (double)mat3[i][j]);
        printf("\n") ;
    //printf("Time taken : \n %lf " , (end - start) );
    }
    return 0 ;
}

一切都很好(几乎):它编译，执行，甚至终止:-D(并且它实际上提供了加速(比较。

但不幸的是，输出看起来像这样:

[tejas@localhost Documents]$ gcc -pedantic -Wall  -std=c99 -fopenmp  par.c
par.c: In function ‘mat_mul’:
par.c:28:30: warning: implicit declaration of function ‘rand’ [-Wimplicit-function-declaration]
    mat1[i][j] = mat2[i][j] = rand() ;
                              ^~~~
[tejas@localhost Documents]$ ./a.out
REsult :
...
(LOTS of blank spaces later)
...
Time taken : 
0.177506 [tejas@localhost Documents]$ 

我做错了什么？ 我正在使用 Fedora 27 、GCC

提前谢谢您。

Answer 1

这是因为 j < k循环条件:

for( j =0 ; j < k ; j ++)

k的值使用 openmp 时有所不同。你可以将 k 替换为 NCB ？