Parallelization of Upwind Scheme

do j=2, maxy-1

do i=2, maxx-1

isgn=1

if (un(i,j).lt.0.d0) then

isgn=-1

end if

jsgn=1

if (vn(i,j),lt.0.d0) then

jsgn=-1

end if

xsgn=isgn

ysgn=jsgn

im1=i-isgn

jm1=j-jsgn

tmp=gy(im1,j) - gy(i,j)

a8=fs(i,j)-fs(im1,j) - fs(i,jm1) + fs(im1,jm1)

a1=gx(im1,j)+gx(i,j) - 2.d0*xsgn*(fs(i,j)-fs(im1,j))

a2=(-a8-xsgn*(gx(i,jm1)-gx(i,j)))*ysgn

・・・・・・・・・・

end do

end do

Compiler can not generate

the optimal data transfer

for fs and gx.

The direction of SHIFT

communication is conditional.

前のｽﾗｲﾄﾞ

次のｽﾗｲﾄﾞ

最初のｽﾗｲﾄﾞに戻る

ｸﾞﾗﾌｨｯｸｽの表示