spmatrix *Chen(spmatrix *a,spmatrix *b)
{
int ano,bno,cno,col;
spmatrix * c;
c=InitSpmatrix(a->m,b->n,(a->m)*(b->n));
c->t=0;
cno=0;
for(col=0;col<a->n;col++){
bno=0;
for(ano=0;ano<a->t;ano++){
if(a->data[ano].j==b->data[bno].i){
c->data[cno].i=a->data[ano].i;
c->data[cno].j=b->data[bno].j;
c->data[cno].v+=a->data[ano].v*b->data[bno].v;
}
bno++;
}
if(c->data[cno].v!=0)
(c->t)++;
cno++;
}
return c;
}
//关键函数,实现矩阵相乘,存入三元矩阵C