accelerate.f90

!Crown Copyright 2012 AWE.
!
! This file is part of CloverLeaf.
!
! CloverLeaf is free software: you can redistribute it and/or modify it under 
! the terms of the GNU General Public License as published by the 
! Free Software Foundation, either version 3 of the License, or (at your option) 
! any later version.
!
! CloverLeaf is distributed in the hope that it will be useful, but 
! WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 
! FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more 
! details.
!
! You should have received a copy of the GNU General Public License along with 
! CloverLeaf. If not, see http://www.gnu.org/licenses/.

!>  @brief Driver for the acceleration kernels
!>  @author Wayne Gaudin
!>  @details Calls user requested kernel

MODULE accelerate_module

CONTAINS

  SUBROUTINE accelerate()

    USE clover_module
    USE accelerate_kernel_module

    IMPLICIT NONE

    INTEGER :: tile

    REAL(KIND=8) :: kernel_time,timer

    IF(profiler_on) kernel_time=timer()
  
    
    IF(use_fortran_kernels) THEN
      DO tile=1,tiles_per_chunk


        CALL accelerate_kernel(chunk%tiles(tile)%t_xmin,                &
          chunk%tiles(tile)%t_xmax,                  &
          chunk%tiles(tile)%t_ymin,                  &
          chunk%tiles(tile)%t_ymax,                  &
          dt,                                     &
          chunk%tiles(tile)%field%xarea,                  &
          chunk%tiles(tile)%field%yarea,                  &
          chunk%tiles(tile)%field%volume,                 &
          chunk%tiles(tile)%field%density0,               &
          chunk%tiles(tile)%field%pressure,               &
          chunk%tiles(tile)%field%viscosity,              &
          chunk%tiles(tile)%field%xvel0,                  &
          chunk%tiles(tile)%field%yvel0,                  &
          chunk%tiles(tile)%field%xvel1,                  &
          chunk%tiles(tile)%field%yvel1                  )


      ENDDO
    ELSEIF(use_C_kernels)THEN
      DO tile=1,tiles_per_chunk
        CALL accelerate_kernel_c(chunk%tiles(tile)%t_xmin,                &
          chunk%tiles(tile)%t_xmax,                  &
          chunk%tiles(tile)%t_ymin,                  &
          chunk%tiles(tile)%t_ymax,                  &
          dt,                                     &
          chunk%tiles(tile)%field%xarea,                  &
          chunk%tiles(tile)%field%yarea,                  &
          chunk%tiles(tile)%field%volume,                 &
          chunk%tiles(tile)%field%density0,               &
          chunk%tiles(tile)%field%pressure,               &
          chunk%tiles(tile)%field%viscosity,              &
          chunk%tiles(tile)%field%xvel0,                  &
          chunk%tiles(tile)%field%yvel0,                  &
          chunk%tiles(tile)%field%xvel1,                  &
          chunk%tiles(tile)%field%yvel1                  )
      ENDDO
    ENDIF

    IF(profiler_on) profiler%acceleration=profiler%acceleration+(timer()-kernel_time)

  END SUBROUTINE accelerate

END MODULE accelerate_module