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