sll_m_particle_mesh_coupling_base_3d.F90

Go to the documentation of this file.

 
module sll_m_particle_mesh_coupling_base_3d
 
  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#include "sll_working_precision.h"
 
  use sll_m_splines_pp, only : &
       sll_t_spline_pp_3d
 
  implicit none
 
  public :: &
       sll_c_particle_mesh_coupling_3d
 
  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_3d
     sll_int32                :: dim = 3
     sll_real64               :: delta_x(3)  
     sll_real64               :: domain(3,2) 
     sll_int32                :: n_total  
     sll_int32                :: spline_degree(3) 
     sll_int32                :: n_cells(3) 
     sll_int32                :: n_dofs(3) 
 
     type(sll_t_spline_pp_3d) :: spline_pp_0
     type(sll_t_spline_pp_3d) :: spline_pp_1
 
   contains
 
     procedure(add_single), deferred           :: add_charge
     procedure(add_array), deferred            :: add_particle_mass
     procedure(add_array_mixed), deferred      :: add_particle_mass_od
     procedure(eval_single), deferred          :: evaluate
     procedure(eval_multiple), deferred        :: evaluate_multiple
     procedure(add_current), deferred          :: add_current
     procedure(add_current_evaluate), deferred :: add_current_evaluate
     procedure(add_update), deferred           :: add_current_update_v_component1
     procedure(add_update), deferred           :: add_current_update_v_component2
     procedure(add_update), deferred           :: add_current_update_v_component3
     procedure(empty), deferred                :: free
 
 
  end type sll_c_particle_mesh_coupling_3d
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine add_single(self, position, marker_charge, degree, rho_dofs)
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_3d
       class(sll_c_particle_mesh_coupling_3d), intent( inout ) :: self 
       sll_real64,                               intent( in )    :: position(3) 
       sll_real64,                               intent( in )    :: marker_charge 
       sll_int32,                                intent( in )    :: degree(3) 
       sll_real64,                               intent( inout ) :: rho_dofs(:)
 
     end subroutine add_single
  end interface
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine add_array(self, position, marker_charge, degree, particle_mass )
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_3d
       class(sll_c_particle_mesh_coupling_3d), intent( inout ) :: self 
       sll_real64,                               intent( in )    :: position(3) 
       sll_real64,                               intent( in )    :: marker_charge 
       sll_int32,                                intent( in )    :: degree(3) 
       sll_real64,                    intent( inout ) :: particle_mass(:, :) 
 
     end subroutine add_array
  end interface
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine add_array_mixed(self, position, marker_charge, degree1, degree2, particle_mass )
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_3d
       class(sll_c_particle_mesh_coupling_3d), intent( inout ) :: self 
       sll_real64,                               intent( in )    :: position(3) 
       sll_real64,                               intent( in )    :: marker_charge 
       sll_int32,                                intent( in )    :: degree1(3), degree2(3) 
       sll_real64,                    intent( inout ) :: particle_mass(:, :) 
 
     end subroutine add_array_mixed
  end interface
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine eval_single(self, position, degree, field_dofs, field_value)
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_3d
       class(sll_c_particle_mesh_coupling_3d), intent( inout ) :: self 
       sll_real64,                              intent( in )    :: position(3) 
       sll_int32 ,                              intent( in )    :: degree(3) 
       sll_real64,                              intent( in )    :: field_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_3d
       class(sll_c_particle_mesh_coupling_3d), intent( inout ) :: self 
       sll_real64,                    intent( in )    :: position(3) 
       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, xdot, j_dofs) 
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_3d
       class(sll_c_particle_mesh_coupling_3d), intent( inout ) :: self 
       sll_real64,                    intent( in )    :: position_old(3) 
       sll_real64,                    intent( in )    :: position_new(3) 
       sll_real64,                    intent( in )    :: xdot(3) 
       sll_real64,                    intent( inout ) :: j_dofs(:) 
 
     end subroutine add_current
  end interface
 
  !---------------------------------------------------------------------------!
  abstract interface
     subroutine add_current_evaluate(self, position_old, position_new, xdot, efield_dofs, j_dofs, efield_val) 
       use sll_m_working_precision
       import sll_c_particle_mesh_coupling_3d
       class(sll_c_particle_mesh_coupling_3d), intent( inout ) :: self 
       sll_real64,                               intent( in )    :: position_old(3) 
       sll_real64,                               intent( in )    :: position_new(3) 
       sll_real64,                               intent( in )    :: xdot(3) 
       sll_real64,                               intent( in )    :: efield_dofs(:) 
       sll_real64,                               intent( inout ) :: j_dofs(:) 
       sll_real64,                               intent( out   ) :: efield_val(3) 
 
     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_3d
       class(sll_c_particle_mesh_coupling_3d), intent(inout) :: self
       sll_real64, intent(in)    :: position_old(3) 
       sll_real64, intent(in)    :: position_new 
       sll_real64, intent(in)    :: marker_charge          
       sll_real64, intent(in)    :: qoverm   
       sll_real64, intent(in)    :: bfield_dofs(:) 
       sll_real64, intent(inout) :: vi(3) 
       sll_real64, intent(inout) :: j_dofs(:) 
 
     end subroutine add_update
  end interface
 
  !---------------------------------------------------------------------------!  
  abstract interface
     subroutine empty(self)
       import sll_c_particle_mesh_coupling_3d
       class(sll_c_particle_mesh_coupling_3d), intent( inout ) :: self 
 
     end subroutine empty
  end interface
 
end module sll_m_particle_mesh_coupling_base_3d