分布式存储环境下矩阵广义特征值问题的并行计算

来源 :中国人民解放军国防科学技术大学 | 被引量 : 0次 | 上传用户:zmatch
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
矩阵广义特征值问题是当前迅速发展的计算机科学和数值代数中的一个非常活跃的研究课题。它在很多应用中扮演非常重要的角色,从数学角度来看,矩阵特征值问题的应用大多来自数学物理方程、差分方程、Markov过程等。目的是为了计算固体、流体、电磁、微观粒子、系统控制等重大问题。在实际的科学研究与工程应用中,比如在建筑工程、航空航天研究、生物科学、计算物理以及石油勘探中,都要涉及到大规模矩阵广义特征值问题的计算。 大规模并行处理系统(MPP)和PC机群为并行求解矩阵广义特征值问题提供了分布式存储环境。本文从理论和实验两方面深入研究了分布式环境下实矩阵广义特征值问题Ax=λBx的并行计算,提出了一些新的算法。本文的主要研究成果概括如下: (1)关于实对称定三对角矩阵广义特征值问题 ⅰ) 提出了一种基于等规模矩阵划分策略的分治算法,子问题的规模之和等于原问题的规模。该算法利用割线法迭代计算广义特征值,具有计算量少,计算速度快的优点,并且算法具有内在并行性。理论分析和数值实验均证明该算法比K.Y.Li算法快,在矩阵规模较大时,计算时间约可以减少10%-20%。当n→∞时,加速比S_p→p,在分布式环境下可以获得很好的线性加速比和稳定的效率。 ⅱ) 提出了一种具有全新矩阵划分策略的分治算法。在这个划分策略中,子问题的规模之和等于原问题的规模减1,文中给出了特征值分割定理及其证明。该算法所提供的特征区间要优于第一种算法,因而在采用相同的迭代方式时,算法具有计算量少、计算时间短的优点。在矩阵规模较大时,计算时间约可以减少30%甚至更多。同时,算法也具有良好的内在并行性,当n→∞时,加速比S_p→p,在分布式环境下可以获得很好的线性加速比和稳定的效率。 (2)关于实对称带状矩阵广义特征值问题 ⅰ) 提出了一种结合多分法的并行分治算法,给出了特征值分割定理及其证明。在该算法中,子问题的规模之和等于原问题的规模。算法利用二分法结合广义Rayleigh商迭代,计算包含在特征区间内的每个特征对。在计算全部特征值和特征向量的情况下,算法的计算复杂性为O(n~2r~2),当r<<n时,优于坛一一一~一续攀遨粼望廷二生掣夔 LApACK的计算复杂性O(n’)。在并行实现时,算法结合并行性好的多分法, 在分布式环境下获得了很高的并行加速比和效率。 ii)提出了另一种结合多分法的并行分治算法,给出了特征值分割定理及其证明。 在该算法中,子问题的规模之和等于原问题的规模减;,其中;为半带宽。在 计算全部特征值和特征向量的情况下,算法的计算复杂性为。(。’:’),并且该 算法的计算时间绝大多数情况都少于前一种分治算法。结合多分法,算法具 有很好的并行性,在分布式环境下,获得了很高的并行加速比和效率。 iii)提出了一种基于分治策略的并行同伦连续方法。通过限制步长、抛弃特征值 束的方法,提高算法效率。在矩阵规模较大时、算法的计算时间明显比 LAPACK少。在算法的并行实现时,充分利用计算与通信的重叠,在分布式 环境下获得了很高的并行加速比和稳定的效率。 (3)提出了一种计算实对称定矩阵广义特征值问题的并行分治算法。该算法是在计 算矩阵B的谱分解和MDR约化的基础上,将问题转化为对称三对角一正定对角 矩阵广义特征值问题,然后利用分治算法结合割线法迭代进行计算。根据扰动 理论,给出了每个特征值的包含区间。数值实验表明,在分布式环境下,该算 法可以获得很高的并行加速比和效率。 (4)提出了一种计算非对称矩阵广义特征值问题的分治算法。该算法在分治策略的 基础上,利用Laguerre迭代计算广义特征值,并采用各种方法尽可能减少特征 值的遗漏。采用收缩方法计算遗漏的特征值。数值实验表明,该算法计算速度 快,并行效果好。 (5)在分布式环境下实现了本文提出的所有算法,并形成了一个符合规范的软件包。
其他文献
随着我国事业单位会计制度的发展,初步完善了事业单位财务处理中的弊端和漏洞,创新了财务处理工作方法。2013年,财政局发布了新的《事业单位会计制度》,进一步规范了事业单位
目的:探讨甲状腺腔镜下行部分切除术的临床效果。方法:选取收治的甲状腺腔镜下部分切除术患者87例,对其临床资料进行回顾性分析。结果:85例成功完成甲状腺手术,术中出血量50~200
运用虚拟现实技术模拟海上交通与控制,可以提供一个用户与计算机系统自然交互的三维人机界面。为了进一步提高用户与地理信息系统交互的深度与广度,突破传统图形用户界面,1995年
目的:探讨甲氧西林耐药金黄色葡萄球菌(MRSA)耐药性。方法:从各临床标本中分离出的不重复MRSA菌株,将其分为社区获得组和医院获得组,并分析其耐药性原因。结果:医院获得组耐药性明
计算机辅助工艺设计(CAPP)作为计算机集成制造系统的的关键技术,是连接计算机辅助设计和计算机辅助制造的桥梁,是建立工程产品生产的全局信息模型的基础和关键。现代企业制造对工
随着互联网技术的发展,XML已经成为数据表达和交换的新标准。各行各业都积极采用XML格式表示其数据和信息,借助于XML,Web正逐步转化为一个巨大的数据库。由于XML信息量的剧增
目的:研究半完壁式乳突病变切除鼓室成形术治疗中耳后天原发性胆脂瘤的疗效。方法:行半完壁式乳突病变切除鼓室成形术治疗中耳后天原发性胆脂瘤21例,同期行开放式乳突切开鼓室
<正> 张国榉是我国著名的凿岩工程专家。30多年以来,他潜心于凿岩钎具的研究。由他研制的各式新型钎具,在大瑶山隧道、四川沪沽铁矿、京九铁路、三峡水电站,以及全国矿山和国
异构数据库由一组互连的、自治的成员数据库组成。这些成员相互通信,以便交换信息和回答查询。一个成员数据库中的对象可能与其他成员数据库中的对象很不相同,使得很难将他们
在分析了生物质燃烧过程的基础上,设计了一种新型户用高效生物质直燃炉。该炉点火方便、操作简单、节能安全,一次风、二次风配风系统分开,对燃料无苛刻要求,配有星型给料机可