c语言数据构造稀疏矩阵

C语言数据构造稀疏矩阵

vl.O可编辑可修改

实验十稀疏矩阵

#include

#defineOKl

#defineERRORO

#defineTRUEl

#defineFALSEO

#defineMAXSIZE100

typedefintStatus;

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,int

row,intcol,ElemType*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;}

StatusFastTransposeSMatrix(constTSMatrix

vl.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(MJ,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",