prgm_lib.F90

Go to the documentation of this file.

 
module opencl_prgm_lib
  use opencl_intf
  use utils, only : neko_error
  use, intrinsic :: iso_c_binding, only : c_ptr, c_null_ptr
  implicit none
  private
 
#ifdef HAVE_OPENCL
 
  type(c_ptr), public, bind(c) :: math_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: mathops_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: dirichlet_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: inflow_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: zero_dirichlet_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: symmetry_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: facet_normal_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: neumann_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: inhom_dirichlet_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: dudxyz_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: cdtp_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: conv1_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: convect_scalar_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: set_convect_rst_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: cfl_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: opgrad_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: gs_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: ax_helm_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: ax_helm_full_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: jacobi_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: rhs_maker_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: pnpn_res_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: pnpn_stress_res_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: euler_res_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: compressible_ops_compute_max_wave_speed_program = c_null_ptr
  type(c_ptr), public, bind(c) :: compressible_ops_compute_entropy_program = c_null_ptr
  type(c_ptr), public, bind(c) :: compressible_ops_update_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: fdm_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: tensor_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: schwarz_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: dong_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: coef_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: scalar_residual_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: lambda2_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: compute_max_wave_speed_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: mapping_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: find_rst_legendre_program = c_null_ptr
 
  type(c_ptr), public, bind(c) :: entropy_viscosity_program = c_null_ptr
 
  public :: opencl_prgm_lib_release
 
contains
 
  subroutine opencl_prgm_lib_release
 
    if (c_associated(math_program)) then
       if(clreleaseprogram(math_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       math_program = c_null_ptr
    end if
 
    if (c_associated(mathops_program)) then
       if(clreleaseprogram(mathops_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       mathops_program = c_null_ptr
    end if
 
    if (c_associated(dirichlet_program)) then
       if(clreleaseprogram(dirichlet_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       dirichlet_program = c_null_ptr
    end if
 
    if (c_associated(inflow_program)) then
       if(clreleaseprogram(inflow_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       inflow_program = c_null_ptr
    end if
 
    if (c_associated(zero_dirichlet_program)) then
       if (clreleaseprogram(zero_dirichlet_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       zero_dirichlet_program = c_null_ptr
    end if
 
    if (c_associated(symmetry_program)) then
       if(clreleaseprogram(symmetry_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       symmetry_program = c_null_ptr
    end if
 
    if (c_associated(facet_normal_program)) then
       if(clreleaseprogram(facet_normal_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       facet_normal_program = c_null_ptr
    end if
 
    if (c_associated(inhom_dirichlet_program)) then
       if(clreleaseprogram(inhom_dirichlet_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       inhom_dirichlet_program = c_null_ptr
    end if
 
    if (c_associated(dudxyz_program)) then
       if(clreleaseprogram(dudxyz_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       dudxyz_program = c_null_ptr
    end if
 
    if (c_associated(cdtp_program)) then
       if(clreleaseprogram(cdtp_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       cdtp_program = c_null_ptr
    end if
 
    if (c_associated(conv1_program)) then
       if(clreleaseprogram(conv1_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       conv1_program = c_null_ptr
    end if
 
    if (c_associated(cfl_program)) then
       if(clreleaseprogram(cfl_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       cfl_program = c_null_ptr
    end if
 
    if (c_associated(opgrad_program)) then
       if(clreleaseprogram(opgrad_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       opgrad_program = c_null_ptr
    end if
 
    if (c_associated(gs_program)) then
       if(clreleaseprogram(gs_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       gs_program = c_null_ptr
    end if
 
    if (c_associated(ax_helm_program)) then
       if(clreleaseprogram(ax_helm_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       ax_helm_program = c_null_ptr
    end if
 
    if (c_associated(ax_helm_full_program)) then
       if(clreleaseprogram(ax_helm_full_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       ax_helm_full_program = c_null_ptr
    end if
 
    if (c_associated(jacobi_program)) then
       if(clreleaseprogram(jacobi_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       jacobi_program = c_null_ptr
    end if
 
    if (c_associated(rhs_maker_program)) then
       if(clreleaseprogram(rhs_maker_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       rhs_maker_program = c_null_ptr
    end if
 
    if (c_associated(pnpn_res_program)) then
       if(clreleaseprogram(pnpn_res_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       pnpn_res_program = c_null_ptr
    end if
 
    if (c_associated(pnpn_stress_res_program)) then
       if(clreleaseprogram(pnpn_stress_res_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       pnpn_stress_res_program = c_null_ptr
    end if
 
    if (c_associated(euler_res_program)) then
       if(clreleaseprogram(euler_res_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       euler_res_program = c_null_ptr
    end if
 
    if (c_associated(compressible_ops_compute_max_wave_speed_program)) then
       if(clreleaseprogram(compressible_ops_compute_max_wave_speed_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       compressible_ops_compute_max_wave_speed_program = c_null_ptr
    end if
 
    if (c_associated(compressible_ops_compute_entropy_program)) then
       if(clreleaseprogram(compressible_ops_compute_entropy_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       compressible_ops_compute_entropy_program = c_null_ptr
    end if
 
    if (c_associated(compressible_ops_update_program)) then
       if(clreleaseprogram(compressible_ops_update_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       compressible_ops_update_program = c_null_ptr
    end if
 
    if (c_associated(fdm_program)) then
       if(clreleaseprogram(fdm_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       fdm_program = c_null_ptr
    end if
 
    if (c_associated(tensor_program)) then
       if(clreleaseprogram(tensor_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       tensor_program = c_null_ptr
    end if
 
    if (c_associated(schwarz_program)) then
       if(clreleaseprogram(schwarz_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       schwarz_program = c_null_ptr
    end if
 
    if (c_associated(dong_program)) then
       if(clreleaseprogram(dong_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       dong_program = c_null_ptr
    end if
 
    if (c_associated(coef_program)) then
       if(clreleaseprogram(coef_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       coef_program = c_null_ptr
    end if
 
    if (c_associated(scalar_residual_program)) then
       if(clreleaseprogram(scalar_residual_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       scalar_residual_program = c_null_ptr
    end if
 
    if (c_associated(lambda2_program)) then
       if(clreleaseprogram(lambda2_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       lambda2_program = c_null_ptr
    end if
 
    if (c_associated(compute_max_wave_speed_program)) then
       if(clreleaseprogram(compute_max_wave_speed_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       compute_max_wave_speed_program = c_null_ptr
    end if
 
    if (c_associated(mapping_program)) then
       if(clreleaseprogram(mapping_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       mapping_program = c_null_ptr
    end if
 
    if (c_associated(find_rst_legendre_program)) then
       if(clreleaseprogram(find_rst_legendre_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       find_rst_legendre_program = c_null_ptr
    end if
 
    if (c_associated(entropy_viscosity_program)) then
       if(clreleaseprogram(entropy_viscosity_program) .ne. cl_success) then
          call neko_error('Failed to release program')
       end if
       entropy_viscosity_program = c_null_ptr
    end if
 
  end subroutine opencl_prgm_lib_release
 
#endif
 
end module opencl_prgm_lib