非对称稀疏矩阵特征问题的并行迭代解法

来源 :中国科学院研究生院 中国科学院大学 | 被引量 : 0次 | 上传用户:hbzhwyf
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
矩阵特征值问题是科学计算的一个重要组成部分,其研究可以追溯到一个半世纪之前.许多应用都会归结为矩阵特征值问题,例如在材料科学和化学中离散偏微分方程后就需要求解矩阵特征值问题,而这些问题的矩阵规模很容易超过百万阶.然而,这类问题的矩阵通常是稀疏的,这使得我们可以用很少的运算量计算矩阵一向量乘积.而且这类问题只计算矩阵的一部分特征值.传统的正交变换方法将会产生稠密矩阵,因此受到问题规模的限制,迭代法被广泛用于计算大规模稀疏矩阵特征值问题.因为矩阵特征值问题经常是计算瓶颈,所以研究快速收敛的迭代方法及其高效的并行实现具有重要意义。   重新启动的Arnoldi算法是求解大规模稀疏矩阵特征值问题的一类重要方法.本文用向量与子空间的夹角分析了Krylov子空间的构造特点,分析了各种带有重新正交Arnoldi过程的并行实现的计算和通信开销.在此基础上,本文提出几种带有选择重新正交的Arnoldi过程的并行实现.本文所给方法减少了不必要的判断是否重新正交的次数,从而使其具有较少的计算和通信开销。   隐式重新启动的Arnoldi算法在求解许多矩阵特征问题上是一个成功的算法.然而对于某些困难问题,如果子空间维数较小,那么隐式重新启动的Arnoldi算法收敛很慢甚至收敛不到我们想要的特征信息.从并行计算的角度看,每个进程只存储矩阵和向量的一部分,并且只承担总的运算量的一部分,因此并行计算适于通过增大子空间维数来提高Arnoldi算法收敛速度.本文提出用流水线模型并行计算隐式重新启动过程中的一系列带位移的QR分解,使得子空间维数增大后,隐式重新启动的Arnoldi算法的并行实现仍然具有较大的加速比。   精化策略具有改善Arnoldi算法收敛性的特点,其核心思想是用精化Ritz向量代替Ritz向量作为特征向量的近似.求解每个精化Ritz向量需要计算与投影矩阵相关的奇异值分解.本文根据投影矩阵具有上Hessenberg矩阵的特点,提出用Givens变换部分代替Householder变换的奇异值分解算法,从而减少了求解精化Ritz向量的运算量.此外,本文利用各个奇异值分解之间不存在数据相关性的特点,提出二级并行方法计算精化Ritz向量,降低了奇异值分解的并行实现的弱扩展性对算法总的并行性能的影响。   分块Arnoldi算法具有浮点运算效率高,对重特征值或特征值簇问题快速收敛的优点.与Arnoldi算法类似,分块Arnoldi算法要求基向量在数值上正交.本文提出一种带有重新正交的分块Gram-Schmidt正交化过程,并将之应用到分块Arnoldi算法.误差分析的结果表明本文所给算法能够保证生成的基向量的正交性达到机器精度。
其他文献
在C2C电子商务中,信用评价至关重要,直接影响交易成功与否。目前的C2C电子商务信用评价系统,评分等级的设置相对简单,不能很好的反映真实的交易状况的问题。从而产生信用诋毁、信
随着经济的发展与人们生活方式的改变,中国已成为世界上移动终端用户最多的国家。能够支持移动接收的多媒体广播服务因庞大的移动终端使用人群而显得日渐必要。国际上现有一
当前,国内水环境污染十分严重,尤其是江河流域普遍遭到污染。水利部对全国700余条共约10万km长的河流开展的水资源质量评价结果表明:水质污染严重而不能用于灌溉(即劣于Ⅴ类)
学位
今天Web应用程序的界面不再是通过Web服务器中的模板生成,而是通过浏览器中的JavaScript生成。用Ajax技术构建Web应用程序,使Web应用程序的架构产生了一次重大变革。但通过Java
视频图像中的运动目标跟踪作为计算机视觉领域的核心研究课题之一,经过近50多年的发展得到了广泛而深入的研究。它融合了图像处理、模式识别、人工智能、自动控制等相关领域的
学位
动态心电图(DCG,Dynamic Electrocardiography)是临床上分析诊断心血管疾病的重要手段。如何从心电信号中有效地提取各种特征并进行分类识别处理,辅助医务人员进行各种心血管疾
随着计算机软硬件技术的飞速进步与发展,分布式虚拟现实系统受到各行各业越来越多的关注,在军事、医学、建筑、娱乐、教育等领域都显示出巨大的经济和社会效益。本文探讨了分布
龙芯系列处理器是中国科学院计算技术研究所自主研发的,兼容MIPS指令集的高性能通用处理器。目前已经广泛应用于各个领域,包括高性能计算、桌面和网络安全等。为了充分发挥龙芯
自然计算(Natural Computation)是表示由自然启发的计算的一般性术语,其研究内容一般包括人工神经网络,遗传算法,免疫算法,蚁群算法和粒子群算法等。由于多数的自然计算模型
近几年随着互联网以及电子商务的飞速发展,互联网中的信息呈现出爆炸性的增长,用户无法从过量无用信息中挖掘出自己需要的物品或者信息,在这种情况下,个性化推荐系统应运而生