HSS矩阵的快速算法及其在球谐展开中的应用研究

来源 :国防科学技术大学 | 被引量 : 0次 | 上传用户:yyk20071999
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
大规模稠密矩阵的数值线性代数是科学与工程计算的核心问题之一,此类问题广泛来源于机器学习的核密度估计、椭圆型边界值问题的边界元方法、数值天气预报的谱模式等诸多应用领域。通过近似分解或构造逼近基等方法对稠密矩阵中的低秩子块进行压缩,可有效地降低存储开销以及矩阵运算的计算复杂度。分层矩阵(Hierarchicalmatrices)可作为某些稠密矩阵的稀疏近似表示,根据矩阵划分方式和数据结构的差异,其可分为H-矩阵、H2-矩阵、HODLR矩阵、HSS矩阵和FMM矩阵这五个子类。目前,分层矩阵已被用于求解线性方程组、奇异值与特征值等问题。本文主要研究了HSS矩阵的基本理论与快速算法,以此出发改进了球谐展开的快速算法,并通过数值实验说明了算法的高效性与数值稳定性。本文的内容主要分为三部分:  一、以秩揭示QR分解、插值分解以及奇异值分解为主要对象,介绍了低秩矩阵构造近似分解的基本方法。本文还介绍了当前研究热点之一的随机低秩逼近算法。  二、详细地介绍了求解对称三对角矩阵特征分解的分而治之(Divide-and-Conquer)算法,以及用HSS矩阵改进后的HSSDC算法。本文还详细地介绍了构造HSS矩阵的随机算法以及HSS矩阵的快速乘法。  三、研究了MarkTygert提出的球谐展开的快速算法,并用HSSDC算法去加速连带勒让德变换。实验结果表明:对称三对角矩阵的特征分解在连带勒让德变换中的计算量比重很大,引入HSSDC算法能显著地加速连带勒让德变换。当截断波数?=18432且矩阵的规模n=18432时,连带勒让德变换的改进算法能达到4.21倍加速比,且算法的数值稳定性很好。
其他文献
本文研究如下多重调和方程组的正径向解的不存在性:{(-△)mu=a(|x|)vp,x∈RN,(1.1)(-△)mv=b(|x|)uq, x∈RN,其中m≥1,N>2m,p,q≥0且a,b是非负函数.我们主要是建立一些适当的条件使得此
该文主要研究泛函微分包含周期解的存在性.利用解的连续选择定理和渐近不动点定理,我们首先给出具非凸右端项无穷时滞泛函微分包含周期解存在的一般性结果.作为应用,我们证明
该文研究了国内外在病原体(虫卵)的自动检测中使用的方法和技术,提出了针对虫卵彩色显微图像的预处理方法和图像分割方法,引入了数学形态学的方法来处理分割后的虫卵图像,提
本文研究如下临界增长的重调和方程的解的存在性:  {△2u=|u|2*-2u+h,x∈Ω,u=▽u=0, x∈(a)Ω,其中Ω(∈) RN(N>4)是一个有界的光滑区域,h∈H-2(Ω),2*=2N/N-4是H2(RN)→L2*(RN)的临
该文我们证明了半局部右Kasch右单-内射环为半完全环.并证明了如下定理:R为右Kasch右A-内射环,当且仅当R为半完全右A-内射环且具有本质右基座,当且仅当R为右有限上生成的右A-
Hamilton-Jacobi(简称H-J)方程出现于最优控制、计算流体力学、计算机图形图像、微分几何、晶体生长、网格生成等许多领域.近年来,越来越多的数学家关注和研究H-J方程数学理
数学形态学是图像与信号分析、处理的非线性理论之一,其最大特点是它的几何直观性以及摒弃了传统的数值建模思想.该文首次将模糊逻辑和完备格相结合,提出一类具有广泛意义的
该文主要针对在开发设计"智能模拟三视图识图测试系统"CAI软件过程中遇到的在CAI软件中如何实现由平面三视图快速重构其三维模型的难点问题,提出了一种基于立方单元体(BOX)的
10月15~19日,深圳华意隆电气股份有限公司作为逆变焊接切割机行业的知名企业之一,在公司高层的大力支持下和海内外客户的热切关注与期待中盛装参加了第112届中国进出口商品交
该文共有五章组成.第一章中,较全面地概述了目前国内外全局优化的发展动态.第二章主要引入后二章构造算法及证明实现算法收敛性所需的数论中的主要结果,其中包括我们所作的一