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 矩阵作为矩阵乘法的结果，但它给出了更大的矩阵和很多垃圾元素它给出的垃圾只有一些值是正确的，但有很多无用的错误元素。

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 function 属性错了吗？我是 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);

并打印输出向量：

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

c中的矩阵乘法

matrix multipication in c

c

void

matrix-multiplication