矩阵方程直接求解算法测试与研究-电磁场与微波技术专业毕业论文.docxVIP

万方数据 万方数据 西安电子科技大学 学位论文独创性（或创新性）声明 秉承学校严谨的学风和优良的科学道德，本人声明所呈交的论文是我个人在 导师指导下进行的研究工作及取得的研究成果。尽我所知，除了文中特别加以标 注和致谢中所罗列的内容以外，论文中不包含其他人已经发表或撰写过的研究成 果；也不包含为获得西安电子科技大学或其它教育机构的学位或证书而使用过的 材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中做了明确的说 明并表示了谢意。 申请学位论文与资料若有不实之处，本人承担一切的法律责任。 本人签名： 日期: 西安电子科技大学 关于论文使用授权的说明 本人完全了解西安电子科技大学有关保留和使用学位论文的规定，即：研究 生在校攻读学位期间论文工作的知识产权单位属西安电子科技大学。学校有权保 留送交论文的复印件，允许查阅和借阅论文；学校可以公布论文的全部或部分内 容，可以允许采用影印、缩印或其它复制手段保存论文。同时本人保证，毕业后 结合学位论文研究课题再撰写的文章一律署名单位为西安电子科技大学。 （保密的论文在解密后遵守此规定） 本学位论文属于保密，在 年解密后适用本授权书。 本人签名： 日期: 导师签名： 日期: 摘要 在计算电磁学中，矩量法产生的大规模稠密矩阵方程的快速求解是矩量法计 算电大问题的关键，而 LU 分解是直接求解线性方程组的一种有效方法。本文主要 测试了一种核外并行 LU 分解程序，该程序以 Linux 系统为平台，大量使用 MPI 库、lapack 库和 scalapack 库中现有的函数，因而计算效率很高。 影响核外并行 LU 分解的因素有四个，分别是矩阵规模的大小、矩阵分块的大 小、进程网格的分配和负载平衡。通过测试发现，该程序中分布式矩阵的分块大 小取 90 左右为宜；进程网格尽量取为接近正方形；随着矩阵规模 N 增大，计算时 间以 O(N3)的速度快速增加。此外，本文还重点对比了并行核外 LU 分解、并行核 内 LU 分解和并行核内 CGNR 三种方法的计算效率。发现求解大规模稠密矩阵方 程时核外并行 LU 算法与迭代算法相比具有很多优势。首先迭代算法的速度要超越 LU 分解方法的前提是矩阵方程的条件数很低，但是大多数情况下矩阵方程的条件 数都很高；其次，核外 LU 分解的速度只比核内 LU 分解稍慢，因此大规模矩阵方 程求解时用硬盘空间代替内存空间不会导致计算时间的急剧增加。 关键词：矩阵方程 LU 分解 并行核外技术 迭代算法 Abstract In computational electromagnetics, the moment method is used to produce large-scale fast solving of dense matrix equation is the key to the problem of moment method for calculation, the LU decomposition is an effective way to solve the linear system of equations directly. In this thesis, we tested a out-of-core parallel LU decomposition procedure, the program based on Linux system platform, extensive use of MPI, lapack libraries and scalapack libraries in the library function, and high computational efficiency. The influencing factors of out-of-core parallel LU decomposition has four, respectively is the size of the matrix, the matrix blocking, processes, and load balancing and the distribution of the grid. By testing found that the distributed block size of the matrix in the program it is advisable to take 90 or so; Try to get as close to the square grid process; As the size of the matrix increasing th