int *rowStart=new int[b.mu+1];
ElemType *temp=new ElemType[b.nu];
for(i=0;i
for(i=1;i<=b.mu;i++) rowStart[i]=rowStart[i-1]+rowSize[i-1];
int current=0,lastInresult=-1;
while(current
for(i=0;i
for(i=rowStart[colA];i
temp[colB]+=data[current].e*b.data[i].e;}
current++;}
for(i=0;i
{lastInresult++;
c.data[lastInresult].ii=rowA;
c.data[lastInresult].jj=i;
c.data[lastInresult].e=temp[i];
}
}
c.mu=mu;c.nu=b.nu;c.tu=lastInresult+1;
delete []rowSize;delete []rowStart;delete []temp;
}}
类的调用如下
[cpp]
#include "stdafx.h"
//稀疏矩阵相关操作的测试
//三元组顺序表的输出
void print(TSMatrix a)
{cout<<" i j e\n";
for(int k=0;k
void main()
{cout<<"运行结果:\n";
int b[N][L],i,j,k=0;
int a[L][N]={0,12,9,0,0,0,0,0,0,0,0,0,0,0,-3,0,0,0,0,14,0,
0,0,24,0,0,0,0,0,18,0,0,0,0,0,15,0,0,-7,0,0,0};
TSMatrix a1,b1,c1;
a1.mu=L;a1.nu=N;a1.tu=0;
cout<<"稀疏矩阵a1:\n";
for(i=0;i
//创建三元组顺序表
for (i=0;i
{ a1.data[k].ii=i;
a1.data[k].jj=j;
a1.data[k].e=a[i][j];
a1.tu++;k++;}
cout<<"原三元组顺序表a1:\n";
print(a1);
cout<<"按任意键继续...\n";getch();
a1.TrMatrix(c1);
for (i=0;i
cout<<"转置矩阵 c1:\n";
for (i=0;i
cout<<"转置矩阵三元组顺序表:\n";
print(c1);
cout<<"按任意键继续...\n";getch();
a1.FastTrMatrix(b1);
for (i=0;i
cout<<"转置矩阵 b1:\n";
for (i=0;i
int p[3][4]={{10,0,5,7},{2,1,0,0},{3,0,4,0}};
int w[4][2]={{2,0},{4,8},{0,14},{3,5}};
TSMatrix aa,bb,cc;
aa.mu=3;aa.nu=4;aa.tu=0;
bb.mu=4;bb.nu=2;bb.tu=0;
//创建三元组顺序表
for (k=0,i=0;i<3;i++)
for (j=0;j<4;j++)
if (p[i][j]!=0)
{ aa.data[k].ii=i;aa.data[k].jj=j;
aa.data[k].e=p[i][j];
aa.tu++;k++;}
//创建三元组顺序表
for (k=0,i=0;i<4;i++)
for (j=0;j<2;j++)
if (w[i][j]!=0)
{ bb.data[k].ii=i;bb.data[k].jj=j;
bb.data[k].e=w[i][j];
bb.tu++;k++;}
aa.mulmatrix(bb,cc);
cout<<"乘积矩阵三元组顺序表:\n";
print(cc);
cout<<"创建的稀疏矩阵:\n";
TSMatrix ff(MROW,MCOL,4);
cout<<"按任意键结束!\n";getch();}
效果如下