[發明專利]一種應用于對稱矩陣與向量乘法的計算方法

有效

• 申請號: 201710671890.2
• 申請日: 2017-08-08
• 公開（公告）號: CN107590106B
• 公開（公告）日: 2020-07-28
• 發明（設計）人: 薛瑞;張浩;范東睿;葉笑春;朱亞濤
• 申請（專利權）人: 北京中科睿芯科技有限公司
• 主分類號: G06F17/16
• 分類號: G06F17/16
• 代理公司: 北京科龍寰宇知識產權代理有限責任公司 11139
• 代理人: 孫皓晨

1.一種應用于對稱矩陣與向量乘法的計算方法，該方法用于計算n1×n1的對稱矩陣與n1維列向量的乘積，其特征在于，所述方法對n1×n1的對稱矩陣進行并行處理，具體包括以下步驟：

S1：讀取n1×n1對稱矩陣中的上三角矩陣中的數據；

S2：以m為邊長對上三角矩陣進行分塊，分塊后的數據塊A_ij排列成z行z列，1≤i≤j≤z，其中，

若n1/m為整數，則z=n1/m，分塊后得到z個m×m三角矩陣塊以及n1（n1-m）/2m²個m×m的普通矩陣塊，

若n1/m不為整數，則z=[n1/m]+1，分塊后得到z-1個m×m三角矩陣塊、1個a×a三角矩陣塊、(z-1)(z-2)/2個m×m普通矩陣塊以及z-1個m×a矩陣塊，其中a= n1-（z-1）×m，[]為取整符號；

S3：分別對S2中得到的z-1個m×m三角矩陣塊A_ii進行微量數據擴展，使之成為對稱矩陣塊；

S4：以m為邊長對n1維列向量進行分塊，分塊后的數據塊B_i1共z行，其中，

若n1/m為整數，則z=n1/m，分塊后得到z個m×1矩陣塊，

若n1/m不為整數，則z=[n1/m]+1分塊后得到z-1個m×1矩陣塊以及1個a×1矩陣塊，其中a= n1-（z-1）×m；

S5：構建一中間數據塊C_z1，根據下式計算C_z1中每一元素的值：

C_i1= A_ii×B_i1+……+ A_iz×B_z1；

S6：計算最終結果向量D_z1，其中：

D₁₁= C₁₁，

D_i1= f（A_1i，B₁₁）+……+ f（A_（i-1）i，B_（i-1）1）+ C_i1 ，其中，2≤i≤z，

;

所述S3步驟之后的上三角矩陣中的數據存儲于第一內存區，且占用第一內存區的存儲空間大小為；所述S4步驟之后的n1維列向量存儲于第二內存區，且占用第二內存區的存儲空間大小為n1×1;中間數據塊C_z1存儲于第三內存區，且占用第三內存區的存儲空間大小為n1×1。

2.根據權利要求1所述的應用于對稱矩陣與向量乘法的計算方法，其特征在于，n1=64，m=8。

3.根據權利要求1所述的應用于對稱矩陣與向量乘法的計算方法，其特征在于，n1=7，m=2。