第6章二维数组应用实例及概念.ppt

两个稀疏矩阵相乘的乘积不一定是稀疏矩阵；反之，即使A［i, k］×B［k, j］不为零，其累加和Q［i, j］也可能为零。累加和是否为零， 只有在求得相应的累加和后才能得知。只有累加和非零，才能将其作为乘积矩阵Q的三元组表中的一个元素。 由于Q中的行号与A是一致的，且A.item是以行序为主序的，因此可对Q进行逐行处理，并设立累计求和的中间变量temp数组来存储Q的一行值，在累加完毕后，将其中的非零元素写入三元组表Q.item。 稀疏矩阵成员函数function Txsjz.mutiply(b: Txsjz): Txsjz 的功能为返回当前矩阵对象与矩阵b的乘积矩阵，其处理过程为： (1) 求矩阵b中每一行非零元素的个数， 将结果存入数组rnum；  (2) 根据数组rnum求矩阵b中每行第一个非零元素在三元组表中的位置， 将结果存入数组rstart；  (3) 依次扫描当前矩阵的三元组表， 按数组rnum和rstart中指出的位置，将原矩阵非零元素与矩阵b中的相应的元素相乘，每处理完一行后， 将非零元素数据填入乘积矩阵的三元组表。 程序清单如下： function Txsjz.mutiply(b: Txsjz): Txsjz; {稀疏矩阵相乘} var i, k, r, p, ca, cb: integer;  rnum, rstart: array ［1..maxitem］ of integer;  temp: array ［1..maxitem］ of elemtp;  begin if C b.R then 显示出错并返回空矩阵  else begin result:= Txsjz.create; {产生返回对象} for i:=1 to b.R do rnum［i］:=0; {数组清空} for i:=1 to b.num do  rnum［b.item［i］.row］:= rnum［b.item［i］.row］+1; {求矩阵B中每一 行非零元素的个数} rstart［1］:=1;  for i:=2 to b.R do {求矩阵B中每行第一个非零元素的开始位置} rstart［i］＝r num［i-1］＋rstart［i-1］;  k:=1; {k用来指示矩阵A中三元组表当前位置} p:=1; {p用来指示乘积矩阵中三元组表的尾部} while k=num do  begin r:= item［k］.row; {取A的三元组表当前行号} for i:=1 to b.C do temp［i］:=0; {清空累加中间数组} while item［k］.row = r do begin ca:= item［k］.col; {取矩阵A当前的列号} for i:= rstart［ca］ to rstart［ca］+rnum［ca］-1 do begin cb:= b.item［i］.col; {取矩阵B当前的列号} temp［cb］:= temp［cb］+item［k］.value*b.item［i］.value;  end; {A的r行与B的cb列相乘}