[发明专利]一种基于辅助树的多层快速多极子并行网格细剖方法

权利要求书

1.一种基于辅助树的多层快速多极子并行网格细剖方法，其特征在于，包括以下步骤：

步骤一、将复杂结构目标剖分为三角形网格，根据三角形网格的顶点编号和顶点坐标计算三角形中心坐标，基于三角形中心，按照多层快速多极子盒子划分原则，生成网格的辅助树结构；

步骤二、根据预设的并行方案和MPI进程数，对网格进行负载分配，确定归属于各进程上的三角形网格编号；

步骤三、对分配后的三角形网格编号、顶点编号和棱边编号进行重排序，并输出对应的存储文件；

对所述三角形网格编号、顶点编号和棱边编号进行重排序具体为：

三角形网格编号：首先对离散到进程上不连续的三角形单元编号按照原编号大小顺序压缩重排为自1开始的连续编号，之后按照进程增序连续编号；

顶点编号：先将所有进程内的顶点编号，然后将公共顶点按顺序重新编号；

棱边编号：对进程内的棱边按原编号压缩重排后按照进程增序编号，然后对所有公共边压缩重排后编号，公共边编号的起始位置在所有进程内边编号之后；

步骤四、并行读取上述步骤三得到的存储文件，对三角形网格进行预设次数的细剖，并对新生成的棱边、顶点和三角形网格重排序，形成完整的细网格信息；

对所述新生成的棱边、顶点和三角形网格重排序，具体为：

三角形网格编号和顶点编号采用所有进程的全局编号；

顶点编号的原有顶点编号值不变，新产生的顶点按照从属边编号压缩后局部顺序重新编号。

2.如权利要求1所述的一种基于辅助树的多层快速多极子并行网格细剖方法，其特征在于，步骤四中，细剖的剖分尺度m与次数按如下方式确定：

m≤m₀×2ⁿ

其中，n是预设的细剖次数，m₀是最终需要达到的平均剖分边长，设置为十分之一波长。