岩土工程亿级单元有限元模型可扩展并行计算

讨论了亿万单元有限元模型的可扩展并行计算方法。从软件和硬件两个方面提出了前处理、并行计算方法、程序算法、后处理实现等核心问题的解决方案。采用网格加密方法生成一亿单元的有限元模型,利用对偶原始有限元撕裂内联法（FETI-DP）求解系统方程。基于图论理论建立了子区域间的通讯拓扑关系,实现了子区域间点对点通讯,避免速度慢、通信量大的全局通讯。在自主开发程序基础上,增加相应模块,采用面向对象编程技术和MPI消息传递库开发程序。对一个一亿多单元的工程实例运用5 000核并行计算,得到了超线性加速比。计算结果在专用图形工作站上进行后处理,显示和交互操作速度良好。研究在两方面实现了突破：一是将模型规模提高到了一亿多单元;二是同时调用了5 000个计算核来并行运算,并得到了很高效率。高分辨率有限元并行模拟研究成果可为岩土工程中结构特别复杂、计算区域特别大、地质情况复杂等模拟提供很好的技术方法和实现手段。     

大地电磁三维数据空间反演并行算法研究

目前大地电磁三维反演实际应用的主要问题是计算效率低.在对大地电磁三维数据空间反演算法进行深入分析的基础上,本文提出了基于频点和矩阵划分的大粒度并行反演方案和具体实现步骤,并在曙光TC5000A高性能计算平台上实现了基于MPI的大地电磁三维数据空间反演并行算法.该算法实现了包括三维正演、灵敏度矩阵、叉积矩阵以及模型改正量的并行执行,不仅计算效率高,而且每个节点机上灵敏度矩阵的存储空间只需原来微机上的2/N(N是参加并行计算的节点机个数),大大地减少了内存开销.通过两个理论模型合成的数据对实现的三维数据空间反演并行算法进行试算,对比分析了多个节点机下程序的执行效率.测试结果表明,所实现的三维数据空间反演并行算法是可行的、高效的,与单机相比,不仅可以提高运行速度,缩短计算时间,而且还可以扩大计算规模,极大地推动了大地电磁三维反演的实用化.     

三维地震动场数值模拟并行计算系统

并行计算系统是进行三维复杂场地地震动场数值模拟的基础。本文作者利用10台微机构建了一套并行计算集群(常称为Beowulf系统),节点微机上运行Windows2000操作系统,采用MPI标准的MPICH1．2．4版本作为并行计算的支撑环境,并以100M bps高速交换式以太网作为互连网络。利用编制的并行计算程序对该并行计算系统性能进行了测试。结果表明：该系统能够达到非常高的并行加速比和并行效率。该系统将用于三维地震动场数值模拟研究。     

基于集群MPI的图层级多边形并行合并算法

在集群环境下,基于MPI并行编程模型和OGC简单要素规范进行并行多边形合并时,需要处理叠加图层间要素的“多对多”映射关系,由于空间上相邻的多边形在要素序列上并不一定连续,导致无法按要素序列为子节点分配任务,给并行任务映射带来了困难。本文以集群环境下的并行多边形合并算法为研究对象,通过比较叠加分析中两种多边形映射关系对算法并行化带来的影响,基于R树空间索引、MySQL精确空间查询,以及MPI通信机制,提出了6种不同的并行任务映射策略;通过实验分析和比较了6种策略的优劣。结果显示：基于R树预筛选的直接合并策略,在各算法中具有最高的串行计算效率和优秀的并行性能表现。虽然MySQL精确空间查询的预筛选过程较为耗时,但可有效地过滤掉不真正相交的多边形,从而提高合并操作的效率。因此,在集群MPI环境下,基于R树和MySQL精确空间查询的预筛选策略是解决并行任务映射难题,实现图层级多边形并行合并算法的有效途径。     

多接收阵合成孔径声纳后向投影成像快速算法

提出了一种内存共享环境下的多子阵合成孔径声纳后向投影快速算法。首先根据后向投影算法的特点设计二维动态数据结构,保存同距离向上的时延和补偿相位信息,省去时延的重复计算;然后在共享内存环境下,将远距离和近距离聚焦运算相结合,实现计算负载平衡;最后通过实测数据的成像试验验证了本文算法的正确性和高效性。     

