并行程序设计基础.pptVIP

计算圆周率的c语言代码段#defineN1000000main(){ doublelocal,pi=0.0,w; longi; w=1.0/N; for(i=0;iN;i++){ local=(i+0.5)*w; pi=pi+4.0/(1.0+local*local); } printf(“piis%f\n”,pi*w);}**第30页，共45页，星期日，2025年，2月5日并行程序设计基础12.1并行程序设计概述12.2进程12.3线程12.4同步12.5通信12.6并行程序设计模型**第31页，共45页，星期日，2025年，2月5日进程进程的基本概念进程的并行执行进程的相互作用**第32页，共45页，星期日，2025年，2月5日并行程序设计基础第1页，共45页，星期日，2025年，2月5日并行程序设计基础12.1并行程序设计概述12.2进程12.3线程12.4同步12.5通信12.6并行程序设计模型**第2页，共45页，星期日，2025年，2月5日并行程序设计概述并行程序设计难的原因并行语言的构造方法并行性问题交互/通信问题五种并行编程风范计算圆周率的样本程序**第3页，共45页，星期日，2025年，2月5日1并行程序设计难的原因技术先行,缺乏理论指导程序的语法/语义复杂,需要用户自已处理任务/数据的划分/分配数据交换同步和互斥性能平衡并行语言缺乏代可扩展和异构可扩展,程序移植困难,重写代码难度太大环境和工具缺乏较长的生长期,缺乏代可扩展和异构可扩展**第4页，共45页，星期日，2025年，2月5日2并行语言的构造方法串行代码段for(i=0;iN;i++)A[i]=b[i]*b[i+1];for(i=0;iN;i++)c[i]=A[i]+A[i+1];(a)使用库例程构造并行程序id=my_process_id();p=number_of_processes();for(i=id;iN;i=i+p)A[i]=b[i]*b[i+1];barrier();for(i=id;iN;i=i+p)c[i]=A[i]+A[i+1];例子:MPI,PVM,Pthreads(b)扩展串行语言my_process_id,number_of_processes(),andbarrier()A(0:N-1)=b(0:N-1)*b(1:N)c=A(0:N-1)+A(1:N)例子:Fortran90(c)加编译注释构造并行程序的方法#pragmaparallel#pragmashared(A,b,c)#pragmalocal(i){#pragmapforiterate(i=0;N;1)for(i=0;iN;i++)A[i]=b[i]*b[i+1];#pragmasynchronize#pragmapforiterate(i=0;N;1)for(i=0;iN;i++)c[i]=A[i]+A[i+1];}例子：SGIpowerC**第5页，共45页，星期日，2025年，2月5日三种并行语言构造方法比较2并行语言的构造方法**第6页，共45页，星期日，2025年，2月5日3并行性问题3.1进程的同构性SIMD:所有进程在同一时间执行相同的指令MIMD:各个进程在同一时间可以执行不同的指令SPMD:各个进程是同构的，多个进程对不同的数据执行相同的代码(一般是数据并行的同义语)常对应并行循环，数据并行结构，单代码MPMD:各个进程是异构的，多个进程执行不同的代码（一般是任务并行，或功能并行，或控制并行的同义语）常对应并行块，多代码要为有1000个处理器的计算机编写一个完全异构的并行程序是很困难的**第7页，共45页，星期日，2025年，2月5日并行块parbeginS1S2S3…….SnparendS1S2S3…….Sn可以是不同的代码并行循环:当并行块中所有进程共享相同代码时parbeginS1S2S3…….SnparendS1S2S3…….Sn是相同代码简化为parfor(i=1;i=n,i++)S(i)进程的同构性3并行性问题**第8页，共45页，星期日，2025年，2