ILUT和最小度算法在大型线性方程组求解中的应用研究

来源 :电子科技大学 | 被引量 : 0次 | 上传用户:Moon_____light
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
科学计算的一个重要课题就是求解大型的线性方程组,在实际应用中,稀疏矩阵占了很大部分,因此高效求解大型稀疏线性方程组就成为我们研究的主要方向之一。其中,不完全分解预条件子因为能够保持线性系统的稀疏性,降低存储的复杂性,减少计算量而受到重视。本文主要研究求解稀疏线性系统的性能算法,特别是构造有效的不完全分解预条件子。  首先介绍了经典的迭代法、预条件技术以及不完全分解;然后对不完全分解预条件方法ILUT(p,τ)进行改进得到算法MILUT(p,τ);最后结合最小度排序思想得到算法 MDILUTP(p,τ)。通过数值实验可以看出,算法对于稀疏矩阵的分解有着较好的效果。  对于稀疏矩阵而言,无论是对对称正定的矩阵分解构造预条件子,还是对一般矩阵分解构造预条件子,为了保证分解后矩阵的稀疏性质不会遭到破坏,通常我们都要设定非零元素控制参数。通过构造合理的非零元素控制参数,得到改进的不完全分解预条件子MILUT(p,τ),取得了一定的效果。针对非对称矩阵,研究基于ILUTP(p,τ)的预处理技术,结合最小度排序思想,在选主元过程中加入列非零元权值参数,使重排序后的矩阵在分解过程中减少填充元的产生,从而降低存储的复杂性,减少计算量,提高运行效率,同时确保矩阵的稀疏性在分解过程中不会遭到破坏。
其他文献
如何较好地处理非线性问题,已经成为各个领域内非常重视的课题.该文主要对在工程、物理等方面应用广泛的非线性积分方程在再生核空间W[a,b]中进行求解.首先对所需再生核空间W
近年来,流体运动所涉及到的方程模型种类极多。具有分数阶扩散的流体方程和不可压缩的MHD方程是流体力学方程研究的重要分支。目前,此类方程已被广泛应用到天气和气候预报,洋流
该文所考虑的一类具有波动算子的非线性Schrodinger方程具有多辛结构,从而我们引入正则动量把它写成多辛Hamiltonian方程组,并发现它有多辛守恒律、局部能量守恒律及局部动量
随着经济的不断建设和发展,人们越来越追求高尚而又有情操的精神卫华生活,音乐就是因此而诞生的产物,古往今来,音乐就一直有着其崇高的地位,现如今随着时代的变迁,教育也开始
该文的目的是试图建立p-adic分析的基本理论.主要内容如下:(1)本文首先建立了非负实数域R与p-adic数域QM(M一个零测度集)之间的可逆映射P.(2)借助于所建立的映射,构造了p-adi
全文共四章,内容如下:第一章介绍了该文工作的背景,概括了全文的研究意义、内容以及需要进一步研究和讨论的问题.第二章运用四元Heisenberg群的热核的具体公式研究相关的几何
纽结的Alexander多项式是纽结论中最重要的多项式不变量之一.该文要介绍它的推广,即扭的Alexander多项式,使得可以获得关于纽结基本群更多的信息.该文采用Wada关于有限表现群
对混沌控制和混沌同步的研究是当今现代混沌驾驭技术领域内的热点.该文是从统一的角度出发,将Lorenz型混沌系统分为确定和不确定系统,利用不同的方法和理论对其进行控制与同
随着机群并行计算技术的逐步成熟,机群已开始走出研究室,进入工程应用。但由于机群系统结构的特殊性,在进行基于机群系统的并行编程时,程序员并不能真正地将机群系统当作一台计算
局部有限性问题在从有限维代数的研究转向无限维的研究中有着很重要的意义.本文在对传统的КуроЩ问题的讨论上,提出新的问题,并对之进行了证明,主要讨论了以下问题:一、