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 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中的矩阵乘法

1 个答案: