linear_medium_to_em_field.F90

Go to the documentation of this file.

!! Copyright (C) 2021 F. Bonafé
!!
!! This program is free software; you can redistribute it and/or modify
!! it under the terms of the GNU General Public License as published by
!! the Free Software Foundation; either version 2, or (at your option)
!! any later version.
!!
!! This program is distributed in the hope that it will be useful,
!! but WITHOUT ANY WARRANTY; without even the implied warranty of
!! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!! GNU General Public License for more details.
!!
!! You should have received a copy of the GNU General Public License
!! along with this program; if not, write to the Free Software
!! Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
!! 02110-1301, USA.
!!
#include "global.h"
 
module linear_medium_to_em_field_oct_m
  use debug_oct_m
  use global_oct_m
  use grid_oct_m
  use interaction_oct_m
  use interaction_partner_oct_m
  use messages_oct_m
  use namespace_oct_m
  use namespace_oct_m
  use profiling_oct_m
  use quantity_oct_m
 
  implicit none
 
  private
  public ::                    &
    linear_medium_to_em_field_t,  &
    single_medium_box_t,       &
    single_medium_box_allocate,&
    single_medium_box_end
 
  type single_medium_box_t
    real(real64), allocatable            :: ep(:)
    real(real64), allocatable            :: mu(:)
    real(real64), allocatable            :: c(:)
    real(real64), allocatable            :: sigma_e(:)
    real(real64), allocatable            :: sigma_m(:)
    integer                       :: points_number
    integer                       :: global_points_number
    integer, allocatable          :: points_map(:)
    logical                       :: has_mapping = .true.
    integer                       :: bdry_number
    integer, allocatable          :: bdry_map(:)
    real(real64), allocatable            :: aux_ep(:,:)
    real(real64), allocatable            :: aux_mu(:,:)
    integer                       :: box_shape
    real(real64)                  :: center(3)
    real(real64)                  :: lsize(3)
    character(len=256)            :: filename
    logical                       :: check_medium_points = .false.
  contains
    procedure :: to_grid => single_medium_box_to_grid
  end type single_medium_box_t
 
  type, extends(interaction_t) :: linear_medium_to_em_field_t
    private
 
    type(grid_t), pointer, public    :: system_gr => null() 
    type(single_medium_box_t), public :: partner_medium_box
    type(single_medium_box_t), public :: medium_box
 
  contains
    procedure :: init => linear_medium_to_em_field_init
    procedure :: calculate => linear_medium_to_em_field_calculate
    procedure :: calculate_energy => linear_medium_to_em_field_calculate_energy
    final :: linear_medium_to_em_field_finalize
  end type linear_medium_to_em_field_t
 
 
  interface linear_medium_to_em_field_t
    module procedure linear_medium_to_em_field_constructor
  end interface linear_medium_to_em_field_t
 
