我试图编写一个程序来使用 C 中的 OpenMP 执行矩阵乘法。每当我尝试将字符串文字作为参数传递给函数时,我都会面临分段错误错误。
#include<stdio.h>
#include<time.h>
#include<omp.h>
struct matrix{
int r;
int c;
int mat[1000][1000];
};
void read_matrix(char* fname, struct matrix* m){
FILE *fp;
fp = fopen(fname,"r");
fscanf(fp, "%d %d",&m->r,&m->c);
for(int i=0;i<m->r;i++){
for(int j=0;j<m->c;j++){
fscanf(fp, "%d", &m->mat[i][j]);
}
}
fclose(fp);
}
void write_matrix(char* fname, struct matrix m){
FILE *fp;
fp = fopen(fname,"w");
fprintf(fp,"%d %d\n",m.r,m.c);
for(int i=0;i<m.r;i++){
for(int j=0;j<m.c;j++){
fprintf(fp,"%d\n",m.mat[i][j]);
}
}
fclose(fp);
}
void main(){
struct matrix m1;
struct matrix m2;
struct matrix res;
read_matrix("m1",&m1);
read_matrix("m2",&m2);
int r1 = m1.r;
int c1 = m1.c;
int c2 = m2.c;
res.r = r1;
res.c = c2;
for(int i=0;i<r1;i++){
for(int j=0;j<c2;j++){
res.mat[i][j] = 0;
}
}
#pragma omp parallel
{
#pragma omp for
for(int i = 0; i < r1; i++){
for(int j = 0; j < c2; j++){
for(int k = 0; k < c1; k++){
#pragma omp atomic update
res.mat[i][j] += m1.mat[i][k]*m2.mat[k][j];
}
}
}
}
write_matrix("res",res);
}
代码Segmentation fault (core dumped)
在运行时显示。在 GDB 上运行时,它显示Program received signal SIGSEGV, Segmentation fault. 0x0000000000400ad9 in main () at mm.c:40 40 read_matrix("m1",&m1);
在第一次 read_matrix() 调用之前我添加了一个 printf 语句printf("check\n");
printf 调用现在开始抛出分段错误。我假设传递字符串文字是错误的原因。代码可能有什么问题?
在第一次
read_matrix()
调用之前,我添加了一个 printf 语句printf("check\n");
printf 调用现在开始抛出一个分段错误。
我认为这可能是因为您struct matrix
在堆栈上分配。每个矩阵大约 4MB,所以堆栈上有 12MB。我猜它会创建站点(堆栈溢出)。
尝试对矩阵使用静态变量或动态分配它们。它可能会解决您的问题。如果没有,永远不要在堆栈上分配 12MB 的结构......
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句