线性代数方程组树机并行贪心算法研究.pdfVIP

线性代数方程组树机并行贪心算法 杨本立 徐永红 李安志 李方军 中国工程物理研究院工学院，四川绵阳 621900 xuyonghong@ lifangjun@163.net 摘 要 本文利用行处理法贪心方法和分治策略，给出一个基于分布式存储 MIMD 二叉树树机模型求解线 性代数方程组的并行迭代算法，证明该算法的收敛性，并讨论算法的通信复杂度。 关键词 线性代数方程组，行处理法贪心方法，二叉树树机模型，分布式并行迭代算法 1． 前言及预备知识 [1] MIMD 分布式算法和分布共享存储算法是 20 世纪 90 年代以来的一个活跃课题 ．本文利用行处理法贪 心方法[2-5] [6] 和分治策略 给出一个基于分布式存储 MIMD 二叉树树机模型求解线性代数方程组的并行迭代算 法。所给算法具有 “对任意相容性线性代数方程组收敛，附加运算量相对较小，容易嵌入加速技术和可移 植性好”等优点（本文不讨论其加速技术和可移植性）． 称并行算法的串行模拟算法为并行解法．按照 “并行解法+并行计算机=并行算法”的模式安排本文内 容． 本文在实数域中讨论．用（·,·）表示内积， T • 表示欧式范数，（·） 表示转置，[·|·]表示矩 2 阵的列分块划分， • 表示只入不舍取整运算，//…//表示注解．用 自然语句和类PASCAL 语句描述算法．语   句“for each s:1≤s≤q parado I”表示对每个子问题 S 执行同一指令序列 I. 将任意的 n×m 阶线性代数方程组 ∂ X b ( , )  i i  ∂ ≠0 (1-1) i 2  i 1, 2, , n  记为 AX b . （1-1） 称 i 为超平面（∂ ,X）=b 在方程组（1-1）中的序号. i i 定义 1.1 设有 n×m 阶线性代数方程组（1-1）及给定正整数 q(1≤q≤n)，称算法∏G 为 （1-1）的适 应于行处理贪心方法并行解法的一个分组方法。 ∏G : (1) s=1,2,…,q (2) 划分（1-1）成为 q 个子方程组 (∂ ,X ) b S :  j j (1-2) j 1, 2, , j max  使得 i=s+(j-1)·q . 基金项目：中国工程物理研究院科学技术基金资助项目 记子方程组（1-2）为 A X b . (1-2) S S 记分组过程为 (S :A X b ) Π (AX b, q) . (1-3) S S G q [ ] [ ] n  n  易知：∪AS | bS