复杂特厚煤层综放工作面煤柱应力分布规律研究-煤炭科学技术.PDF

科研信息化技术与应用 2 0 1 4 , 5 ( 3 ) : 4 5 – 4 9 TECHNOLOGY ·技术 /  · LU 一种特殊稀疏矩阵的 分解加速技巧 聂宁明，刘芳，周纯葆 中国科学院计算机网络信息中心, 北京 100190 摘  要： 随着待求解问题的维数的增加以及网格的加密，在数值计算过程中产生的稀疏线性方程组的求解 便成 了计算过程 中最消耗计算资源的模块。在本文中我们根据 问题本身所产生的矩阵的特殊结 构，设计了一种精巧的 PR (Permutation and Reduction) 置换约减算法，可以有效地减少矩阵的维 数，减少矩阵分解过程中的填充，加快了分解的速度。 关键词： 稀疏矩阵；LU 分解；填充；置换；约减 doi: 10.11871/j.issn.1674-9480.2014.03.006 A New Technique for Accelerating the LU Factorization of a Special Class of Sparse Matrix Nie Ningming, Liu Fang, Zhou Chunbao Computer Network Information Center, Chinese Academy Sciences, Beijing 100190, China Abstract: With the increase of the dimension of the computational problems, solving the sparse linear equations is the most time-consuming. In this paper, a new method is presented for a class of matrices with special sparse structure. Based on this, we propose a new method called PR(Permutation and Reduction) method. After performing the permutation and reduction, the original matrix is reduced to a smaller matrix which has the same sparse structure as the corresponding submatrix of the original matrix. The numerical experiments also show the efficiency of the new method. Keywords: Sparse Matrix; LU factorization; fill-in; permutation; reduction 引言 但是在数值计算过程中的核心问题—线性方程组的求 解却一直没有高效的算法来求解。现在学术界和工程 虽然计算机的计算能力在随着摩尔定律增长， 界应用最多的是两大类方法，一类是直接分解法， 基金项目：国家 自然科学基金 91130019) 45 科研信息化技术与应用，2014, 5 (3) [1 2] 如 LU 分解、QR 分解等 - ，另一类是迭代求解法， 1，其余元素为 0 。给定一个 m 个元素的置换，π:{1,