有关矩阵相乘的问题

[wmjcom]

void CSelfTuningPID::Matrix_Multiply(double *A, double *B, int m, int n, int k, double *C)//矩阵A、B相乘算法:C=A*B
{
int i,j,l;
for(i=0;i<m;i++)
{
for(j=0;j<k;j++)
{
*(C+i*k+j)=0.0;
for(l=0;l<n;l++)
{
*(C+i*k+j)=*(C+i*k+j)+(*(A+i*n+l))*(*(B+l*k+j));
}
}

}
}
能否帮忙看一下有什么问题？总出错，没找到问题。多谢！

[upuper]

up

[carolbaby]

矩阵相乘的公式我也忘差不多了，楼主的程序也没仔细看，不过，感觉传入的参数好像不太对劲？是不是应该是double**? 不可能3个矩阵都是一维的吧。

[wmjcom]

就是为了采用一维指针实现二维数组运算，多谢，请各位帮忙！