超级计算机上矩阵乘的并行计算与实现

数值预报系统中经常要用到矩阵乘运算。在分布式超级计算机 (如IBM-SP) 上, 矩阵乘的并行计算需要较多的数据移动, 有效的数据传输对矩阵乘的实现至关重要。该文讨论了两种矩阵乘的并行算法, 一种是基于矩阵的列-行划分方式, 一种是基于矩阵的网格划分方式。在IBM-SP计算机上的实验结果表明, 网格划分的矩阵乘并行算法通讯开销更小, 并行效率更高, 其并行加速比较列-行并行算法改善约10 %。     

一种在异构系统中实现负载平衡的方法

提出了在异构系统实现负载平衡的区域分解算法和实现负载平衡的计算方法，利用它的负反馈性质解决了异构系统处理机计算速度测量误差造成的负载测量不准问题，并对处理机速度变化，速度测量误差、处理机数量、网格点计算量的分布等因素的影响进行了计算，结果表明本方法具有很强的平衡负载能力和较强的适应性；根据计算结果提出了解决模式网格点计算量不易测量问题的解决方案，并用扩散方程和模拟物理过程进行试验，试验表明这种方法是可行的，平衡负载的效果十分显著。     

基于GPU的GeoEye-1影像正射纠正

在介绍和分析GPU并行编程方法的基础上,给出了在GPU上进行Geo Eye-1影像正射纠正的算法过程。利用NVIDIA Quadro 600显卡对31 985像素×15 251像素的3波段Geo Eye-1影像进行正射纠正,其总体速度比没有加速的算法提升40%以上。试验结果表明,采用GPU并行处理技术可以提升遥感影像处理效率。     

共享内存环境下的分块最小不连续相位解缠算法

提出了一种共享内存环境下的分块最小不连续相位解缠方法。首先,利用多核协同计算整块相位质量图;然后,将干涉图和相位质量图分割为规则小块,采用质量引导与最小不连续相合成的方法同时对子块干涉图进行相位解缠;最后,主线程在不同解缠相位块以及位于边界的低质量区域上再次进行最小不连续优化,获取最终合并后的解缠结果。仿真和真实干涉相位图的分块解缠实验结果表明,所提算法在保持解缠精度的同时,加速比分别达到3.98和2.26。     

基于CUDA的高阶旋转交错网格有限差分法的弹性波正演模拟

应用高阶旋转交错网格的有限差分法,求解了基于一阶的速度-应力弹性波方程,推导了时间二阶精度空间2M阶精度的高阶有限差分离散格式.为了提高弹性波正演模拟效率,论文探讨了常规的CPU串行计算、基于GPU平台的节点并行计算、基于GPU平台的分量、节点同时并行的3种正演模拟计算模式,给出了不同算法的流程示意图.通过对比3种算法的耗时及后两种算法的加速比,可清晰判定,基于GPU平台的节点并行算法较CPU平台的串行算法有明显加速,最大加速比可达到60倍,而基于分量、节点同时并行算法的最大加速比最大可到180左右,加速效果更为明显.最后应用Marmousi模型的15000次时间步迭代的正演计算,基于分量、节点GPU并行算法大约是节点并行计算耗时的1/3,说明了该并行设计能有效提高弹性波正演模拟的效率,节约计算耗时.     

集群环境下的复合最小不连续相位解缠算法

提出了一种集群环境下的复合最小不连续相位解缠算法。首先主线程根据计算资源数将原始缠绕相位分为规则小块,并将未解缠相位块发送至空闲计算节点进行解缠。单块缠绕相位图解缠时,先计算相位质量图,并将缠绕相位分为高低质量区域,然后采用质量引导与最小不连续相结合的复合相位解缠策略进行解缠,最后将解缠结果和区域分割结果发送回主线程。完成所有分块缠绕相位解缠后,主线程在不同解缠相位块边界及其与边界相邻的低质量区域进行最小不连续优化来获取最终的解缠相位。通过集群环境下的并行相位解缠试验,验证了所提算法的正确性和高效性。     

网格环境下缓冲区分析的并行计算

提出了基于网格环境下缓冲区分析的并行计算方法,用以提高缓冲区分析的效率。考虑到大多数地理信息系统中缓冲区分析数据的特点,提出了基于图层和地理空间区域的任务分解方法。利用已经构建的网格计算环境,实现了武汉城市地图中一条道路缓冲区分析的图层式并行计算,从实验结果的分析,可以看出网格环境下缓冲区分析的并行计算在效率方面具有明显的优势。     

大型GNSS观测网的并行计算

CGCS2000精化及其全球拓展需要采用最新的模型,处理长期积累的大型GNSS观测网数据,大型GNSS网联合、快速和协同解算是空间基准精化、维持与服务的重要技术方向.在大数据技术背景下,以并行计算、云计算为代表的高性能计算技术逐渐成为大规模数据处理的首选方法.针对海量、多源、异构GNSS数据在解算处理与平差分析等方面面...     

地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.