并行代数多重网格算法可扩展性能分析-计算物理.pdf

【文章编号】 1001．246x1200r7)04．0387．08 并行代数多重网格算法可扩展性能分析 徐小文1’2， 莫则尧2 (1．中国工程物理研究院研究生部，北京100明8；2．北京应用物理与计算数学研究所高性能计算中心。北京1000镐) [摘 要】 对当今求解大型稀疏线性代数方程组最有效的迭代方法之一——代数多重网格(AMG)算法的并行 计算进行可扩展性能分析．给出一套并行计算可扩展性能分析方法，用于分析和指导并行迭代算法及实现技术 的设计与优化并应用于并行AMc算法．分析表明，网格算子的平均模式大小和迭代过程的算法效率分别制约了 AMG算法启动阶段和迭代求解阶段并行性能的发挥，成为该类算法急需解决的两个关键问题． 【关键词】 代数多重网格(AMG)；迭代方法；并行计算；可扩展性分析 [中图分类号】 0241．6；’TP338．6 【文献标识码】 A O 引言 随着科学计算需要解决的问题渐趋复杂和并行计算机的迅速发展，并行计算已成为解决大规模科学计 算问题必不可少的手段．在计算数学和计算物理等科学与工程计算领域，很多问题最终都归结为求解稀疏线 性代数方程组．因此，在并行机上高效求解大规模稀疏线性代数方程组是当前科学计算领域一项重要任务． 由于受舍入误差、计算机内存和计算复杂度的限制，对大规模问题，直接求解该类方程组几乎是不可能 的，通常采用迭代方法．迭代法的主要思想是通过构造有效的迭代格式，在有限步数内收敛于方程的精确解， 算法一般结构如下： 给定初始值菇o∈彤； 对于I|}=O，l，2，…，直到{菇‘}收敛： 茗“1=∥(髫‘)； 掘：=詹+l； 其中弗。为初始近似解，G‘为迭代格式，玩为迭代次数．并记单次迭代计算量为c．． 在数值计算领域，人们一直在寻求具有最优计算复杂度的算法，即算法的计算量线性地依赖于物理量个 数【l】．代数多重网格(AMG)方法[2’31就具有这种性质，现已成为数值计算领域一种重要的加速迭代收敛技术． 近年来，在大规模并行机上实现AMG算法成为研究热点H叫训．然而，大量数值实验和实际应用表明[6’9’1驯，与 单机(串行)情形良好的计算复杂度相比，AMG并行计算性能不理想．由于其算法本身的特点，无法事先对其 计算量进行估计，给可扩展性分析带来困难，目前尚无相关文献．较早涉及到AMG并行计算的文献【6‘12盘1虽 然都很强调可扩展性，但只停留在描述层面和基于实验结果的定性分析上．最近的文献(文[9，10，12—14]) 基本上是通过提出不同的并行粗化算法，希望能改善并行AMG算法的算子复杂度和收敛因子． 对一般迭代法而言，算法的最优性可通过迭代次数和单次迭代的计算复杂度来衡量．对于并行迭代算 法，影响其性能的因素还包括单次迭代的并行效率、浮点峰值性能发挥比例(Cache、指令级流水线并行度的 影响)等．因此，如何揭示算法在上述各方面的能力、合理评价其并行性能，对于指导并行迭代算法的设计和 优化有重要作用．现有的性能评价方法大多侧重于并行实现某些方面的性能分析¨“19】，缺乏全面揭示算法 实际性能的能力． 本文提出一套可扩展性能分析方法，目的是用于分析和指导并行迭代算法及实现技术的设计与优化．首 [收稿日期】2∞6一03一13；【修回日期】2006—09一∞ [基金项目】国家杰出青年基金(604252晒)、国家9r73项目(2005cB321702)和NsF(No资助项目 【作者简介】徐小文(1978一)，男，湖南。博士生，从事代数多重网格方法和并行计算方面的研究． 388 计 算 物 理 第24卷 先，在文[18]基础上，给出一套性能评价方法．该方法将算法计算复杂度的评价从原有的准则中分离出来，将 传统的并行计算效率进行分解，这对于在当今科学计算领域中强调算法本身的可扩展性具有重要意义．然 后，基于具有最优计算复杂度的串行算法，提出最优并行计算时间的概念，并采用该时间被放大的倍数，即时 间扩展因子，来刻画并行迭代算法及其实现的可扩展能力．类似于效率分解，时间扩展因子也可进行相应地 分解．基于此，可全面