高性能计算教案-第九讲-last.ppt

ordered 示例 void work(int k) { #pragma omp ordered printf( %d, k); } …… #pragma omp for ordered schedule(dynamic) for (i = lb; i ub; i += st) work(i); threadprivate 指导语句 使一个全局文件作用域的变量在并行域内变成每个线程私有。 #pragma omp threadprivate (list) newline 每个线程对该变量复制一份私有拷贝。 threadprivate 示例 int alpha[10], beta[10], i; #pragma omp threadprivate(alpha) int main () { // 第一个并行区域 #pragma omp parallel private(i,beta) for (i=0; i 10; i++) alpha[i] = beta[i] = i; // 第二个并行区域 #pragma omp parallel printf(alpha[3]= %d and beta[3]=%d\n, alpha[3], beta[3]); } 数据域属性子句 用于指定变量的作用域范围。 private子句； shared子句； default子句； firstprivate子句； lastprivate子句； copyin子句； reduction子句； private子句 表示它列出的变量（可能局部）对于每个线程是局部的。 语句格式为：private(list)。 ? 数据类型 位置 持久性 扩充性 初始化 private 变量 在域的开始或共享任务单元 否 只是词法的，除非作为子程序 的参数而传递 使用firstprivate threadprivate 变量 在块或整个文件区域的例程定义上 是 动态的 使用copyin private和threadprivate区别 shared子句 表示它所列出的变量被线程组中所有的线程共享所有线程都能对它进行读写访问。 语句格式为：shared (list)。 default子句 用户自行规定在一个并行域的静态范围中所定义的变量的缺省作用范围。 default (shared | none)。 int x, y, z[1000]; #pragma omp threadprivate(x) void fun(int a) { const int c = 1; int i = 0; #pragma omp parallel default(none) private(a) shared(z) { int j = omp_get_num_thread(); // 正确！因为j声明在并行区域内 a = z[j]; // 正确！因为a和z分别被列在private和shared子句中 x = c; // 正确！因为x是threadprivate，c是常量限定类型 z[i] = y; // 错误！不能在这里引用i或y } } firstprivate子句 是private子句的超集，用于对变量做原子初始化。 firstprivate (list)。 incr = 0; #pragma omp parallel for firstprivate(incr) for ( i = 0; i MAX; i++) { if ( (i%2) == 0) incr++; a[i] = incr; } lastprivate子句 也是private子句的超集，用于将变量从最后的循环迭代或段复制给原始的变量。 lastprivate (list)。 void sq2(int n, double *last) { double x; #pragma omp parallel #pragma omp for lastprivate(x) for (int i = 0; i n; i++) { x = a[i] * a[i] + b[i] * b[i]; b[i] = sqrt(x); } last = x; } copyin子句 用来为线程组中所有线程的threadprivate变量赋相同的值，主线程该变量的值作为初始值。 语句格式为：copyin(list)。 reduction子句 使用指定的操作对其列表中出现的变量进行归; 初始时，每个线程都保留一份私有拷贝，在结构尾部根据指定的操作对线程中的相应变量进行归约，并更新该变量的全局值。 语句格式为：reduction (operator: list) #pragma omp pa