sll_m_particle_mesh_coupling_base_1d.F90

Go to the documentation of this file.

 
module sll_m_particle_mesh_coupling_base_1d
 
  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#include "sll_working_precision.h"
 
  implicit none
 
  public :: &
       sll_p_collocation, &
       sll_p_galerkin, &
       sll_c_particle_mesh_coupling_1d
 
  private
  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
  ! Define parameters to set if Galerkin or collocation scaling should be used in accumulation routines
  sll_int32, parameter :: sll_p_galerkin = 0
  sll_int32, parameter :: sll_p_collocation = 1
 
  type, abstract :: sll_c_particle_mesh_coupling_1d
     sll_int32              :: dim
     sll_real64             :: delta_x 
     sll_int32              :: n_cells 
     sll_int32              :: n_dofs  
     sll_real64             :: domain(2) 
     sll_int32              :: spline_degree 
 
   contains
 
     procedure(add_single), deferred        :: add_charge
     procedure(add_array), deferred         :: add_particle_mass
     procedure(add_array), deferred         :: add_particle_mass_full
     procedure(eval_single), deferred       :: evaluate
     procedure(eval_multiple), deferred     :: evaluate_multiple
     procedure(add_update) , deferred       :: add_current_update_v
     procedure(add_current), deferred       :: add_current
     procedure(add_current_evaluate), deferred :: add_current_evaluate
     procedure(empty), deferred             :: free
 
 
  end type sll_c_particle_mesh_coupling_1d
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine add_single(self, position, marker_charge, rho_dofs) 
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_1d
       class(sll_c_particle_mesh_coupling_1d), intent( inout ) :: self 
       sll_real64,                    intent( in )    :: position(self%dim) 
       sll_real64,                    intent( in )    :: marker_charge 
       sll_real64,                    intent( inout ) :: rho_dofs(self%n_dofs) 
 
     end subroutine add_single
  end interface
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine add_array(self, position, marker_charge, particle_mass) 
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_1d
       class(sll_c_particle_mesh_coupling_1d), intent( inout ) :: self 
       sll_real64,                    intent( in )    :: position(self%dim) 
       sll_real64,                    intent( in )    :: marker_charge 
       sll_real64,                    intent( inout ) :: particle_mass(:, :) 
 
     end subroutine add_array
  end interface
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine eval_single(self, position, field_dofs, field_value)
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_1d
       class(sll_c_particle_mesh_coupling_1d), intent( inout ) :: self 
       sll_real64,                    intent( in )    :: position(self%dim) 
       sll_real64,                    intent( in )    :: field_dofs(self%n_dofs) 
       sll_real64, intent(out) :: field_value 
     end subroutine eval_single
  end interface
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine eval_multiple(self, position, components, field_dofs, field_value)
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_1d
       class(sll_c_particle_mesh_coupling_1d), intent( inout ) :: self 
       sll_real64,                    intent( in )    :: position(self%dim) 
       sll_int32,                     intent(in)      :: components(:)
       sll_real64,                    intent( in )    :: field_dofs(:,:) 
       sll_real64,                    intent(out)     :: field_value(:) 
     end subroutine eval_multiple
  end interface
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine add_current(self, position_old, position_new, marker_charge, j_dofs) 
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_1d
       class(sll_c_particle_mesh_coupling_1d), intent( inout ) :: self 
       sll_real64,                    intent( in )    :: position_old(self%dim) 
       sll_real64,                    intent( in )    :: position_new(self%dim) 
       sll_real64,                    intent( in )    :: marker_charge 
       sll_real64,                    intent( inout ) :: j_dofs(self%n_dofs) 
 
     end subroutine add_current
  end interface
 
   !---------------------------------------------------------------------------!
  abstract interface
     subroutine add_current_evaluate(self, position_old, position_new, marker_charge, vbar, field_dofs, j_dofs, field)
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_1d
       class(sll_c_particle_mesh_coupling_1d), intent( inout ) :: self 
       sll_real64,                    intent( in )    :: position_old(self%dim) 
       sll_real64,                    intent( in )    :: position_new(self%dim) 
       sll_real64,                    intent( in )    :: marker_charge 
       sll_real64,                               intent( in )    :: vbar 
       sll_real64,                               intent( in    ) :: field_dofs(self%n_dofs) 
       sll_real64,                    intent( inout ) :: j_dofs(self%n_dofs) 
       sll_real64,                               intent( out   ) :: field 
       
     end subroutine add_current_evaluate
  end interface
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine add_update (self, &
          position_old, &
          position_new, &
          marker_charge, &
          qoverm, &
          bfield_dofs, &
          vi, &
          j_dofs)
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_1d
       class(sll_c_particle_mesh_coupling_1d), intent(inout) :: self
       sll_real64, intent(in)    :: position_old(self%dim) 
       sll_real64, intent(in)    :: position_new(self%dim) 
       sll_real64, intent(in)    :: marker_charge          
       sll_real64, intent(in)    :: qoverm   
       sll_real64, intent(in)    :: bfield_dofs(:) 
       sll_real64, intent(inout) :: vi(:) 
       sll_real64, intent(inout) :: j_dofs(:) 
 
     end subroutine add_update
  end interface
 
  !---------------------------------------------------------------------------!  
  abstract interface
     subroutine empty(self)
       import sll_c_particle_mesh_coupling_1d
       class(sll_c_particle_mesh_coupling_1d), intent( inout ) :: self 
 
     end subroutine empty
  end interface
 
end module sll_m_particle_mesh_coupling_base_1d