基于异构集群计算的地统计面插值并行算法

地统计面插值算法在空间统计分析中有广泛应用,其目的是通过一组面要素的某已知属性值估算另一组面要素的属性值。地统计面插值算法多是基于克里金（Kriging）插值及其衍生算法。克里金插值算法考虑属性在空间位置上的变异性,需计算要素之间的协方差,是典型的计算密集型算法。本文分析了基于克里金插值的地统计算法计算过程,该算法中面要素间协方差计算相互独立,可作为并行计算单元划分。另外,面要素间协方差计算可使用快速傅里叶变换（FFT）快速计算,而FFT是一种非常适合并行处理的计算密集型算法。本文根据算法特征设计了基于异构集群计算的并行算法,并使用MPI+CUDA实现了该算法。实验结果表明,本文实现的算法比使用MPI实现的CPU集群的算法有更好的性能,具备良好的可扩展性,并且随着插值精度提高表现出更好的性能。     

GRAPES全球格点模式的并行计算负载平衡策略

随着高性能计算机技术的发展和应用,并行计算已成为保证数值天气预报模式业务运行时效的关键技术之一.目前高性能计算机计算能力已达到每秒千万亿次浮点计算.系统中处理器数量也早已达十万甚至更多,如此巨大的计算资源对应用软件系统的设计也提出了挑战.数值天气预报软件系统要充分利用高性能计算机提供的计算资源,必须依靠并行计算方法,这包括适合计算问题的可扩展并行算法的设计、合适的数据分配方案以及良好的任务负载平衡方案.作为中国新一代数值天气预报格点模式,GRAPES(Global and Regional Assimilation and PrEdiction System)设计的最终目标是一个科研/业务通用,区域/全球通用模式.作为一个格点模式,GRAPES的并行计算具有与欧洲中期数值顶报研究中心谱模式并行计算不同的特点,GRAPES的并行计算采用了经典的水平网格数据划分.但对于全球的GRAPES模式,由于采用拉格朗日差分方案,模式极地及附近区域格点与格点之间距离的减小.使得模式并行计算在采用简单的经纬网格划分方式实现时,必须考虑极地区域并行计算跨越多个处理器时导致的频繁通讯解决途径.本研究提出了利用消息传递组通讯实现全球格点模式并行计算的一种方法,其核心思想是将极点附近一定区域内的处理器按纬向划归不同的处理器组.文中还给出了该实现方法的任务分配算法,提出了改进的任务分配负载平衡方案.在中国气象局高性能计算机IBM-cluster1600上的测试表明,算法具有较好的可扩展性,其负载平衡方案改善了计算的绝对墙钟时间,使并行计算效率提高10%以上.模式的准业务运行结果表明计算墙钟时间基本可以满足数值预报业务的实时性要求.     

基于OpenMP的Li-Openshaw算法并行计算方法

随着网络地图、移动地图的发展,对地图自动综合效率的要求不断提高;并行计算技术已经在遥感图像处理、计算几何等方面显示出优势。本文尝试将地图综合算法与并行计算方法相结合,探讨地图综合算法效率提高的方法;并通过对Li-Openshaw线要素简化算法的分析,探讨了在多核处理器环境下该算法的并行计算策略,基于OpenMP实现了Li-Openshaw算法的并行计算,最后通过实验验证了基于OpenMP的Li-Openshaw线要素化简算法在效率上的提高。     

弹性波传播数值模拟的区域分裂法

本文针对弹性波波场传播的数学模型,提出了一种基于有限元和有限差分耦合的区域分裂方法．它有灵活的网格剖分,克服了单纯用差分方法对区域的依赖性,可以适用于曲边地表;达到同样精度所需的计算量比有限元方法小;并易于实现并行计算．数值实验表明了算法的有效性．     

OpenMP在MM5中尺度模式中的应用试验

简要介绍了一种新的并行编程标准——OpenMP通过与其他几种并行编程标准的比较,了解其优越与不足之处。详细介绍了OpenMP在PSU/NCAR的中尺度模式MM5中的运用。分别用OpenMP并行方式与MPI并行方式在SGI Origin 2000上进行了并行效能测试。测试结果表明：OpenMP是一种简单、高效的并行编程标准,非常适合于气象科研与业务领域应用。