c語(yǔ)言數(shù)據(jù)構(gòu)造稀疏矩陣

C語(yǔ)言數(shù)據(jù)構(gòu)造稀疏矩陣vl.O可編輯可修改實(shí)驗(yàn)十稀疏矩陣#include#defineOKl#defineERRORO#defineTRUEl#defineFALSEO#defineMAXSIZE100typedefintStatus;typedeffloatElemType;typedefstruct{returnERROR;)p=l;=row;M->data[p].j=col;M->data[p].e=e;M->tu++;returnOK;}for(t=l;ttu;t++)&&(col>=M->data[t].j))p++;if(row==M->data[t-l].i&&col==M->data[t-l].j){=e;returnOK;}for(i=M->tu;i>=p;i-){=M->data[i].i;M->data[i+l].j=M->data[i].j;M->data[i+l].e=M->data[i].e;}=row;M->data[p].j=col;M->data[p].e=e;M->tu++;returnOK;}StatusFindElem(constTSMatrix*M,introwjntcol,日emType*e)(==row&&M->data[p].j==col)(*e=M->data[p].e;returnTRUE;)returnFALSE;}StatusTransposeSMatrix(constTSMatrix*M,TSMatrix*T){==col){T->data[q].i=M->data[p].j;T->data[q].j=M->data[p].i;T->data[q].e=M->data[p].e;q++；))returnOK;}StatusFastTransposeSMatrixfconstTSMatrixvl.O可編輯可修改*M,TSMatrix*T)(n")；exit(O);)for(col=l;colnu;col++)num[col]=0;];cpot[l]=l;;q=cpot[col];T->data[q].i=M->data[p].j;T->data[q].j=M->data[p].i;T->data[q].e=M->data[q].e;++cpot[col];}n")；returnERROR;}Q->mu=M->mu;Q->nu=T->nu;Q->tu=0;P=l；for(i=l;imu;i++){for(j=l;jnu;j++){s=0;for(k=l;knu;k++){if(FALSE==FindElem(M,i,k,&m))continue;if(FALSE==FindElem(T,k,j,&t))continue;s+=m*t;}if(s!=O){=i;Q->data[p].j=j;Q->data[p].e=s;P++；Q->tu++;)))returnOK;)voidPrintSMatrix(constTSMatrix*M){&&j==M->data[p].j)printf("%g\t",