适合分布式并行计算的PBiCOR方法

技术领域

本发明涉及并行机通讯技术领域，特别是适合分布式并行计算的PBiCOR方法。

背景技术

随着并行计算机的发展和应用，复杂物理系统的高分辨率数值模拟已被提上议事日程，并行计算机按存储方式划分，可分为共享内存和分布式内存两种基本的并行计算机存储方式，共享内存的并行计算机，各个处理单元通过对共享内存的访问来交换信息、协调各处理器对并行任务的处理，对这种共享内存的编程，实现起来相对简单，但共享内存往往成为性能特别是扩展性的重要瓶颈，分布式内存的并行计算机，各个处理单元都拥有自己独立的局部存储器，由于不存在公共可用的存储单元，因此各个处理器之间通过消息传递来交换信息、协调和控制各个处理器的执行，不难看出，通讯对分布式内存并行计算机的性能有重要的影响，复杂的消息传递语句的编写成为在分布式并行计算机上进行并行程序设计的难点所在，但是，对于分布式并行计算机，由于它有很好的扩展性和很高的性能，它的应用非常广泛。

对于分布式存储并行机，由于矩阵和向量被分布存储在各处理器上，因此，即使矩阵运算通过并行运算可有效地实现，全局通讯，即内积计算所引起的所有处理器间的通讯，我们仍不能有效避免，Yang等分别提出了并行化的CGS和Bi-CGSTAB方法，Sturler等提出了如何降低GMRES和CG算法中全局通讯影响的方法，Gu和Zuo等分别给出了并行化的BICGSTAB(2)、GPBiCG和IICGS方法，Collignon等给出了并行IDR(s)方法；同时，基于区域分解，Gu，Liu和Mo给出了一种不需整体内积计算的CG类方法，即多搜索方向共轭梯度(MSD-CG)方法，该方法将CG方法中的内积计算用小线性方程组来代替，从而完全消除了全局通讯。

发明内容

首先，在超融合的服务器上，存储跟计算节点融合外，我们采用sas3为存储交换网络，这跟传统的分布式不太一样，传统的分布式它的架构是每个节点都会有自带的硬盘，或者是所谓的固态盘，也就是数据，在这个节点上，数据跟其它节点之间是独立的。发明内容发明目的本实用新型的目的是为了克服现有技术中的不足，提供一种中心和四周双层固定、曝气均匀、抗撕裂、使用寿命长、曝气面积大、利用率高的大面积膜片式微孔曝气器。同时，我们在存储网络基础上，又跟传统分布式存储采取不一样的技术与方式，传统分布式存储与计算因为透过网络交换，它的数据独立在各个节点上，所以在数据交换的过程中是采取副本的方式来保护数据，我们知道副本的方式保护数据，多一个副本对磁盘阵列空间来说存储空间就直接减一半。

一种用于检测管桩倾斜程度的装置，该管桩倾斜检测装置包括主控单元、命令输入单元、传感器单元、信号调理单元、a/d转换单元(c8051f350)、存储单元、显示单元、通信接口，其特征在于命令输入单元与主控单元相连，传感器单元依次经过信号调理单元及 a/d转换单元后与主控单元相连，存储单元与主控单元直接相连，显示单元与主控单元相连，通信接口与主控单元相连，存储单元包括sd卡以及flash存储器。专利摘要本实用新型公开了一种用于检测管桩倾斜程度的装置，该管桩倾斜装置包括主控单元、命令输入单元、传感器单元、信号调理单元、a/d转换单元、存储单元、显示单元、通信接口，命令输入单元与主控单元相连，传感器单元与信号调理单元相连，传感器单元中的电子罗盘直接与主控单元相连，传感器单元中的加速度计与信号调理单元相连，电子罗盘直接与主控单元相连，温度传感器直接与主控单元相连，传感器单元依次经过信号调理单元及a/d转换单元后再与主控单元相连，存储单元包括sd卡以及flash存储器，均与主控单元直接相连，显示单元与主控单元相连，通信接口与主控单元相连。a/d转换单元110与信号调理单元108及主控单元102相连，存储单元112包括存储测量数据的sd卡128以及存储汉字字库及临时缓存用的flash存储器130 (flash)，均与主控单元102直接相连。

假定采用SPMD模型，即所有的处理单元执行相同的程序，且一个处理单元需要远程处理单元的数据时，必须通过消息传递完成；

P为处理器台数，N为稀疏矩阵A的维数，nz为矩阵A每行(列)非零元素的平均个数，tfl为一个浮点操作的时间，ts为一个消息启动时间，tw为一个字在两台处理器间传递所需的时间，一个内积或一个向量校正的计算时间为2tflN/P，其中N/P为分配给每个处理器的未知量个数，(2nz-1)tflN/P为稀疏矩阵向量乘积的计算时间；

