基于服务器并行求解三角形方程组的设计与实现.pdfVIP

计算机工程与设计 ComputerEngineeringandDesign 2010，31(18) 4005 ·软件与算法 · 基于服务器并行求解三角形方程组的设计与实现 刘青昆， 邢 芳 (辽宁师范大学 计算机与信息技术学院，辽宁 大连 l16029) 摘 要：针对在 实际应用 中超出内存容量的大规模矩 阵不可调入的问题 ，设计并实现 了一种面向A矩 阵连续划分的、子 文件大小可调节的求解三角形方程组 的并行方案 。理论上只要在不超过硬盘容量的前提 下，该算法可以在小规模 内存 条件下求解任意规模的矩阵。最后给 出在Altix3700服务器上基于ABEEMtrn模型的数值实验 ，实验结果表明了该算法 的有效性。 关键词：ABEEMcrn模型；分布式存储；三角形方程组；并行算法；MPI；矩阵划分 中图法分类号：TP316 文献标识码 ：A 文章编号：1000．7024(2olo)18—4005．04 Designingandbuildingalgorithm forparallelsolutionoftriangularsystemsbasedonserver LIUQing—kun．XINGFang (SchoolofComputerandInformationTechnology,LiaoningNormalUniversity,Dalian116029，China) Abstract：Aparallelalgorithmisproposed，inwhichit’SAmatrixisdividedintoseveraladjustablesmallermatrixfilestosolvethe problem thatthematrix’Sscaleislargerthanthememory’Scapacity．Theoretically, thisalgorithm maysolvetherandom scalematrix underthesmallconditionofmemory aslongasitdoesnotsurpasstheharddisk’Scapacity．Finally,numericalexperimentsaboutthe algorithm basedonABEEMcrnmodelonAltix3700arecarriedon，anditobtainsthecorrectresultsandconfirmsthealgoritmh ’Svalidity． Keywords：ABEEMcrnmodel；distributed—memory；triangularsolver；parallelalgorithm ；MPI；divisionofMatrix 置量，它表示在每个子文件中，每个处理器计算的行数，np代 O 引 言 表处理器个数，q是划分后的子文件个数，S是控制子文件的循 三角方程组求解是许多现代科研领域计算 问题 的核心， 环变量，在 np：3，q=3，kb=2时行卷帘存储后的情况如图1 因而它是一个吸引众多研究者的课题，近几年来国内外学者 所示 。数据的划分及分配情况如 图2所示 。 针对三角形方程组的并行算法 已经有 了一些显著的成果 [1-5] 2 算法思想 鉴于以上多种算法的出现，使得三角方程的求解速度大幅度 提升，但同时现代科学又不断提出了更高的要求，许多超大规 对超大规模 A阵连续拆分成 q个小文件后存储于某一 模的系数矩阵应运而生，而大部分情况下我们设计的并行计 台节点机上，每次调入 内存的只是其子文件，所 以只要子文 算算法都不去考虑 内存情况，但实际内存却有限，以至于我们 件 的大小在空闲内存的范围内，算法就可以成功运行，以图 虽精心设计了算法却无法实现它，文献[6】提出了一种面向实 1的情形来说明算法 的思想，完整 的算法将在下一节给 出，具 际应用