经典迭代的并行算法.PPTVIP

经典迭代的并行算法 一般形式的Jacobi迭代 正方网格上的Jacobi迭代 超松弛迭代 红黑排序 多色排序 一般形式的Jacobi迭代 考虑线性方程组 正方网格上的Jacobi迭代 考虑二维区域[0,1]×[0,1]上的Laplace方程 正方网格上的Jacobi迭代(续) 正方网格上的Jacobi迭代(续) 对每个偏微分项采用中心差分离散，则有 正方网格上的Jacobi迭代(续) 以n=4为例，A与b分别为 正方网格上的Jacobi迭代(续) 用点Jacobi迭代时， 超松弛迭代 对i=1,2,…,m，块超松弛迭代的公式为 红黑排序 红黑排序(续) 在这种排序下，对应的线性方程组可以写成如下形式 红黑排序(续) 对应的SOR迭代为 红黑排序(续) 正方网格上SOR迭代的串行算法 红黑排序(续) 共享存储并行算法设计方式之一 红黑排序(续) 共享存储并行算法设计方式之二 红黑排序(续) 分布存储并行算法 多色排序 红黑排序下SOR能并行计算的根本原因：网格点进行了分类，任何一类网格点上值的更新只依赖于其它类网格点的值，而同一类网格点间的最新迭代值之间互不依赖 对许多稀疏矩阵，其非零元结构对应的图要比前面所述的情形复杂得多，单纯将网格点分为两类无法达到上述目的 多色排序(续) 多色排序可按下列方式进行： 第一步，置起始颜色值l=1，将结点1的颜色值置为l ； 第二步，记图中与已着色结点相邻的结点集合为S，从S中任意选取一个结点，如果该结点的邻结点颜色值取遍了所有1到l，则置l=l+1，且置该结点的颜色值为l，否则置该结点的颜色值为邻结点尚未取到的最大颜色值； 第三步，如果所有结点都已着色，则停机，否则转第二步继续进行。 多色排序(续) 多色排序(续) 在多色排序下，得到的线性方程组的形式为 多色排序(续) 多色排序下SOR迭代的形式为 * 国防科学技术大学并行与分布处理重点实验室 * 国防科学技术大学并行与分布处理重点实验室 其块Jacobi迭代为 在满足f(x,0)=2x；f(x,1)=2x+1；f(0,y)=y；f(1,y)=2+y 时的求解 在每个方向上取n+2个网格点：0, h, 2h, …, nh, 1，其中h=1/(n+1) 记fi, j=f(ih, jh)，就可以将上述方程组简写为 代入边界条件f0, j＝jh, fn+1, j＝2+jh, fi, 0＝2ih, fi, n+1＝1+2ih,就得到了一个稀疏线性方程组Ax=b 当对某个小于i的j，Ai,j 不为0时，xi的最新迭代值的计算只有在xj的最新迭代值计算出来之后，才能进行 并行化的方法：利用排序技术，使排序后矩阵可以分块为对角块是对角矩阵的形式 SOR2D(f, n, omega) fac1 := 1-omega; fac2 := omega/4; for(k := 1; k = n; k++) { for(j := 1; j = n; j++) fk,j := fac1 * fk,j + fac2 * (fk,j-1 + fk-1,j + fk+1,j + fk,j+1); } SOR2D_omp1(f, n, omega) c1 := 1-omega; c2 := omega/4; #omp parallel for shared(f,n,c1,c2) schedule(static) for(k := 1; k = n; k++) { for(j := 1; j = n; j++) if( (k+j) % 2 == 0) fk,j := c1 * fk,j + c2 * (fk,j-1+ fk-1,j + fk+1,j + fk,j+1); } #omp parallel for shared(f,n,c1,c2) schedule(static) for(k := 1; k = n; k++) { for(j := 1; j = n; j++) if( (k+j) % 2 == 1) fk,j := c1 * fk,j + c2 * (fk,j-1+ fk-1,j + fk+1,j + fk,j+1); } SOR2D_omp2(f, n, omega) c1 := 1-omega; c2 := omega/4; #omp parallel for shared(f,n,c1,c2) schedule(dynamic) for(kj := 2; kj = 2*n; kj = kj+2) { for(k := max(1, kj-n); k = min(n, kj-1); k++) { j := kj – k; fk,j := c1 * fk,j + c2