diff options
author | upstream source tree <ports@midipix.org> | 2015-03-15 20:14:05 -0400 |
---|---|---|
committer | upstream source tree <ports@midipix.org> | 2015-03-15 20:14:05 -0400 |
commit | 554fd8c5195424bdbcabf5de30fdc183aba391bd (patch) | |
tree | 976dc5ab7fddf506dadce60ae936f43f58787092 /gcc/testsuite/gfortran.dg/vect/vect-gems.f90 | |
download | cbb-gcc-4.6.4-15d2061ac0796199866debe9ac87130894b0cdd3.tar.bz2 cbb-gcc-4.6.4-15d2061ac0796199866debe9ac87130894b0cdd3.tar.xz |
obtained gcc-4.6.4.tar.bz2 from upstream website;upstream
verified gcc-4.6.4.tar.bz2.sig;
imported gcc-4.6.4 source tree from verified upstream tarball.
downloading a git-generated archive based on the 'upstream' tag
should provide you with a source tree that is binary identical
to the one extracted from the above tarball.
if you have obtained the source via the command 'git clone',
however, do note that line-endings of files in your working
directory might differ from line-endings of the respective
files in the upstream repository.
Diffstat (limited to 'gcc/testsuite/gfortran.dg/vect/vect-gems.f90')
-rw-r--r-- | gcc/testsuite/gfortran.dg/vect/vect-gems.f90 | 58 |
1 files changed, 58 insertions, 0 deletions
diff --git a/gcc/testsuite/gfortran.dg/vect/vect-gems.f90 b/gcc/testsuite/gfortran.dg/vect/vect-gems.f90 new file mode 100644 index 000000000..66e878d3d --- /dev/null +++ b/gcc/testsuite/gfortran.dg/vect/vect-gems.f90 @@ -0,0 +1,58 @@ +! { dg-do compile } +! { dg-require-effective-target vect_double } + +MODULE UPML_mod + +IMPLICIT NONE + +PUBLIC UPMLupdateE + +PRIVATE + +real(kind=8), dimension(:,:,:), allocatable :: Dx_ilow + +real(kind=8), dimension(:), allocatable :: aye, aze +real(kind=8), dimension(:), allocatable :: bye, bze +real(kind=8), dimension(:), allocatable :: fxh, cxh + +real(kind=8) :: epsinv +real(kind=8) :: dxinv, dyinv, dzinv + +integer :: xstart, ystart, zstart, xstop, ystop, zstop + +CONTAINS + +SUBROUTINE UPMLupdateE(nx,ny,nz,Hx,Hy,Hz,Ex,Ey,Ez) + +integer, intent(in) :: nx, ny, nz +real(kind=8), intent(inout), & + dimension(xstart:xstop+1,ystart:ystop+1,zstart:zstop+1) :: Ex, Ey, Ez +real(kind=8), intent(inout), & + allocatable :: Hx(:,:,:), Hy(:,:,:), Hz(:,:,:) + +integer :: i, j, k +real(kind=8) :: Dxold, Dyold, Dzold + +do k=zstart+1,zstop + do j=ystart+1,ystop + do i=xstart+1,0 + + Dxold = Dx_ilow(i,j,k) + + Dx_ilow(i,j,k) = aye(j) * Dx_ilow(i,j,k) + & + bye(j) * ((Hz(i,j,k )-Hz(i,j-1,k))*dyinv + & + (Hy(i,j,k-1)-Hy(i,j,k ))*dzinv) + + Ex(i,j,k) = aze(k) * Ex(i,j,k) + & + bze(k) * (cxh(i)*Dx_ilow(i,j,k) - fxh(i)*Dxold) * epsinv + end do + end do +end do + +END SUBROUTINE UPMLupdateE + +END MODULE UPML_mod + +! { dg-final { scan-tree-dump-times "vectorized 1 loops" 1 "vect" } } +! { dg-final { cleanup-tree-dump "vect" } } +! { dg-final { cleanup-modules "upml_mod" } } |