分裂与并行计算作业.docxVIP

线性代数方程组串行与并行求解方法 刘争光 2013：L：L27：L 1.问题 考虑一般的线性方程组Au=bf最普遍的方法是先对A做厶〃分解，其中L为 单位下三角阵，U是上三角阵。于是方程组化为LUx = b,率先求解Ly=h这 个三角形方程组，然后再求解t/x = yo 串行算法: (1) A的LU分解 采用部分选主元Gauss消去法进行列消元，使得L是单位下三角阵。算法中 △k表示A的第K行。 算法: for j=O,n-2 选主元，找到 二 maxL.. ,z = j, n 1} 如果厶那么交换儿和￡ Hjj =aij/ajj,i = j + 1,……n-l for k=j+l,n-l %k =aik —QijXa从,i = j + l,……n-\ end for end for (3)计算Ly=bfUx = y解法类似。 算法： for i=0,n-l 儿=bilu for j=i+l,n-l 5 =切- 1卩x为 end for end for 并行算法 (1)LU分解算法： 矩阵A采用卷帘方式存储，即把矩阵A的第i列存放在PinilxJp中o 匕嗣的算法如下： icol - 0 forj=0, n?2 如果 myid =mod(j,p),那么 找到列主元素所在行 选主元，找到/：|a/Jcol| = max仏扁异=j, 1} 如果/ H久那么交换勺?和a】 ai,icoi =⑷问 /勺问，i = j + 1，……〃 _ 1 /尸叽，W + 1，?……1 把I和f广播给其他进程，为以后进行行交换 icol=icol+l 所有进程进行行交换 即:如果Zf,那么交换仏和A for k=icol,m-l %k =aik -fiXaA J = j + 1,……n-\ end for end for (2)流程图 结果(A的LU分解) 分析：由于核数为2时，通信时间占据大部分并行时间，故效率很低。 大小 100*100 160*160 200*200 串行 8.0*10-3 1.6*10-2 3.2*10-2 时间 效率 时间 效率 时间 效率 2 7.9*10-3 0.506 1.2*10-2 0.66 2.8*10-2 0.57 4 2.2*10-3 0.91 4.3*10-3 0.93 8.8*10-3 0.91 5.附录 (1)串行程序: !The program is to solve the equation AX二F !This is a serial program !subrout inc in program named LU is to solve LU decomposition for A ! program main implicit none interface subroutine LU(n, A, F) implicit none integer n real::A(n, n), F(n) real: :x(n), y (n) real maxx, b, c end subroutine end interface !define A is a n*n matrix integer, parameter::n=80 rea 1:: A(n, n), AA(n, n), F(n) integer i, j real:: x (n), y (n) real tl, t2 (22)(23) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) I when i=j, A(i, j)=0, else A(i, j)=l ! then we choose F(i)二nT !so of course the exact solution is x=(l) open(10, file二AAAA. txt) open(ll, file二〃FFFF. txt) read (10, *) AA read(l 1,*)F do i=l, n do j=l, n A(j, i)二AA(i, j) end do encl do call cpu_time(tl) call LU(n, A, F) cal 1 cpu_time(t2) Iget A二LU {firstly, solve Ly二F !of course L is a under unit tr