mapping_fctry.f90

 
!
submodule(mapping) mapping_fctry
  use utils, only: neko_type_error
  use linear_mapping, only : linear_mapping_t
  use pde_filter, only: pde_filter_t
  use ramp_mapping, only: ramp_mapping_t
  use heaviside_mapping, only: heaviside_mapping_t
  use simp_mapping, only: simp_mapping_t
  use json_utils, only : json_get
  use utils, only : concat_string_array, neko_error
  implicit none
 
  ! List of all possible types created by the factory routine
  character(len=20) :: MAPPING_KNOWN_TYPES(5) = [character(len=20) :: &
       "linear", &
       "PDE_filter", &
       "RAMP", &
       "SIMP", &
       "heaviside_mapping"]
 
contains
 
  module subroutine mapping_factory(object, json, coef)
    class(mapping_t), allocatable, intent(inout) :: object
    type(json_file), intent(inout) :: json
    type(coef_t), intent(inout) :: coef
    character(len=:), allocatable :: type_name
 
    call json_get(json, "type", type_name)
 
    select case (trim(type_name))
    case ("linear")
       allocate(linear_mapping_t::object)
    case ("PDE_filter")
       allocate(pde_filter_t::object)
    case ("RAMP")
       allocate(ramp_mapping_t::object)
    case ("SIMP")
       allocate(simp_mapping_t::object)
    case ("heaviside_mapping", "heaviside_projection")
       allocate(heaviside_mapping_t::object)
    case default
       call neko_type_error("Mapping function", type_name, mapping_known_types)
    end select
 
    ! Initialize
    call object%init(json, coef)
 
  end subroutine mapping_factory
 
end submodule mapping_fctry