Question

我正在尝试使用2个文件作为参数（mat1.txt和mat2.txt）来实现矩阵乘法程序。我决定将结果加载到一个文件中，但是当我执行程序时：./m mat1.txt mat2.txt 20 20 20 20（文件名和每个矩阵的行数和列数），我得到并且错误：分段错误（核心倾销）。问题在于创建或写入文件结果，你能帮帮我吗？代码如下：

 /* Memory management */
mem_mat1 = (int *) malloc(M1 * N1 * sizeof(int));
if (mem_mat1 == NULL) {
    fprintf(stderr,"Error: malloc mem_mat1\n");
    /*MPI_Finalize();*/
    return (3);
}

mat1 = (int **) malloc(M1 *sizeof(int *));
if (mat1 == NULL) {
    fprintf(stderr,"Error: malloc mat1\n");
    /*MPI_Finalize();*/
    return (3);
} 

for (i=0; i<M1; i++) {
    mat1[i] = mem_mat1+(i*N1);
}

mem_mat2 = (int *) malloc(M2 * N2 * sizeof(int));
if (mem_mat2 == NULL) {
    fprintf(stderr,"Error: mem_mat2\n");
    /*MPI_Finalize();*/
    return (3);
}

mat2 = (int **) malloc(M2 *sizeof(int *));
if (mat2 == NULL) {
    fprintf(stderr,"Error: malloc mat2\n");
    /*MPI_Finalize();*/
    return (3);
} 

for (i=0; i<M2; i++) {
    mat2[i] = mem_mat2+(i*N2);
}

mem_matR = (int *) malloc(M1 * N2 * sizeof(int));
if (mem_matR == NULL) {
    fprintf(stderr,"Error: malloc mem_matR\n");
    /*MPI_Finalize();*/
    return (3);
}

matR = (int **) malloc(M1 *sizeof(int *));
if (matR == NULL) {
    fprintf(stderr,"Error: malloc matR\n");
    /*MPI_Finalize();*/
    return (3);
}

for (i=0; i<M1; i++) {
    matR[i] = mem_matR+(i*N2);
}

/* SEQUENTIAL PRODUCT MATRIX */

/*Open and Read file 1: mat1.txt*/
fmat1 = fopen(argv[1],"rb");
fread(mat1, M1 *sizeof(int *), N1 *sizeof(int *), fmat1);

/*Open and Read file 2: mat2.txt*/
fmat2 = fopen(argv[2],"rb");
fread(mat2, M2 *sizeof(int *), N2 *sizeof(int *), fmat2);

/*Create a file to write the result fmatR.txt*/
fdR = creat("matR.txt", "w");
if (fdR < 0) {
    fprintf(stderr, "Error create result file\n");
    return(2);
}

for (i=0;i<M1;i++) {
    for (j=0;j<N2;j++) {
        sum=0;
        for (k=0;k<N1;k++) {
            sum+=mat1[i][k]*mat2[k][j];
            matR[i][j]=sum;
        }
    }
 }
 write(fdR, &matR, M1*N2 *sizeof(int));

Answer 1

我认为你在使用matR编写随机内存方面做得很糟糕......最后一行也是错误的。

修正：

1）带走 mem_matR 以及与之相关的所有代码

2）让 matR 成为指向矩阵的第一个元素（你的缓冲区）的指针，让它成为M1xN2元素的简单一维数组

matR = (int *) malloc( M1 * N2 * sizeof(int));

3）在for循环中以这种方式写matR

matR[j+i*N2]=sum;

4）最后以这种方式保存数据

write(fdR, (void *) matR, M1*N2 *sizeof(int));

-

我也认为fread调用是错误的（即使这不会导致问题）

fread(mat1, M1 *sizeof(int *), N1 *sizeof(int *), fmat1);

如果我知道你有M1 x N1元素的整数。所以...

fread(mat1, sizeof(int *), M1 * N1, fmat1);

同样适用于secon电话

如何将内存中的结果放到文件* .txt中

1 个答案: