并行算法实践.pptVIP

国家高性能计算中心（合肥） * * 共享任务结构 共享任务结构将它所包含的代码划分给线程组的各成员来执行 并行for循环 并行sections 串行执行 国家高性能计算中心（合肥） * * for编译制导语句 for语句指定紧随它的循环语句必须由线程组并行执行； 语句格式 #pragma omp for [clause[[,]clause]…] newline [clause]= Schedule(type [,chunk]) ordered private (list) firstprivate (list) lastprivate (list) shared (list) reduction (operator: list) nowait 国家高性能计算中心（合肥） * * for编译制导语句 schedule子句描述如何将循环的迭代划分给线程组中的线程 如果没有指定chunk大小，迭代会尽可能的平均分配给每个线程 type为static，循环被分成大小为 chunk的块，静态分配给线程 type为dynamic,循环被动态划分为大小为chunk的块，动态分配给线程 国家高性能计算中心（合肥） * * Sections编译制导语句 sections编译制导语句指定内部的代码被划分给线程组中的各线程 不同的section由不同的线程执行 Section语句格式： #pragma omp sections [ clause[[,]clause]…] newline { [#pragma omp section newline] … [#pragma omp section newline] … } 国家高性能计算中心（合肥） * * Sections编译制导语句 clause= private (list) firstprivate (list) lastprivate (list) reduction (operator: list) nowait 在sections语句结束处有一个隐含的路障，使用了nowait子句除外 国家高性能计算中心（合肥） * * Sections编译制导语句 #include omp.h #define N 1000 int main (){ int i; float a[N], b[N], c[N]; /* Some initializations */ for (i=0; i N; i++) a[i] = b[i] = i * 1.0; #pragma omp parallel shared(a,b,c) private(i) { #pragma omp sections nowait { #pragma omp section for (i=0; i N/2; i++) c[i] = a[i] + b[i]; #pragma omp section for (i=N/2; i N; i++) c[i] = a[i] + b[i]; } /* end of sections */ } /* end of parallel section */ } 国家高性能计算中心（合肥） * * single编译制导语句 single编译制导语句指定内部代码只有线程组中的一个线程执行。 线程组中没有执行single语句的线程会一直等待代码块的结束，使用nowait子句除外 语句格式： #pragma omp single [clause[[,]clause]…] newline clause= private(list) firstprivate(list) nowait 国家高性能计算中心（合肥） * * 组合的并行共享任务结构 parallel for编译制导语句 parallel sections编译制导语句 国家高性能计算中心（合肥） * * parallel for编译制导语句 Parallel for编译制导语句表明一个并行域包含一个独立的for语句 语句格式 #pragma omp parallel for [clause…] newline clause= if (scalar_logical_expression) default (shared | none) schedule (type [,chunk]) shared (list) private (list) firstprivate (list) last