Optimized code

!hpf$ independent,new(i,isgn,jsgn,xsgn,ysgn,im1,jm1,tmp,

!hpf$& a1,a2,a3,a4,a5,a6,a7,a8,xx,yy)

do j=2, maxy-1

do i=2, maxx-1

isgn=1

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

isgn=-1

endif

jsgn=1

fs_im1jm1 = fs(im1, j-1)

fs_ijm1 = fs(i , j-1)

gy_ijm1 = gy(i , j-1)

gx_ijm1 = gx(i , j-1)

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

jsgn=-1

fs_im1jm1 = fs(im1, j+1)

fs_ijm1 = fs(i , j+1)

gy_ijm1 = gy(i , j+1)

gx_ijm1 = gx(i , j+1)

endif

xsgn=isgn

ysgn=jsgn

im1=i-isgn

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

a8=fs(i,j)-fs(im1,j)-fs_ijm1+fs_im1jm1

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

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

・・・・・・・・・・

Copy the data to be shifted into

scalar temporary variable.

HPF/JA: no code modification needed

!hpf$ refrect fs,gx,gy,vn,un

!hpf$ independent,new(i,isgn,jsgn,xsgn,ysgn,im1,jm1,tmp,

!hpf$& a1,a2,a3,a4,a5,a6,a7,a8,xx,yy)

do j=2, maxy-1

do i=2, maxx-1

!hpf$ on home(un(i,j)), local begin

・・・・・・・…

loop body

・・・・・・・・・・

!hpf$ end on

end do

end do

!hpf$ refrect fn,gx1,gy1

前のｽﾗｲﾄﾞ

次のｽﾗｲﾄﾞ

最初のｽﾗｲﾄﾞに戻る

ｸﾞﾗﾌｨｯｸｽの表示