灵活预处理GPBiCG和并行BiCG方法

来源 :中国工程物理研究院 | 被引量 : 0次 | 上传用户:hou189
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在科学工程计算中,许多大规模科学计算问题的数值模拟最终可以归结为大型稀疏线性代数方程组的求解。而代数方程组的求解时间往往在整体数值模拟时间中占有非常大的比重,以致成为整体数值模拟的瓶颈。因此设计高效的代数解法器是求解这类问题的关键所在,同时也是设计相关高性能软件的基石。目前,在给定Krylov子空间迭代算法前提下,有两种途径可以降低大型稀疏代数方程组的求解时间:一种是并行计算,另一种是使用预处理技术。  本文主要研究了两项内容:一是利用FLAME框架重新推导BiCG方法,重排内积顺序,降低全局通信次数,提高方法的并行可扩展性;二是对于非对称线性方程组,给出了灵活预处理GPBiCG方法。  首先,介绍如何利用FLAME框架形式地推导Krylov子空间法。然后考虑BiCG方法,利用算法中向量序列的性质,即剩余向量的正交性和搜索方向的双共轭性,推导出几种不同但是数学上等价的BiCG算法,选出适合于在分布式存储处理机上并行计算的那些算法,进行了理论分析。并且用并行试验来表明新算法的确降低了全局通信时间,提高了算法的并行可扩展性。  其次,我们给出可变预处理形式的GPBiCG方法,在方法的每一步中可使用不同的预处理子。特别地,可变预处理子的灵活性是可用任何一种迭代法得到。例如,标准的GPBiCG方法自身可以作为预处理子,其他的Krylov子空间法或是分裂迭代法也可以。对于可变预处理形式的GPBiCG方法,我们还进行了一些数值试验,这些算例表明了可变预处理迭代法的收敛性和健壮性。
其他文献
线性码是一类非常重要的分组码,是讨论各种码的基础。LDPC 码是一大类可以用稀疏校验矩阵定义的线性纠错码,因其校验矩阵为稀疏的,故称其为低密度的,LDPC 码是目前编码理论中的热
随着我国社会主义市场经济体制的建立和完善,市场体系的日益健全,上市公司日渐增多,投资风险也不断增大。投资者要准确及时地掌握上市公司的经营状况和财务状况,做出合理的投资决