contains
 
  ! ---------------------------------------------------------
  function linear_medium_to_em_field_constructor(partner) result(this)
    class(interaction_partner_t), target, intent(inout) :: partner
    class(linear_medium_to_em_field_t),      pointer       :: this
 
    push_sub(linear_medium_to_em_field_constructor)
 
    allocate(this)
 
    this%label = "linear_medium_to_em_field"
    this%partner => partner
 
    this%couplings_from_partner = [permittivity, permeability, e_conductivity, m_conductivity]
    this%intra_interaction = .false. ! This interaction does not support intra-interactions
 
    pop_sub(linear_medium_to_em_field_constructor)
  end function linear_medium_to_em_field_constructor
 
 
  subroutine linear_medium_to_em_field_init(this, gr)
    class(linear_medium_to_em_field_t), intent(inout) :: this
    type(grid_t), target, intent(in)         :: gr
 
    push_sub(linear_medium_to_em_field_init)
 
    this%system_gr => gr
    ! allocate medium box
    call single_medium_box_allocate(this%medium_box, gr%np)
    this%medium_box%has_mapping = .false.
 
    pop_sub(linear_medium_to_em_field_init)
  end subroutine linear_medium_to_em_field_init
 
  ! ---------------------------------------------------------
  subroutine linear_medium_to_em_field_finalize(this)
    type(linear_medium_to_em_field_t), intent(inout) :: this
 
    push_sub(linear_medium_to_em_field_finalize)
 
    pop_sub(linear_medium_to_em_field_finalize)
 
  end subroutine linear_medium_to_em_field_finalize
 
 
  ! ---------------------------------------------------------
  subroutine linear_medium_to_em_field_calculate(this)
    class(linear_medium_to_em_field_t), intent(inout) :: this
 
    push_sub(linear_medium_to_em_field_calculate)
 
    pop_sub(linear_medium_to_em_field_calculate)
  end subroutine linear_medium_to_em_field_calculate
 
  ! ---------------------------------------------------------
  subroutine linear_medium_to_em_field_calculate_energy(this)
    class(linear_medium_to_em_field_t), intent(inout) :: this
 
    push_sub(linear_medium_to_em_field_calculate_energy)
 
    this%energy = m_zero
 
    pop_sub(linear_medium_to_em_field_calculate_energy)
  end subroutine linear_medium_to_em_field_calculate_energy
 
 
  ! ---------------------------------------------------------
  subroutine single_medium_box_allocate(medium_box, n_points)
    type(single_medium_box_t),   intent(inout)    :: medium_box
    integer,                     intent(in)       :: n_points
 
 
    push_sub_with_profile(medium_box_allocate)
    safe_allocate(medium_box%aux_ep(1:n_points,1:3))
    safe_allocate(medium_box%aux_mu(1:n_points,1:3))
    safe_allocate(medium_box%c(1:n_points))
    safe_allocate(medium_box%ep(1:n_points))
    safe_allocate(medium_box%mu(1:n_points))
    safe_allocate(medium_box%sigma_e(1:n_points))
    safe_allocate(medium_box%sigma_m(1:n_points))
    safe_allocate(medium_box%points_map(1:n_points))
    medium_box%points_map = 0
    medium_box%aux_ep(:,1:3) = m_zero
    medium_box%aux_mu(:,1:3) = m_zero
    medium_box%ep(:) = m_zero
    medium_box%mu(:) = m_zero
    medium_box%c(:) = m_zero
    medium_box%sigma_e(:) = m_zero
    medium_box%sigma_m(:) = m_zero
    medium_box%points_number = n_points
    pop_sub_with_profile(medium_box_allocate)
 
  end subroutine single_medium_box_allocate
 
  ! ---------------------------------------------------------
  subroutine single_medium_box_end(medium_box)
    type(single_medium_box_t),   intent(inout)    :: medium_box
 
 
    push_sub_with_profile(medium_box_end)
 
    safe_deallocate_a(medium_box%points_map)
    safe_deallocate_a(medium_box%bdry_map)
    safe_deallocate_a(medium_box%aux_ep)
    safe_deallocate_a(medium_box%aux_mu)
    safe_deallocate_a(medium_box%c)
    safe_deallocate_a(medium_box%ep)
    safe_deallocate_a(medium_box%mu)
    safe_deallocate_a(medium_box%sigma_e)
    safe_deallocate_a(medium_box%sigma_m)
 
    pop_sub_with_profile(medium_box_end)
 
  end subroutine single_medium_box_end
 
  ! ---------------------------------------------------------
  ! return a medium box with all functions mapped to the parent grid
  function single_medium_box_to_grid(medium_box, grid_out) result(medium_box_out)
    class(single_medium_box_t), intent(in) :: medium_box
    type(grid_t),               intent(in) :: grid_out
    class(single_medium_box_t), pointer :: medium_box_out
 
    integer :: ip, ip_out
 
    push_sub_with_profile(medium_box_to_grid)
 
    safe_allocate(medium_box_out)
    call single_medium_box_allocate(medium_box_out, grid_out%np)
    medium_box_out%has_mapping = .false.
 
    do ip = 1, medium_box%points_number
      ip_out = medium_box%points_map(ip)
      medium_box_out%aux_ep(ip_out,1:3) = medium_box%aux_ep(ip,1:3)
      medium_box_out%aux_mu(ip_out,1:3) = medium_box%aux_mu(ip,1:3)
      medium_box_out%ep(ip_out) = medium_box%ep(ip)
      medium_box_out%mu(ip_out) = medium_box%mu(ip)
      medium_box_out%c(ip_out) = medium_box%c(ip)
      medium_box_out%sigma_e(ip_out) = medium_box%sigma_e(ip)
      medium_box_out%sigma_m(ip_out) = medium_box%sigma_m(ip)
    end do
 
    pop_sub_with_profile(medium_box_to_grid)
  end function single_medium_box_to_grid
 
end module linear_medium_to_em_field_oct_m
 
!! Local Variables:
!! mode: f90
!! coding: utf-8
!! End: