基于HSS结构矩阵的快速算法与并行实现

来源 :中国科学院大学 | 被引量 : 0次 | 上传用户:zxqqqzxq
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
基于多层半可分(HSS)结构矩阵的快速算法可有效降低具有数值低秩属性的稠密线性方程组求解的复杂度。随机取样已经被证明是一种高效率的近似分解低秩矩阵的技术。采用随机取样和保结构秩显(SPRR)分解相结合的方法替代秩显QR(RRQR)分解可以快速构造HSS结构矩阵。该方法将压缩构造HSS结构矩阵转换成小矩阵计算,减少存储和通信开销,使构造HSS结构矩阵的时间复杂度进一步降低。  在分布式机群上采用ScaLAPACK的二维块循环分布方式存储各矩阵块,将HSS树和处理器网格进行映射。构造HSS结构矩阵的并行算法包括对矩阵的多层块压缩,数据交换和重分布,然后结合并行ULV分解和并行三角求解实现快速并行求解。文中分析了该并行算法的复杂度。同时文中以二维电磁散射问题为例,结果表明该算法不仅比直接LU分解快一个数量级,而且具有较好的可扩展性。  基于HSS结构的算法也为稀疏矩阵的结构分解奠定了一个良好的基础。在结构多波前算法中,快速计算Schur补能有效提高线性稀疏系统的求解。计算Schur补包括将波前阵构造成HSS形式,部分ULV分解和快速计算Schur补。这种基于HSS结构矩阵的计算Schur的算法保持了多波前算法的可扩展性,能用于大规模稀疏线性方程组的求解。同时数值实验也表明基于HSS结构计算Schur补的算法具有一定的可扩展性。
其他文献
随着通信网络和计算机网络的发展,网络融合是下一代网络发展的必然趋势。VoIP (Voice over IP)是在P网络上传送具有一定服务质量的语音的业务,是近几年内发展起来的一种新的
近年来,信息技术的普及和硬件技术的快速发展,为大数据产生与存储提供了先决条件。在商业上、科研机构、政府部门等都存储着大量的数据。而如何从这些大量的数据集中提取有用
随着互联网的不断发展,个人和小型团队开发者如雨后春笋般不断涌现出来,应用数目呈现高速增长趋势,传统自给自足的研发方式难以满足需求。平台即服务(PaaS)为应用程序提供基础运
焊接是机械工业化社会中一门重要的工艺,为保证焊接质量,需要进行缺陷检测。焊接缺陷若在钢管的内部,则需要采用无损检测的方式,较常用的就是X射线获取。早期主要是靠人工来
数字集成电路在制造过程中受到制造工艺、技术以及外围环境等因素的影响,可能会导致芯片中存在缺陷。在芯片上市之前,必须找出这些存在缺陷的芯片。但随着芯片规模及复杂度的不
该文提出了一个在Client/Server环境下,基于现有的OLTP系统来分析、设计实现OLAP系统的整体解决方案,意在花费最小的代价,通过对OLTP系统的功能进行有效的扩充来加速企业OLAP
随着移动通信和互联网的不断发展,电信运营商一方面需要满足人们日益增长的业务需求,另一方面需要抵抗互联网对于传统增值业务的冲击,这就要求电信运营商推出更加丰富多彩的
网络虚拟化技术作为解决当前网络僵化问题的新兴技术,获得了众多研究机构的广泛关注。本文在此背景下,重点关注网络虚拟化环境下的资源分配问题。资源分配是网络虚拟化环境下
近年来,VoIP应用得到了越来越广泛的应用,越来越多的企业和个人选择网络电话来进行实时通信。随着互联网中联网主机的数目日渐增加,很多主机都会处于NAT之后,网络地址转换技
Web服务的松耦合架构及其开发便捷性使得它的应用范围越来越广‘,越来越多的企业和个人热衷于开发新的Web服务并将其放在网络上,如何对这些Web服务进行识别成为Web服务研究中