Subversion Repositories openrisc_me

! { 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" } }