有部 templ=0 ∥当前行各元素累加器清零 Q. rposlarow=Qtu+l for(p=M. rposlarow p<M rpos arow+1; ++p)i /对当前行中每一个非零元 处brow= M. datalp] j IE if(brow Nnu )t-N.rposlbrowtlf3' M else t=Ntu+1 的每 for(q= N rpos[brow]; q<t;++qi ccO l= N datal,∥乘积元素在Q中列号 ctemplccol] +=M.datalp) e*N datale }∥forq 行}∥求得Q中第 crow =arow )行的非零元 for(ccol=; ccok=Q nu; ++ccol) if (ctemplccollt if(++Q tu> MAXSIZE) return ERROR Q data[Q tu]-=arow, ccol, temp[ccol 南京航空航天大学数据结构课题组版权所有ctemp[] = 0; // 当前行各元素累加器清零 Q.rpos[arow] = Q.tu+1; for (p=M.rpos[arow]; p<M.rpos[arow+1];++p) { //对当前行中每一个非零元 brow=M.data[p].j; if (brow < N.nu ) t = N.rpos[brow+1]; else { t = N.tu+1 } for (q=N.rpos[brow]; q< t; ++q) { ccol = N.data[q].j; // 乘积元素在Q中列号 ctemp[ccol] += M.data[p].e * N.data[q].e; } // for q } // 求得Q中第crow( =arow)行的非零元 for (ccol=1; ccol<=Q.nu; ++ccol) if (ctemp[ccol]) { if (++Q.tu > MAXSIZE) return ERROR; Q.data[Q.tu] = {arow, ccol, ctemp[ccol]}; } // 处 理 的 每 一 行 M