c - c 中的矩阵乘法

Question

我知道有类似的问题，但我尝试了那里的建议，但它仍然不起作用，这里是我的完整代码；

#include <stdio.h>
#include <math.h>
void ludec(int N, double lu[N][N]){
    double temp;
    int k,i,j;
    for(k=0;k<N;k++){
        for(i=k+1;i<N;i++){
            temp=lu[i][k]/=lu[k][k];
            for(j=k+1;j<N;j++){
                lu[i][j]-=temp*lu[k][j];
            }
        }
    }
}
void matmul(int u,int m,int n,double p[u][m],double bb[m][n]){
    double mul[u][n];

    int ii,jj,kk;
    for(ii=0;ii<u;ii++){
        for(jj=0;jj<n;jj++){
            mul[ii][jj]=0.0;
        }
    }
    for(ii=0;ii<u;ii++){
        for(jj=0;jj<n;jj++){
            for(kk=0;kk<m;kk++){
                mul[ii][jj]=mul[ii][jj]+p[ii][kk]*bb[kk][jj];
                printf("la %f ii %d jj%d kk%d\n",mul[ii][jj],ii,jj,kk);
            }
        }
    }

}

int main(void){
    double lu[3][3]={{4.0,3.0,-2.0},{-1.0,-1.0,3.0},{2.0,-1.0,5.0}};
    double a[3][3]={{4.0,3.0,-2.0},{-1.0,-1.0,3.0},{2.0,-1.0,5.0}};
    double y[3][1];
    double b[3][1]={{9.0},{-4.0},{6.0}};
    double x[3][1];
    int i,j,N=3;
    double f=0;
    ludec(N, lu);
    y[0][0]=b[0][0];
    for(i=1;i<3;i++){
        for(j=0;j<i;j++){
             f+=lu[i][j]*y[j][0];

         }
         y[i][0]=(b[i][0]-f);
         f=0;
     }
     f=0;



     x[2][0]=y[2][0]/lu[2][2];
     for(i=1;i>=0;i--){

        for(j=1+i;j<3;j++){
             f+=lu[i][j]*x[j][0];

        }
         x[i][0]=(y[i][0]-f)/lu[i][i];
         f=0;

     }
     printf("%f\n %f\n %f\n",x[0][0],x[1][0],x[2][0]);      
    for(i=0;i<3;i++){
        for(j=0;j<3;j++){
            printf("%f %d %d\n",lu[i][j],i,j);
        }
    }



    double mul[3][1];
    matmul(3,3,1,a,x);
    for(i=0;i<3;i++){
        for(j=0;j<2;j++){
            printf("%f %d %d\n",mul[i][j],i,j);
        }
    }

    return 0.0;
}

这里是乘法负责的部分；问题是
我必须得到 3*1 矩阵作为矩阵乘法的结果，但它给出了带有大量垃圾元素的更大矩阵，我在 matmul 函数中通过 printf 监视矩阵乘法，它似乎有效，但是当我在 matmul 函数之外执行相同操作时它只会产生垃圾，只有一些值是正确的，而有很多无用的错误元素。

void matmul(int u,int m,int n,double p[u][m],double bb[m][n]){
        double mul[u][n];

    int ii,jj,kk;
    for(ii=0;ii<u;ii++){
        for(jj=0;jj<n;jj++){
            mul[ii][jj]=0.0;
        }
    }
    for(ii=0;ii<u;ii++){
        for(jj=0;jj<n;jj++){
            for(kk=0;kk<m;kk++){
                mul[ii][jj]=mul[ii][jj]+p[ii][kk]*bb[kk][jj];
                printf("la %f ii %d jj%d kk%d\n",mul[ii][jj],ii,jj,kk);
            }
        }
    }

}

问题是 matmul 函数，我不知道为什么它不起作用，我已经用零初始化了矩阵。我使用 void 函数属性是否错误？我对 C 很陌生，因此我可能做错了什么。整个代码用于 lu 分解，我想通过将结果与初始矩阵相乘来查看代码是否有效。

this is the out put

Answer 1

mul 是 matmul 中的局部变量，因此结果对调用函数不可见。您应该修改 matmul 函数以将结果矩阵作为参数:

void matmul(int u, int m, int n,
            double mul[u][n], const double p[u][m], const double bb[m][n])
{
    for (int i = 0; i < u; i++) {
        for (int j = 0; j < n; j++) {
            mul[i][j] = 0.0;
            for (int k = 0; k < m; k++) {
                mul[i][j] += p[i][k] * bb[k][j];
            }
        }
    }
}

从main调用为

    double mul[3][1];
    matmul(3, 3, 1, mul, a, x);

并打印输出 vector :

for (i = 0; i < 3; i++) {
    for (j = 0; j < 1; j++) {
        printf("%f %d %d\n", mul[i][j], i, j);
    }
}