并行处理总复习–130421.docVIP

《并行处理》总复习 一、 1.1、用向量语言VFORTRAN，把一维数组A（0：1000）的所有元素均置成PAI； A（0：1000）= PAI 1.2、用向量语言VFORTRAN，把一维数组A（-1：1000）的第偶数号元素减少1； 1.3、用向量语言VFORTRAN，把一维数组A（0：1000）的第奇数号元素置成其原值立方的正弦； A（1：1000:2） = sin( A（1：1000:2）**3) 1.4、已知一维数组A（－1000：999）。对于每个元素，倘若其正切值不超过立方值，则令原值减半；否则，令原值加倍。用向量语言VFORTRAN的语句，描述之。 1.5、在下面的 PFORTRAN 主程序段里，标记出并行语句和并行进程名，并指出可能的最大并行数目。 1 program matmul 2 external init_matrix, mul_matrix, prt_matrix 3 shared_distributed c : block 4 parameter(SIZE=32768) 5 integer a(SIZE,SIZE), b(SIZE,SIZE), c(SIZE,SIZE) 6 call init_matrix(a,b) 7 call m_set_procs(SIZE÷32) 8 call mg_fork(mul_matrix(a,b)) 9 call m_kill_procs() 10 call prt_matrix(a,b,c) 11 end 1.6、在分辨率为SIZE*SIZE(SIZE为2的幂)的黑白带灰度的显示屏幕上：左上角为坐标原点，Ｘ轴向下，Ｙ轴向右；０为白，１为黑。用NUMNODES(亦为2的幂，且SIZE/4)个并行节点把现有的图象(位图表示)map[0:SIZE-1,0:SIZE-1]加以变换。 Map的初值：map[i][j]=(i+j)/(2*(SIZE-1))……0≦i,jSIZE 主进程： for (I=0,row=0;INUMNODES;I++,row+=SIZE/NUMNODES) send(row,Pi); for (I=0;ISIZE**2;I++) { recv(oldrow,oldcol,newvalue,Pany); temp_map[oldrow][oldcol] = newvalue;} for (I=0;ISIZE;I++) for (j=0;jSIZE;j++) map[I][j]=temp_map[I][j]; 从进程： recv(row,Pmaster); for (oldrow=row;oldrow(row+SIZE/NUMNODES);oldrow++) for (oldcol=0;oldcolSIZE;oldcol++) {newrow=oldrow÷4*4; newcol=oldcol÷4*4; /* ÷ 是整除运算符*/ newvalue=map[newrow][newcol]; send (oldrow,oldcol,newvalue,Pmaster);} 6.1、在卷面上标记出消息发送语句； 6.2、在卷面上标记出消息接收语句，以及与之配对的消息发送语句。 6.3、指出新图上灰度为１的点； 6.4、指出新图上灰度为０的点； 6.5、给出新图自左上至右下对角线的图象； 6.6、描述所进行的图形／图象变换； 6.7、跟例2.14.1相比，哪个变换更接近原图？ 二、限定语句视图。画出优先关系图，判别可并行性。若可，给出目标语序，以及 VFORTRAN 目标程序。 2.1、 FOR ( I = M ; I .LT. N ; I++ ) A(I) = A(I+1)**3 2.2、 FOR ( I = M ; I .LT. N ; I++ ) A(I) = A(I-1)**6 2.3、 FOR ( I = M ; I .GE. N ; I-- ) C(I) = COS (C(I-2)**2 ) 2.4、 FOR ( J = M-5 ; J .LT. N+5 ; J++ ) { #A(J+1)=A(J+1) #A a B(2*J) = #A(j+1) B #a A(J-1) = B(2*J+2) A b } 解:结点命名如上.同名数组出现对A,a:下标式之差(J-1)-(J+1)=-2,异号,a-A B,b:下标式之差(2*J)-(2*J+2)=-2,异号,b-B 优先图 无矛盾回路,可并行;目标语序1,2 和2,1; 挑2,1,存在