[发明专利]一种复数矩阵的优化方法

权利要求书

1.一种在龙芯上的复数矩阵的优化方法，其特征在于操作步骤如下：

第一步：利用展开粒度公式K_AK_B/（K_A+K_B）=m/n计算最优取数比，其中K_A和KB分别为矩阵A和矩阵B的展开粒度，m/n为最优运算取数比，令两个矩阵展开粒度的乘积除以两个矩阵的展开粒度之和所得的值与最优取数比相同，在龙芯上调整两个矩阵的展开粒度KA和KB，使得在循环结构中保证矩阵A和矩阵B的运算取数比m∶n为2∶1，则所得的展开粒度即为在运算中一次循环需要展开的运算次数；

第二步：在矩阵运算时，对于矩阵内数据大于500的大规模矩阵计算，对每一分块进行运算时在循环中采用4×4的展开，每个展开的运算对应龙芯的一个发射槽，通过发射槽发射运算指令给龙芯运算器进行乘加运算；根据矩阵分块公式24×(nb)²≤NK，其中NK为龙芯的一级数据缓存cache的大小，将矩阵分块大小nb选取为一个小于52并使得24乘以矩阵分块大小nb的平方（nb）²小于龙芯处理器一级数据缓存64千字节的最大的值，得到在龙芯上最优的理想矩阵分块大小；

第三步：对于在复数矩阵运算中出现的普通复数乘加运算，利用复数经典算法中两个复数之间的乘法算法，设第一复数为A=x+yi，第二复数为B=m+ni，第三复数为C=p+qi，计算第一复数A乘以第二复数B后再加第三复数C的值，将此值赋给第三复数作为第三复数的新值，算法流程如下：计算第一复数A的实部x乘以第二复数B的实部m的值为第一实数e←x×m，计算第一复数A的虚部y乘以第二复数B的虚部n的值为第二实数f←y×n，计算第三复数C的实部p加第一实数e后减去第二实数f的值，并将此值付给第三复数C的实部p←p+e-f，计算第一复数A的实部x加第一复数A的虚部y的值为第三实数k←x+y，计算第二复数B的实部m加第二复数B的虚部n的值为第四实数t←m+n，计算第三复数C的虚部q加第三实数k乘以第四实数t的积再减去第一实数e和第二实数f的值，并将此值付给第三实数k的虚部q←q+k×t-e-f，则得到第三复数C的新值，通过龙芯发射槽传给运算器参与分块矩阵的运算。