我正在研究使用分而治之方法执行矩阵乘法的代码。我的代码如下：
#include <stdio.h>
#include <time.h>
#include <stdlib.h>
#define SIZE 8 //Able to change size
void add(int **a, int **b, int size, int **c)
{
int i, j;
for(i=0; i<size; i++)
{
for(j=0; j<size; j++)
{
c[i][j]=a[i][j]+b[i][j];
}
}
}
void subtract(int **a, int **b, int size, int **c)
{
int i, j;
for(i=0; i<size; i++)
{
for(j=0; j<size; j++)
{
c[i][j]=a[i][j]-b[i][j];
}
}
}
void multiply(int **c, int **d, int size, int **result)
{
if(size==1)
{
result[0][0]=c[0][0]*d[0][0]; /*Could this base case calculation be wrong?*/
}
else
{
int i, j;
int new_size=size/2;
int **c11=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
c11[i]=(int *)malloc(new_size*sizeof(int));
}
int **c12=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
c12[i]=(int *)malloc(new_size*sizeof(int));
}
int **c21=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
c21[i]=(int *)malloc(new_size*sizeof(int));
}
int **c22=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
c22[i]=(int *)malloc(new_size*sizeof(int));
}
int **d11=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
d11[i]=(int *)malloc(new_size*sizeof(int));
}
int **d12=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
d12[i]=(int *)malloc(new_size*sizeof(int));
}
int **d21=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
d21[i]=(int *)malloc(new_size*sizeof(int));
}
int **d22=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
d22[i]=(int *)malloc(new_size*sizeof(int));
}
int **temp1=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
temp1[i]=(int *)malloc(new_size*sizeof(int));
}
int **temp2=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
temp2[i]=(int *)malloc(new_size*sizeof(int));
}
int **temp3=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
temp3[i]=(int *)malloc(new_size*sizeof(int));
}
int **temp4=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
temp4[i]=(int *)malloc(new_size*sizeof(int));
}
int **temp5=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
temp4[i]=(int *)malloc(new_size*sizeof(int));
}
int **temp6=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
temp6[i]=(int *)malloc(new_size*sizeof(int));
}
int **temp7=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
temp7[i]=(int *)malloc(new_size*sizeof(int));
}
int **temp8=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
temp8[i]=(int *)malloc(new_size*sizeof(int));
}
int **res1=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
res1[i]=(int *)malloc(new_size*sizeof(int));
}
int **res2=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
res2[i]=(int *)malloc(new_size*sizeof(int));
}
int **res3=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
res3[i]=(int *)malloc(new_size*sizeof(int));
}
int **res4=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
res4[i]=(int *)malloc(new_size*sizeof(int));
}
for(i=0; i<new_size; i++)
{
for(j=0; j<new_size; j++)
{
c11[i][j]=c[i][j];
c12[i][j]=c[i][j+new_size];
c21[i][j]=c[i+new_size][j];
c22[i][j]=c[i+new_size][j+new_size];
d11[i][j]=d[i][j];
d12[i][j]=d[i][j+new_size];
d21[i][j]=d[i+new_size][j];
d22[i][j]=d[i+new_size][j+new_size];
}
}
multiply(c11, d11, new_size, temp1);
multiply(c12, d21, new_size, temp2);
multiply(c11, d12, new_size, temp3);
multiply(c12, d22, new_size, temp4);
multiply(c21, d11, new_size, temp5);
multiply(c22, d21, new_size, temp6);
multiply(c21, d12, new_size, temp7);
multiply(c22, d22, new_size, temp8);
add(temp1, temp2, new_size, res1);
add(temp3, temp4, new_size, res2);
add(temp5, temp6, new_size, res3);
add(temp7, temp8, new_size, res4);
int x=0;
int y=0;
int z=0;
int w=0;
for(i=0; i<2*new_size; i++)
{
for(j=0; j<2*new_size; j++)
{
if(i<new_size && j<new_size)
{
result[i][j]=res1[i][j];
}
if(i<new_size && j>=new_size)
{
x=j-new_size;
result[i][j]=res2[i][x];
}
if(i>=new_size && j<new_size)
{
y=i-new_size;
result[i][j]=res3[y][j];
}
if(i>=new_size && j>=new_size)
{
z=i-new_size;
w=j-new_size;
result[i][j]=res4[z][w];
}
}
}
free(c11);
free(c12);
free(c21);
free(c22);
free(d11);
free(d12);
free(d21);
free(d22);
free(temp1);
free(temp2);
free(temp3);
free(temp4);
free(temp5);
free(temp6);
free(temp7);
free(temp8);
free(res1);
free(res2);
free(res3);
free(res4);
}
}
int main(void)
{
int i, j;
int **A=(int **)malloc(SIZE*sizeof(int *));
for(i=0; i<SIZE; i++)
{
A[i]=(int *)malloc(SIZE*sizeof(int));
}
int **B=(int **)malloc(SIZE*sizeof(int *));
for(i=0; i<SIZE; i++)
{
B[i]=(int *)malloc(SIZE*sizeof(int));
}
int **C=(int **)malloc(SIZE*sizeof(int *));
for(i=0; i<SIZE; i++)
{
C[i]=(int *)malloc(SIZE*sizeof(int));
}
int count=0;
int sum=0;
srand(time(NULL));
for(i=0; i<SIZE; i++)
{
for(j=0; j<SIZE; j++)
{
A[i][j]=rand()%1000;
}
}
for(i=0; i<SIZE; i++)
{
for(j=0; j<SIZE; j++)
{
B[i][j]=rand()%1000;
}
}
for(i=0; i<SIZE; i++)
{
for(j=0; j<SIZE; j++)
{
C[i][j]=0;
}
}
for(i=0; i<SIZE; i++)
{
for(j=0; j<SIZE; j++)
{
printf("%5d", A[i][j]);
}
printf("\n");
}
printf("\n");
for(i=0; i<SIZE; i++)
{
for(j=0; j<SIZE; j++)
{
printf("%5d", B[i][j]);
}
printf("\n");
}
multiply(A, B, SIZE, C);
for(i=0; i<SIZE; i++)
{
for(j=0; j<SIZE; j++)
{
printf("%5d", C[i][j]);
}
printf("\n");
}
free(A);
free(B);
free(C);
}
该代码用于生成随机的两个n * n矩阵，将大矩阵拆分为一半，然后按您可能知道的方式递归地乘以子矩阵。计算后，我将显示乘法结果。我将计算出的子矩阵放在结果矩阵中。
我尝试使用调试器，但是程序只是在片刻后退出。根据我以前的经验，我怀疑变量的初始化可能会出现一些问题，但是目前看来对我来说还不错。
谁能指出我错过了什么？对不起，中间有多余的malloc。发现代码正常工作后，我将创建一个函数。预先非常感谢。
答案 0 :(得分：0)
您的代码包含：
int **temp5=(int **)malloc(new_size*sizeof(int *));
for(i=0; i<new_size; i++)
{
temp4[i]=(int *)malloc(new_size*sizeof(int));
}
temp4是一个错字； temp5[i]的内容从未设置。