该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
哪位大佬可以看一下矩阵相乘的代码,这个看不懂,帮忙注释一下,谢谢!
int MultSMatrix(TSMatrix M, TSMatrix N, TSMatrix *Q)//矩阵的相乘
{ int arow, tp,p,brow,t,q,ccol,k;
int ctemp[1000];
if(M.nu != N.mu)
{ printf("对不起,这两个矩阵不能相乘\n"); return 0; }
Q->mu = M.mu; Q->nu = N.nu; Q->tu = 0;//Q矩阵初始化
if( M.tu * N.tu != 0) {
for( arow=1; arow<=M.mu; ++arow) {
for(k=0; k<1000; k++) {
ctemp[k] = 0; }
//Q->rpos[arow] = Q->tu + 1;
if( arow < M.mu ) {
tp = M.rpos[arow+1] ;
}
else {
tp = M.tu + 1;
}
for(p=M.rpos[arow]; p
brow = M.data[p].j;
if(brow < N.mu) t = N.rpos[brow+1];
else { t= N.tu + 1;}
for(q=N.rpos[brow]; q
ccol = N.data[q].j;
ctemp[ccol] += M.data[p].e * N.data[q].e; } }
for(ccol = 1; ccol<=Q->nu; ++ccol)
{ if( ctemp[ccol]) {
++Q->tu;
Q->data[Q->tu].i = arow;
Q->data[Q->tu].j = ccol;
Q->data[Q->tu].e = ctemp[ccol]; }
} } }
return 1;
}