openmp在cfd上的应用说明.docxVIP

openMP 在 CFD 上的应用简介 1. 编译方法 windows编译：采用pgi visual fortran 2008, 打开OpenM选项编译并行版，关闭 OpenM选项编译串行版 linux 编译： ifort -openmp -o sample *.f90 加-openmp编译并行版，不加-openmp编译串行版(另须删除 omp_set_num_threads等OpenM定义的函数) 2. 程序举例 根据CFD§序的特点，主要是循环处理需要并行 ，以下通过举例说明do循环的并行实现方法，实际上仅这一点基本可以满足 需求。OpenM还有其它一些功能，可详见 OpenM教程。可执行代码见sample.f90 。 program main implicit none integer i, j, a(100), b(100), c(100), d(100,50), n1 integer nThread, tid integer OMP_GET_THREAD_NUM !声明函数 OMP_GET_THREAD_NU该函数为 OpenM定义函数,用于 ！获取当前线程编号 write(unit=6, fmt=*) please input threads number: read(unit=5, fmt=*) nThread if (nThread 1) nThread = 1 call omp_set_num_threads(nThread) !该函数为OpenM定义函数，用于设置并行线程数(核) !!!!!!!!!!!!!!!! 举例 1!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! write(*, *) sample 1: i = 0 !$omp parallel private(tid) ! 并行区域开始 ,private 用于指定局部变量(即每一线程会对 ！该变量做一拷贝，互不影响) tid = OMP_GET_THREAD_NUM() write(unit=6, fmt=100) tid !$omp master write(unit=6, fmt=*) master thread operate ! 只有主线程执行 !$omp end master !$omp critical i = i + 1 !$omp end critical !$omp end parallel ! 并行区域结束 write(*, *) i 100 format(current thread is:, i6) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! * 例注释*：并行区内代码段会被所有线程并行执行一遍。 ! !$omp master 语句 , 只允许主线程执行某段代码。 ! !$omp critical 语句,可控制数据有关联时的执行冲突，但在 CFDi序中不建议使用，最好主动控制数据修改时不关联 !!!!!!!!!!!!!!!! 举例 2!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! write(*, *) sample 2: !$omp parallel private(i, n1)! 指定 i,n1 为局部变量，其它变量（如 a,b,c ）默认为全 ！局变量，被所有线程共享存储 !$omp parallel private(i, n1) !$omp do 将100个循环按一定规则分配给 nThread 个线程并行执行do i = 1, 100 a(i) = i n1 将100个循环按一定规则分配给 nThread 个线程并行执行 !$omp do do i = 1, 99, 2 ! 将50个循环按一定规则分配给 nThread 个线程并行执行 c(i) = a(i) + b(i+1) enddo !$omp end do !$omp end parallel write(*, *) a(1), b(2), c(1), a(99), b(100), c(99) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! * 例注释 * ：并行区内可包括多个循环，依次执行 循环分配规则可选动态、静态和随即等，一般不必指定，取默认即可 !!!!!!!!!!!!!!!! 举例 3!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!