一种基于动态序列的单边Jacobi方法.pdfVIP

第 31卷第4期 苏 州 大 学 学 报 (工 科 版) V01．31 No．4 2011年 8月 JOURNALOFSOOCHOW UNIVERSITY(ENGINEERINGSCIENCEEDITION) Aug．2011 文章编号：1673—047X(2011)一04—0016—07 一 种基于动态序列的单边 Jacobi方法 郭 强，赵 雷 (苏州大学计算机科学与技术学院，江苏 苏州215006) 摘 要：采用Jacobi方法并行求解矩阵奇异值有多种数据交换序列，在双边 Jacobi方法中，采用动 态序列要比静态循环序列更加高效，可以将其应用到单边Jacobi方法中。为了在每一次迭代开始时 动态生成数据交换序列，首先计算矩阵子块间的谱范数，然后对这些谱范数形成的完全图应用最大 权完美匹配算法，最终结果作为各计算节点传递数据的依据。实验表明谱范数可以很好地表示矩阵 列对之间的正交程度，将其应用在求解动态序列的过程中，使得单边Jacobi方法计算矩阵奇异值分 解更加高效。 关键词：单边Jacobi算法；奇异值分解；谱范数；动态序列；最大权完美匹配 中图分类号：TP311 文献标识码 ：A 0 引 言 假设计算节点个数为P，将矩阵A按列划分成 =2p块 ，各个节点分别对 自身获得的两个列块A。，和A 进行正交变换，然后按照文献[1]给出的一些交换序列传递数据，迭代 f或Z一1次后，矩阵中所有列与列彼此 之间相互正交一次，完成一轮变换。经过若干轮变换后，当矩阵满足收敛条件时，整个计算过程结束。文献 [2]称上述交换序列为静态序列。为了加快奇异值分解的速度，文献[2]提出一种新的交换序列，称为动态序 列，并将其应用在双边Jacobi方法求解矩阵奇异值分解中，取得了良好的实验效果。结合单边Jacobi法与双 边Jacobi法之间的联系，文献[3—4]首次将动态序列引入到单边分块Jacobi方法求解矩阵奇异值分解中，作 者提出采用列块间所有列对夹角余弦值之和或者最大列对夹角余弦值来表示列块间偏离正交的程度，并将这 些数值作为完全图中的边权重，然后应用最大权完美匹配算法生成数据交换序列。但是由于每一个节点仅包 含两个列块，要想在每次正交变换之前得到所有列块之间的偏离程度，节点间需要交换数据，通信时间增大， 对此文献[3—4]没有给出最终的解决方案 。 本文采用矩阵列块间的谱范数来表示完全图中边的权重，并且设计一种并行幂法来计算上述谱范数，新 算法不需在节点间传递整个列块数据，可以节省大量通信时间。最后通过实验验证，该方法用于求解动态序 列效果 良好 ，可以加速矩阵列块间的正交化过程。 1 背景知识 在求解奇异值分解之前，如果将A 和A 作如下变换写成： 日兰(＼日21日22)／=fA三【AI：A三TAAARR]J 对A 和A 进行正交变换的过程就变为隐式消去 日中非对角块 日。与日：。的过程。这时对A采用单边Jacobi 收稿 日期：2011—05—03 作者简介：郭 强(1985一)，男，硕士研究生，主要研究方向为并行计算。 基金项目：国家 自然科学基金资助项 目(编。 第 4期 郭 强，赵 雷：一种基于动态序列的单边Jacobi方法 方法求解奇异值分解的过程变为对 A采用双边Jacobi方法求解特征值分解的过程，因此将文献 [2]中的动 态序列引入到单边Jacobi方法中。同文献[2]将ATA中所有子列块的罗贝尼乌斯范数I1日iilIF(，=1，…，z，i )作为完全图的初始数据，但是各节点需要相互传递 自身的A 和A 进行矩阵乘法，通信量大。通过观察 知道当对矩阵进行正交变换若干轮后，各子块的弗罗贝尼乌斯范数l1日ijIIF=~／+…+；— ，即l1日iiII一 or，而I1日ijIl：=or，因此用谱范数J1日ijJ：『替换 I1日ijIl作为求解动态序列时完全图的初始数据。 当求解出所有 l1日iiIl：( )之后，就会得到一个拥有z个节点、z(z一1)／2条边的无向完全图G=(V，E)。 如图1所示，其中节点 ∈V表示列块i，无向边( ，)∈E的权重为 l1日ilI。为了加快