Fortran+ OpenMP实现实例

PROGRAM parallel_01

    USE omp_lib
    IMPLICIT NONE

    INTEGER :: i,j
    INTEGER() :: time_begin, time_end, time_rate
    REAL, DIMENSION(:,:) :: f, g
    REAL :: k

    WRITE(*,*) '开始进行串行计算'
!>@ 、通过串行计算获得两个矩阵的初始化计算
    CALL system_clock(time_begin,time_rate)
    DO i = ,
        DO j = ,
            f(i,j) = i*j
            k = k +
        END DO
    END DO

    DO i = ,
        DO j = ,
            g(i,j) = i*j +
            k = k +
        END DO
    END DO

    CALL system_clock(time_end,time_rate)
    WRITE(*,*) 'The value of k after serial computing is: ', k
    WRITE(*,*) 'The time wasted on serial computing is: ',(time_end - time_begin)/time_rate

    WRITE(*,*)
    WRITE(*,*)
    WRITE(*,*) '开始进行第一类串行计算—SECTIONS'

!>@ 、通过块并行计算获得两个矩阵的初始化计算
    k =  ! 重新初始化k的值
    CALL system_clock(time_begin,time_rate)
    CALL omp_set_num_threads()
    !$omp parallel
    !$omp sections private(i,j,k)
        !$omp section
             DO i = ,
                DO j = ,
                    f(i,j) = i*j
                    k = k +
                END DO
            END DO
            WRITE(*,*) 'The value of k after parallel computing is: ', k,', and it comes from the thread of ',omp_get_thread_num()

        !$omp section
            DO i = ,
                DO j = ,
                    g(i,j) = i*j +
                    k = k +
                END DO
            END DO
            WRITE(*,*) 'The value of k after parallel computing is: ', k,', and it comes from the thread of ',omp_get_thread_num()
    !$omp end sections
    !$omp end parallel

    CALL system_clock(time_end,time_rate)
    WRITE(*,*) 'The time wasted on the first class parallel computing is: ',(time_end - time_begin)/time_rate
    WRITE(*,*)
    WRITE(*,*)
    WRITE(*,*) '开始进行第二类并行计算—DO'

!>@ 、通过DO循环实现两个矩阵的初始化计算
    k =  ! 重新初始化k的值
    CALL system_clock(time_begin,time_rate)
    !$omp parallel private(k,j,i)
    !$omp do
        DO i = ,
            DO j = ,
                f(i,j) = i*j
                k = k +
                ! 去掉注释后，可现实每一次循环所在的线程ID
 !               WRITE(*,*) 'The value of k after parallel computing is: ', k,', and it comes from the thread of ',omp_get_thread_num()
            END DO
        END DO
    !$omp end do
    !$omp end parallel

    !$omp parallel private(k,j,i)
    !$omp do
        DO i = ,
            DO j = ,
                g(i,j) = i*j
                k = k +
!                WRITE(*,*) 'The value of k after parallel computing is: ', k,', and it comes from the thread of ',omp_get_thread_num()
            END DO
        END DO
    !$omp end do
    !$omp end parallel
    CALL system_clock(time_end,time_rate)
    WRITE(*,*) 'The time wasted on the first class parallel computing is: ',(time_end - time_begin)/time_rate
END PROGRAM

---------------------

本文来自 B325帅猫-量子前沿技术研究所 的CSDN 博客 ，全文地址请点击：https://blog.csdn.net/dengm155/article/details/78837408?utm_source=copy