find_target_edges_by_group Module Subroutine

module subroutine find_target_edges_by_group(domain, i_material, target_edges)

Arguments

Type IntentOptional Attributes Name
type(type_domain), intent(in) :: domain
integer(kind=int32), intent(in) :: i_material
integer(kind=int32), intent(inout), allocatable :: target_edges(:,:)

Calls

proc~~find_target_edges_by_group~~CallsGraph proc~find_target_edges_by_group find_target_edges_by_group get_group get_group proc~find_target_edges_by_group->get_group get_order get_order proc~find_target_edges_by_group->get_order interface~allocate_array allocate_array proc~find_target_edges_by_group->interface~allocate_array log_warning log_warning proc~find_target_edges_by_group->log_warning proc~get_num_sides type_domain%get_num_sides proc~find_target_edges_by_group->proc~get_num_sides to_string to_string proc~find_target_edges_by_group->to_string proc~allocate_rank1_int16 allocate_rank1_int16 interface~allocate_array->proc~allocate_rank1_int16 proc~allocate_rank1_int32 allocate_rank1_int32 interface~allocate_array->proc~allocate_rank1_int32 proc~allocate_rank1_int64 allocate_rank1_int64 interface~allocate_array->proc~allocate_rank1_int64 proc~allocate_rank1_int8 allocate_rank1_int8 interface~allocate_array->proc~allocate_rank1_int8 proc~allocate_rank1_logical1 allocate_rank1_logical1 interface~allocate_array->proc~allocate_rank1_logical1 proc~allocate_rank1_logical4 allocate_rank1_logical4 interface~allocate_array->proc~allocate_rank1_logical4 proc~allocate_rank1_logical8 allocate_rank1_logical8 interface~allocate_array->proc~allocate_rank1_logical8 proc~allocate_rank1_real128 allocate_rank1_real128 interface~allocate_array->proc~allocate_rank1_real128 proc~allocate_rank1_real32 allocate_rank1_real32 interface~allocate_array->proc~allocate_rank1_real32 proc~allocate_rank1_real64 allocate_rank1_real64 interface~allocate_array->proc~allocate_rank1_real64 proc~allocate_rank2_int16 allocate_rank2_int16 interface~allocate_array->proc~allocate_rank2_int16 proc~allocate_rank2_int32 allocate_rank2_int32 interface~allocate_array->proc~allocate_rank2_int32 proc~allocate_rank2_int64 allocate_rank2_int64 interface~allocate_array->proc~allocate_rank2_int64 proc~allocate_rank2_int8 allocate_rank2_int8 interface~allocate_array->proc~allocate_rank2_int8 proc~allocate_rank2_logical1 allocate_rank2_logical1 interface~allocate_array->proc~allocate_rank2_logical1 proc~allocate_rank2_logical4 allocate_rank2_logical4 interface~allocate_array->proc~allocate_rank2_logical4 proc~allocate_rank2_logical8 allocate_rank2_logical8 interface~allocate_array->proc~allocate_rank2_logical8 proc~allocate_rank2_real128 allocate_rank2_real128 interface~allocate_array->proc~allocate_rank2_real128 proc~allocate_rank2_real32 allocate_rank2_real32 interface~allocate_array->proc~allocate_rank2_real32 proc~allocate_rank2_real64 allocate_rank2_real64 interface~allocate_array->proc~allocate_rank2_real64 proc~error_message error_message proc~allocate_rank1_int16->proc~error_message proc~allocate_rank1_int32->proc~error_message proc~allocate_rank1_int64->proc~error_message proc~allocate_rank1_int8->proc~error_message proc~allocate_rank1_logical1->proc~error_message proc~allocate_rank1_logical4->proc~error_message proc~allocate_rank1_logical8->proc~error_message proc~allocate_rank1_real128->proc~error_message proc~allocate_rank1_real32->proc~error_message proc~allocate_rank1_real64->proc~error_message proc~allocate_rank2_int16->proc~error_message proc~allocate_rank2_int32->proc~error_message proc~allocate_rank2_int64->proc~error_message proc~allocate_rank2_int8->proc~error_message proc~allocate_rank2_logical1->proc~error_message proc~allocate_rank2_logical4->proc~error_message proc~allocate_rank2_logical8->proc~error_message proc~allocate_rank2_real128->proc~error_message proc~allocate_rank2_real32->proc~error_message proc~allocate_rank2_real64->proc~error_message log_error log_error proc~error_message->log_error

Called by

proc~~find_target_edges_by_group~~CalledByGraph proc~find_target_edges_by_group find_target_edges_by_group interface~find_target_edges_by_group find_target_edges_by_group interface~find_target_edges_by_group->proc~find_target_edges_by_group proc~initialize_type_bc_thermal_adiabatic initialize_type_bc_thermal_adiabatic proc~initialize_type_bc_thermal_adiabatic->interface~find_target_edges_by_group proc~initialize_type_bc_thermal_dirichlet initialize_type_bc_thermal_dirichlet proc~initialize_type_bc_thermal_dirichlet->interface~find_target_edges_by_group interface~initialize_type_bc_thermal_adiabatic type_bc_thermal_adiabatic%initialize_type_bc_thermal_adiabatic interface~initialize_type_bc_thermal_adiabatic->proc~initialize_type_bc_thermal_adiabatic interface~initialize_type_bc_thermal_dirichlet type_bc_thermal_dirichlet%initialize_type_bc_thermal_dirichlet interface~initialize_type_bc_thermal_dirichlet->proc~initialize_type_bc_thermal_dirichlet

Source Code

    module subroutine find_target_edges_by_group(domain, i_material, target_edges)
        implicit none
        type(type_domain), intent(in) :: domain
        integer(int32), intent(in) :: i_material
        integer(int32), allocatable, intent(inout) :: target_edges(:, :)

        ! Local variables
        integer(int32) :: i, num_total_segments, current_segment_idx
        integer(int32) :: side_group, num_sides, side_order

        ! 1. Count segments
        num_total_segments = 0
        num_sides = domain%get_num_sides()
        do i = 1, num_sides
            if (allocated(domain%Sides(i)%s)) then
                side_group = domain%Sides(i)%s%get_group()
                if (side_group == i_material) then
                    side_order = domain%Sides(i)%s%get_order()
                    if (side_order == 1) then
                        num_total_segments = num_total_segments + 1
                    elseif (side_order == 2) then
                        num_total_segments = num_total_segments + 2
                    end if
                end if
            end if
        end do

        if (num_total_segments == 0) then
            if (allocated(target_edges)) deallocate (target_edges)
            call global_logger%log_warning(message="No target edges found for material ID: "//to_string(i_material))
            return
        end if

        ! 2. Allocate and fill the array
        call Allocate_Array(target_edges, 2_int32, num_total_segments)
        current_segment_idx = 0
        do i = 1, num_sides
            if (allocated(domain%Sides(i)%s)) then
                side_group = domain%Sides(i)%s%get_group()
                if (side_group == i_material) then
                    side_order = domain%Sides(i)%s%get_order()
                    if (side_order == 1) then
                        current_segment_idx = current_segment_idx + 1
                        target_edges(:, current_segment_idx) = domain%Sides(i)%s%connectivity([1, 2])
                    elseif (side_order == 2) then
                        current_segment_idx = current_segment_idx + 1
                        target_edges(:, current_segment_idx) = domain%Sides(i)%s%connectivity([1, 3])
                        current_segment_idx = current_segment_idx + 1
                        target_edges(:, current_segment_idx) = domain%Sides(i)%s%connectivity([3, 2])
                    end if
                end if
            end if
        end do

    end subroutine find_target_edges_by_group