[发明专利]基于异构众核处理器的格点量子色动力学并行加速方法

权利要求书

1.一种基于异构众核处理器的格点量子色动力学并行加速方法，其特征在于包括有下列步骤：

步骤一，异构众核处理器的从核矩阵位置初始化；

由于异构众核处理器中存在有多个从核，因此需要对从核按照从核标识号进行位置划分，并记录出每个从核的矩阵位置；

从核集CPEs＝{cpe₁,cpe₂,…,cpe_A}按照8×8的矩阵进行位置排序，得到从核集位置矩阵add^CPEs，任意一个从核位置记为d_p,q：

步骤二，主核读取费米子场量和规范场量；

步骤201，主核MPE读取数据信息，将读取的所有数据信息以集合形式表达为SPM＝{S₁,S₂,…,S_g,…,S_G}；

步骤202，主核将所述的SPM＝{S₁,S₂,…,S_g,…,S_G}数据信息中的费米子场量按读取的先后顺序存入一个8×8×8×8的格点矩阵DAA^MPE中，并将DAA^MPE保存到存储器中；

表示第一个数据信息S₁在四维坐标点上的费米子场量值；表示S₁的费米子场量；

表示第二个数据信息S₂在四维坐标点上的费米子场量值；表示S₂的费米子场量；

表示任意一个数据信息S_g在四维坐标点上的费米子场量值；表示S_g的费米子场量；

表示最后一个数据信息S_G在四维坐标点上的费米子场量值；表示S_G的费米子场量；

步骤203，主核将所述的SPM＝{S₁,S₂,…,S_g,…,S_G}数据信息中的规范场量按读取的先后顺序存入一个4×8×8×8×8的格点链接矩阵DBB^MPE中，并将DBB^MPE保存到存储器中；

表示第一个数据信息S₁在四维坐标点的方向上的规范场量值；表示S₁在四维坐标点的方向；表示S₁的规范场量；

表示第二个数据信息S₂在四维坐标点的方向上的规范场量值；表示S₂在四维坐标点的方向；表示S₂的规范场量；

表示任意一个数据信息S_g在四维坐标点的方向上的规范场量值；表示S_g在四维坐标点的方向；表示S_g的规范场量；

表示最后一个数据信息S_G在四维坐标点的方向上的规范场量值；表示S_G在四维坐标点的方向；表示S_G的规范场量；

步骤三，从核基于自身的行号列号读取数据信息实现数据分割；

步骤301，任意从核cpe_A依据步骤一的从核矩阵位置d_p,q，按照Z轴、T时间轴方向将DAA^MPE矩阵中cpe_A负责的数据信息部分读入局部存储空间，记为则所有从核读入的数据信息可记为：

步骤302，任意从核cpe_A依据步骤一的从核矩阵位置d_p,q，按照Z轴、T时间轴方向将DBB^MPE矩阵中cpe_A负责的数据信息部分读入局部存储空间，记为则所有从核读入的数据信息可记为：

步骤四，对任意一个从核中任意一个格点的数据信息进行运算；

步骤401，任意一个数据信息S_g的格点费米子场量记为执行步骤403；

步骤402，任意一个数据信息S_g的规范场量记为执行步骤403；

步骤403，从任意一个格点的数据信息S_g的x,y,z,t四个维度获取相邻8个格点的数据信息，然后获取相邻8个格点的格点费米子场量和规范子场量；执行步骤404；

相邻8个格点的数据信息分别记为S₁、S₂、S₃、S₄、S₅、S₆、S₇和S₈，位于相邻8个格点的中心格点是S_g，则格点费米子场量分别记为和所述S_g的格点费米子场量记为

则规范场量分别记为和

第一个数据信息S₁的费米子场量记为

第一个数据信息S₁的规范场量记为

第二个数据信息S₂的费米子场量记为

第二个数据信息S₂的规范场量记为

第三个数据信息S₃的费米子场量记为

第三个数据信息S₃的规范场量记为

第四个数据信息S₄的费米子场量记为

第四个数据信息S₄的规范场量记为

第五个数据信息S₅的费米子场量记为

第五个数据信息S₅的规范场量记为

第六个数据信息S₆的费米子场量记为

第六个数据信息S₆的规范场量记为

第七个数据信息S₇的费米子场量记为

第七个数据信息S₇的规范场量记为

第八个数据信息S₈的费米子场量记为

第八个数据信息S₈的规范场量记为

步骤404，将相邻8个格点的费米子场量和规范子场量进行矩阵相乘；执行步骤405；

步骤405，以相邻8个格点的矩阵相乘量更新中心格点是S_g的格点费米子场量，更新后属于S_g的格点费米子场量，记为且即中的数据信息S_g被更新为执行步骤五；

步骤五，每个从核对其局部存储空间中的每个格点的数据信息都进行步骤四的并行处理，从而得到更新后的所有格点的费米子场量，即获得执行步骤六；

步骤六，更新完成后，迭代次数加1；并计算格点费米子场量的残差值；

迭代次数记为U，最大迭代次数记为U_max，且U_max取值为1000，当前迭代次数记为U_当前；若U_当前＜U_max，则执行步骤四；若U_当前≥U_max，则执行步骤七；

格点费米子场量的残差记为R，格点费米子场量的残差阈值记为R_min，且R_min的取值为1.0×10^－12；若R＞R_min，则执行步骤四；若R≤R_min，则执行步骤七；

步骤七，将更新完成后的格点矩阵输出给存储器，作为文件保存；

将传给存储器，更新DAA^MPE得到将保存并写入文件。