我使用以下两个代码块串行和并行计算矩阵乘法。
串行-
double** ary1 = new double*[in];
double** ary2 = new double*[in];
double** result = new double*[in];
for (int i=0;i<in;i++){
for (int j=0;j<in;j++){
result[i][j] = 0;
for(int k = 0;k<in; k++){
result[i][j] += ary1[i][k]*ary2[k][j];
}
}
}
并行-
double** ary1 = new double*[in];
double** ary2 = new double*[in];
double** resultsP = new double*[in];
#pragma omp parallel for
for(int i=0;i<size;i++){
int raw = i/in;
int column = i%in;
double sum =0;
for(int k = 0; k < in; k++){
resultsP[raw][column] += ary1[raw][k]*ary2[k][column];
}
resultsP[raw][column] = sum;
}
我在四核计算机上运行了这两个程序,但得到了相同的结果。为什么我不能通过模仿来提高性能?在并行循环中访问ary1、ary2、resultsP共享数组是否会导致它们串行运行?
由于编译代码时未包含'-fopenmp'标志,因此发生了这种情况。通过添加它解决了问题。