具有P台处理器的矩阵(或网状)分布式并行式，并假设通讯是二叉树方式进行的，则一个内积的整体聚集和广播的时间为2logP(ts+tw)，而同时聚集和广播K个内积的时间为2logP(ts+ktw)，稀疏矩阵已被映射到各处理器，使得矩阵向量乘积只需相邻处理器间的局部通讯，矩阵向量乘积所需的通讯是交换边界数据，即从其他处理器接收边界数据和向其他处理器发送边界数据，假设每个处理器必须接收和发送nm个信息，每个处理器上有nb边界数据单元个数，因此，每个处理器必须进行通讯的总字长为2(2nb+nm)；

由于不需通讯，一个向量校正的时间为：

tvec_upd＝2tflN/P(2)

只需一次全局通讯的K个内积的时间为：

tinn_prod(k)＝2ktflN/P+2logP(ts+ktw)(3)

只需局部通讯的一个矩阵向量乘积的时间为：

tmat_vec＝(2nz-1)tflN/P+2nmts+2(2nb+nm)tw(4)

一个矩阵向量乘积的通讯时间为2nmts+2(2nb+nm)tw。

本发明实现全局通讯次数减少，使得PBiCOR方法的并行通讯性能得到了有效提升，具有更好的并行性和扩展性，提高了通讯性能改进比率，势必为油藏模拟、天气预报、计算流体力学、计算电磁学、天体物理和惯性约束聚变(ICF)二维数值模拟等领域提供了理论基础。

附图说明

图1为已知文献中的双共轭A-正交的BiCOR方法。

图2为本发明的适合分布式并行计算的PBiCOR方法分布式存储并行机。

图3为BiCOR方法与PBiCOR方法每次迭代的计算量及所需全局通讯次数的比较情况分布式存储并行机。

具体实施方式

以下结合附图，对本发明的具体实施方式作进一步详细说明。

现结合图1至图3所示，本发明适合分布式并行计算的PBiCOR方法，假设并行机具有P台处理机，每台处理机分别具有局部存储模块和相应的处理单元，处理单元之间通过互联网相连；

假定采用SPMD模型，即所有的处理单元执行相同的程序，且一个处理单元需要远程处理单元的数据时，必须通过消息传递完成；

P为处理器台数，N为稀疏矩阵A的维数，nz为矩阵A每行(列)非零元素的平均个数，tfl为一个浮点操作的时间，ts为一个消息启动时间，tw为一个字在两台处理器间传递所需的时间，一个内积或一个向量校正的计算时间为2tflN/P，其中N/P为分配给每个处理器的未知量个数，(2nz-1)tflN/P为稀疏矩阵向量乘积的计算时间；

具有P台处理器的矩阵(或网状)分布式并行式，并假设通讯是二叉树方式进行的，则一个内积的整体聚集和广播的时间为2logP(ts+tw)，而同时聚集和广播K个内积的时间为2logP(ts+ktw)，稀疏矩阵已被映射到各处理器，使得矩阵向量乘积只需相邻处理器间的局部通讯，矩阵向量乘积所需的通讯是交换边界数据，即从其他处理器接收边界数据和向其他处理器发送边界数据，假设每个处理器必须接收和发送nm个信息，每个处理器上有nb边界数据单元个数，因此，每个处理器必须进行通讯的总字长为2(2nb+nm)；

由于不需通讯，一个向量校正的时间为：

tvec_upd＝2tflN/P(2)

只需一次全局通讯的K个内积的时间为：

tinn_prod(k)＝2ktflN/P+2logP(ts+ktw)(3)

只需局部通讯的一个矩阵向量乘积的时间为：

tmat_vec＝(2nz-1)tflN/P+2nmts+2(2nb+nm)tw(4)

一个矩阵向量乘积的通讯时间为2nmts+2(2nb+nm)tw。

考虑迭代求解大型稀疏复对称线性代数方程组

Ax＝b,(1)

其中A∈CN×N，A≠A-T,A＝AT，x,b∈CN，记x1和r1＝b-Ax1分别为初始猜测和初始残差，在分布式存储并行机上，假设系数矩阵和相应的向量，已经按行(列)被映射到各个处理器上，且负载平衡，由此可知，矩阵向量乘积只需局部通讯，向量内积计算需要全局通讯，而全局通讯是影响Krylov迭代方法并行计算效率的主要瓶颈。

由图1可知，第4)和第8)步骤中需要内积计算，而其他计算步骤又和内积计算密切相关，且从BiCOR算法中可以看出，第4)和第8)步的内积计算具有紧密的数据关联性，因此，BiCOR算法的每一步迭代中就需要两次全局通讯，这对其并行执行效率具有非常严重的影响。

本文来自电脑杂谈，转载请注明本文网址：
https://www.pc-fly.com/a/jisuanjixue/article-101886-1.html