localreduce-operator.F90

Go to the documentation of this file.

module localreduce_operator_mod
  use datadefn_mod, only : default_precision, string_length
  use configuration_parser_mod, only : io_configuration_type, data_values_type, get_data_value_by_field_name
  use data_utils_mod, only : get_action_attribute_string
  use collections_mod, only : hashmap_type, list_type, map_type, c_add_string  
  implicit none

#ifndef TEST_MODE
  private
#endif

  public perform_localreduce_operator, localreduce_operator_get_required_fields
contains

  subroutine perform_localreduce_operator(io_configuration, field_values, action_attributes, source_monc_location, &
       source_monc, operator_result_values)
    type(io_configuration_type), intent(inout) :: io_configuration
    type(hashmap_type), intent(inout) :: field_values
    type(map_type), intent(inout) :: action_attributes
    integer, intent(in) :: source_monc_location, source_monc
    real(kind=DEFAULT_PRECISION), dimension(:), allocatable, intent(inout) :: operator_result_values

    character(len=STRING_LENGTH) :: field_to_reduce, reduction_operator
    type(data_values_type), pointer :: field_local_values

    allocate(operator_result_values(1))

    field_to_reduce=get_action_attribute_string(action_attributes, "field")
    reduction_operator=get_action_attribute_string(action_attributes, "operator")
    field_local_values=>get_data_value_by_field_name(field_values, field_to_reduce)
    operator_result_values=do_local_reduction(field_local_values%values, reduction_operator)
  end subroutine perform_localreduce_operator

  real(kind=DEFAULT_PRECISION) function do_local_reduction(data, reduction_operator)
    real(kind=DEFAULT_PRECISION), dimension(:), intent(in) :: data
    character(len=STRING_LENGTH), intent(in) :: reduction_operator

    if (reduction_operator .eq. "max") then
      do_local_reduction=maxval(data)
    else if (reduction_operator .eq. "min") then
      do_local_reduction=minval(data)
    else if (reduction_operator .eq. "sum") then
      do_local_reduction=sum(data)
    else if (reduction_operator .eq. "mean") then
      do_local_reduction=sum(data)/size(data)
    end if    
  end function do_local_reduction  

  type(list_type) function localreduce_operator_get_required_fields(action_attributes)
    type(map_type), intent(inout) :: action_attributes

    character(len=STRING_LENGTH) :: field_to_reduce

    field_to_reduce=get_action_attribute_string(action_attributes, "field")
    call c_add_string(localreduce_operator_get_required_fields, field_to_reduce)
  end function localreduce_operator_get_required_fields
end module localreduce_operator_mod