scalar_residual_cpu.f90

Go to the documentation of this file.

 
module scalar_residual_cpu
  use scalar_residual, only : scalar_residual_t
  use ax_product, only : ax_t
  use field, only : field_t
  use coefs, only : coef_t
  use space, only : space_t
  use mesh, only : mesh_t
  use num_types, only : rp
  use field, only : field_t
  use mesh, only : mesh_t
  use ax_product, only : ax_t
  use space, only : space_t
  use coefs, only : coef_t
  implicit none
  private
 
  type, public, extends(scalar_residual_t) :: scalar_residual_cpu_t
   contains
     procedure, nopass :: compute => scalar_residual_cpu_compute
  end type scalar_residual_cpu_t
 
contains
 
  subroutine scalar_residual_cpu_compute(Ax, s, s_res, f_Xh, c_Xh, msh, Xh, &
       lambda, rho_cp, bd, dt, n)
    class(ax_t), intent(in) :: Ax
    type(mesh_t), intent(inout) :: msh
    type(space_t), intent(inout) :: Xh
    type(field_t), intent(inout) :: s
    type(field_t), intent(inout) :: s_res
    type(field_t), intent(in) :: f_Xh
    type(coef_t), intent(inout) :: c_Xh
    type(field_t), intent(in) :: lambda, rho_cp
    real(kind=rp), intent(in) :: bd
    real(kind=rp), intent(in) :: dt
    integer, intent(in) :: n
    integer :: i
 
    !OCL NORECURRENCE, NOVREC, NOALIAS
    !DIR$ CONCURRENT
    !DIR$ IVDEP
    !GCC$ ivdep
    !$omp parallel do
    do i = 1, n
       c_xh%h1(i,1,1,1) = lambda%x(i,1,1,1)
       c_xh%h2(i,1,1,1) = rho_cp%x(i,1,1,1) * bd / dt
    end do
    !$omp end parallel do
 
    c_xh%ifh2 = .true.
 
    call ax%compute(s_res%x, s%x, c_xh, msh, xh)
 
    !OCL NORECURRENCE, NOVREC, NOALIAS
    !DIR$ CONCURRENT
    !DIR$ IVDEP
    !GCC$ ivdep
    !$omp parallel do
    do i = 1, n
       s_res%x(i,1,1,1) = (-s_res%x(i,1,1,1)) + f_xh%x(i,1,1,1)
    end do
    !$omp end parallel do
 
  end subroutine scalar_residual_cpu_compute
 
end module scalar_residual_cpu