24小时热门版块排行榜    

查看: 610  |  回复: 1

jiedati123

新虫 (初入文坛)

[求助] 请问下面这段代码怎么进行openmp并行 已有1人参与

DO JJ=2,52
                J=JT2-JJ
                PT(ISTF)=0.
                QT(ISTF)=F(ISTF,J,N)
                DO I=3,152
                    DENOM=AP(I,J)-PT(I-1)*AIM(I,J)
                    PT(I)=AIP(I,J)/DENOM
                    TEMP=CON(I,J)+AJP(I,J)*F(I,J+1,N)+AJM(I,J)*F(I,J-1,N)
                    QT(I)=(TEMP+AIM(I,J)*QT(I-1))/DENOM
                ENDDO
                DO II=3,52
                    I=IT1-II
                    F(I,J,N)=F(I+1,J,N)*PT(I)+QT(I)
                ENDDO
            ENDDO
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cntsinghua

银虫 (初入文坛)

【答案】应助回帖

!$OMP PARALLEL
DO JJ=2,52
                J=JT2-JJ
                PT(ISTF)=0.
                QT(ISTF)=F(ISTF,J,N)

!$OMP DO PRIVATE(DENOM,TEMP)
                DO I=3,152
                    DENOM=AP(I,J)-PT(I-1)*AIM(I,J)
                    PT(I)=AIP(I,J)/DENOM
                    TEMP=CON(I,J)+AJP(I,J)*F(I,J+1,N)+AJM(I,J)*F(I,J-1,N)
                    QT(I)=(TEMP+AIM(I,J)*QT(I-1))/DENOM
                ENDDO
你试一下
!$OMP END DO

!$OMP DO PRIVATE(I)
                DO II=3,52
                    I=IT1-II
                    F(I,J,N)=F(I+1,J,N)*PT(I)+QT(I)
                ENDDO
!$OMP END DO
            ENDDO
!$OMP END PARALLEL
2楼2018-02-09 10:24:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jiedati123 的主题更新
信息提示
请填处理意见