sll_m_time_propagator_pic_vm_3d3v_cef_trafo Module Reference

Description

Particle pusher based on Hamiltonian splitting for 3d3v Vlasov-Maxwell with coordinate transformation.

Author

Benedikt Perse, IPP

MPI parallelization by domain cloning. Periodic boundaries. Spline DoFs numerated by the point the spline starts.

Derived types and interfaces
type	sll_t_time_propagator_pic_vm_3d3v_cef_trafo
	Hamiltonian splitting type for Vlasov-Maxwell 3d3v. More...

Functions/Subroutines
subroutine	strang_splitting_pic_vm_3d3v_cef_trafo (self, dt, number_steps)
	Finalization. More...
subroutine	lie_splitting_pic_vm_3d3v_cef_trafo (self, dt, number_steps)
	Lie splitting. More...
subroutine	lie_splitting_back_pic_vm_3d3v_cef_trafo (self, dt, number_steps)
	Lie splitting (oposite ordering) More...
subroutine	operatorhp_pic_vm_3d3v_cef_trafo (self, dt)
	Push H_p: Equations to be solved $\Xi^{n+1}=\Xi^n+\Delta t DF^{-1}(\bar{\Xi}) V^n$ $M_1 e^n= M_1 e^n- \int_{t^n}^{t^{n+1}} \mathbb{\Lambda}^1(\Xi(\tau))^\top d\tau \mathbb{W}_q DF^{-1}(\bar{\Xi}) V^n$. More...
subroutine	operatorhb_pic_vm_3d3v_cef_trafo (self, dt)
	Push H_B: Equations to be solved $(\mathbb{I}-\Delta \frac{\Delta t q}{2 m} DF^{-\top} \mathbb{B}(\Xi^n,b^n) DF^{-1}) V^{n+1}=(\mathbb{I}+ \frac{\Delta t q}{2 m} DF^{-\top} \mathbb{B}(\Xi^n,b^n) DF^{-1}) V^n$ $M_1 e^{n+1}=M_1 e^n+\Delta t C^\top M_2 b^n$. More...
subroutine	operatorhe_pic_vm_3d3v_cef_trafo (self, dt)
	Push H_E: Equations to be solved $V^{n+1}=V^n+\Delta t\mathbb{W}_{\frac{q}{m}} DF^{-\top}(\Xi^n) \mathbb{Lambda}^1(\Xi^n) e^n$ $b^{n+1}=b^n-\Delta t C e^n$. More...
subroutine	initialize_pic_vm_3d3v_cef_trafo (self, maxwell_solver, particle_mesh_coupling, particle_group, phi_dofs, efield_dofs, bfield_dofs, x_min, Lx, map, boundary_particles, iter_tolerance, max_iter, betar, electrostatic, rhob, control_variate)
	Constructor. More...
subroutine	initialize_file_pic_vm_3d3v_cef_trafo (self, maxwell_solver, particle_mesh_coupling, particle_group, phi_dofs, efield_dofs, bfield_dofs, x_min, Lx, map, filename, boundary_particles, betar, electrostatic, rhob, control_variate)
	Constructor. More...
subroutine	delete_pic_vm_3d3v_cef_trafo (self)
	Destructor. More...

Function/Subroutine Documentation

delete_pic_vm_3d3v_cef_trafo()

subroutine sll_m_time_propagator_pic_vm_3d3v_cef_trafo::delete_pic_vm_3d3v_cef_trafo ( class(sll_t_time_propagator_pic_vm_3d3v_cef_trafo), intent(inout)  self )

private

Destructor.

Parameters

[in,out]

self

time propagator object

Definition at line 687 of file sll_m_time_propagator_pic_vm_3d3v_cef_trafo.F90.

initialize_file_pic_vm_3d3v_cef_trafo()

subroutine sll_m_time_propagator_pic_vm_3d3v_cef_trafo::initialize_file_pic_vm_3d3v_cef_trafo ( class(sll_t_time_propagator_pic_vm_3d3v_cef_trafo), intent(out)  self, class(sll_c_maxwell_3d_base), intent(in), target  maxwell_solver, class(sll_c_particle_mesh_coupling_3d), intent(in), target  particle_mesh_coupling, class(sll_t_particle_array), intent(in), target  particle_group, real(kind=f64), dimension(:), intent(in), target  phi_dofs, real(kind=f64), dimension(:), intent(in), target  efield_dofs, real(kind=f64), dimension(:), intent(in), target  bfield_dofs, real(kind=f64), dimension(3), intent(in)  x_min, real(kind=f64), dimension(3), intent(in)  Lx, type(sll_t_mapping_3d), intent(inout), target  map, character(len=*), intent(in)  filename, integer(kind=i32), intent(in), optional  boundary_particles, real(kind=f64), dimension(2), intent(in), optional  betar, logical, optional  electrostatic, real(kind=f64), dimension(:), intent(in), optional, target  rhob, class(sll_t_control_variates), intent(in), optional, target  control_variate  )

private

Constructor.

Parameters

[out]	self	time propagator object
[in]	maxwell_solver	Maxwell solver
[in]	particle_mesh_coupling	Particle mesh coupling
[in]	particle_group	Particle group
[in]	phi_dofs	array for the coefficients of the scalar potential
[in]	efield_dofs	array for the coefficients of the efields
[in]	bfield_dofs	array for the coefficients of the bfield
[in]	x_min	Lower bound of x domain
[in]	lx	Length of the domain in x direction.
[in,out]	map	Coordinate transformation
[in]	filename	file name
[in]	betar	reciprocal plasma beta
	electrostatic	true for electrostatic simulation
[in]	rhob	charge at the boundary
[in]	control_variate	Control variate (if delta f)

Definition at line 518 of file sll_m_time_propagator_pic_vm_3d3v_cef_trafo.F90.

initialize_pic_vm_3d3v_cef_trafo()

subroutine sll_m_time_propagator_pic_vm_3d3v_cef_trafo::initialize_pic_vm_3d3v_cef_trafo ( class(sll_t_time_propagator_pic_vm_3d3v_cef_trafo), intent(out)  self, class(sll_c_maxwell_3d_base), intent(in), target  maxwell_solver, class(sll_c_particle_mesh_coupling_3d), intent(in), target  particle_mesh_coupling, class(sll_t_particle_array), intent(in), target  particle_group, real(kind=f64), dimension(:), intent(in), target  phi_dofs, real(kind=f64), dimension(:), intent(in), target  efield_dofs, real(kind=f64), dimension(:), intent(in), target  bfield_dofs, real(kind=f64), dimension(3), intent(in)  x_min, real(kind=f64), dimension(3), intent(in)  Lx, type(sll_t_mapping_3d), intent(inout), target  map, integer(kind=i32), intent(in), optional  boundary_particles, real(kind=f64), intent(in), optional  iter_tolerance, integer(kind=i32), intent(in), optional  max_iter, real(kind=f64), dimension(2), intent(in), optional  betar, logical, optional  electrostatic, real(kind=f64), dimension(:), intent(in), optional, target  rhob, class(sll_t_control_variates), intent(in), optional, target  control_variate  )

private

Constructor.

Parameters

[out]	self	time propagator object
[in]	maxwell_solver	Maxwell solver
[in]	particle_mesh_coupling	Particle mesh coupling
[in]	particle_group	Particle group
[in]	phi_dofs	array for the coefficients of the scalar potential
[in]	efield_dofs	array for the coefficients of the efields
[in]	bfield_dofs	array for the coefficients of the bfield
[in]	x_min	Lower bound of x domain
[in]	lx	Length of the domain in x direction.
[in,out]	map	Coordinate transformation
[in]	boundary_particles	particle boundary conditions
[in]	iter_tolerance	iteration tolerance
[in]	max_iter	maximal number of iterations
[in]	betar	reciprocal plasma beta
[in]	rhob	charge at the boundary
[in]	control_variate	Control variate (if delta f)

Definition at line 421 of file sll_m_time_propagator_pic_vm_3d3v_cef_trafo.F90.

lie_splitting_back_pic_vm_3d3v_cef_trafo()

subroutine sll_m_time_propagator_pic_vm_3d3v_cef_trafo::lie_splitting_back_pic_vm_3d3v_cef_trafo ( class(sll_t_time_propagator_pic_vm_3d3v_cef_trafo), intent(inout)  self, real(kind=f64), intent(in)  dt, integer(kind=i32), intent(in)  number_steps  )

private

Lie splitting (oposite ordering)

Parameters

[in,out]	self	time propagator object
[in]	dt	time step
[in]	number_steps	number of time steps

Definition at line 158 of file sll_m_time_propagator_pic_vm_3d3v_cef_trafo.F90.

lie_splitting_pic_vm_3d3v_cef_trafo()

subroutine sll_m_time_propagator_pic_vm_3d3v_cef_trafo::lie_splitting_pic_vm_3d3v_cef_trafo ( class(sll_t_time_propagator_pic_vm_3d3v_cef_trafo), intent(inout)  self, real(kind=f64), intent(in)  dt, integer(kind=i32), intent(in)  number_steps  )

private

Lie splitting.

Parameters

[in,out]	self	time propagator object
[in]	dt	time step
[in]	number_steps	number of time steps

Definition at line 141 of file sll_m_time_propagator_pic_vm_3d3v_cef_trafo.F90.

operatorhb_pic_vm_3d3v_cef_trafo()

subroutine sll_m_time_propagator_pic_vm_3d3v_cef_trafo::operatorhb_pic_vm_3d3v_cef_trafo ( class(sll_t_time_propagator_pic_vm_3d3v_cef_trafo), intent(inout)  self, real(kind=f64), intent(in)  dt  )

private

Push H_B: Equations to be solved $(\mathbb{I}-\Delta \frac{\Delta t q}{2 m} DF^{-\top} \mathbb{B}(\Xi^n,b^n) DF^{-1}) V^{n+1}=(\mathbb{I}+ \frac{\Delta t q}{2 m} DF^{-\top} \mathbb{B}(\Xi^n,b^n) DF^{-1}) V^n$ $M_1 e^{n+1}=M_1 e^n+\Delta t C^\top M_2 b^n$.

Parameters

[in,out]	self	time propagator object
[in]	dt	time step

Definition at line 266 of file sll_m_time_propagator_pic_vm_3d3v_cef_trafo.F90.

Here is the call graph for this function:

operatorhe_pic_vm_3d3v_cef_trafo()

subroutine sll_m_time_propagator_pic_vm_3d3v_cef_trafo::operatorhe_pic_vm_3d3v_cef_trafo ( class(sll_t_time_propagator_pic_vm_3d3v_cef_trafo), intent(inout)  self, real(kind=f64), intent(in)  dt  )

private

Push H_E: Equations to be solved $V^{n+1}=V^n+\Delta t\mathbb{W}_{\frac{q}{m}} DF^{-\top}(\Xi^n) \mathbb{Lambda}^1(\Xi^n) e^n$ $b^{n+1}=b^n-\Delta t C e^n$.

Parameters

[in,out]	self	time propagator object
[in]	dt	time step

Definition at line 337 of file sll_m_time_propagator_pic_vm_3d3v_cef_trafo.F90.

operatorhp_pic_vm_3d3v_cef_trafo()

subroutine sll_m_time_propagator_pic_vm_3d3v_cef_trafo::operatorhp_pic_vm_3d3v_cef_trafo ( class(sll_t_time_propagator_pic_vm_3d3v_cef_trafo), intent(inout)  self, real(kind=f64), intent(in)  dt  )

private

Push H_p: Equations to be solved $\Xi^{n+1}=\Xi^n+\Delta t DF^{-1}(\bar{\Xi}) V^n$ $M_1 e^n= M_1 e^n- \int_{t^n}^{t^{n+1}} \mathbb{\Lambda}^1(\Xi(\tau))^\top d\tau \mathbb{W}_q DF^{-1}(\bar{\Xi}) V^n$.

Parameters

[in,out]	self	time propagator object
[in]	dt	time step

Definition at line 178 of file sll_m_time_propagator_pic_vm_3d3v_cef_trafo.F90.

Here is the call graph for this function:

strang_splitting_pic_vm_3d3v_cef_trafo()

subroutine sll_m_time_propagator_pic_vm_3d3v_cef_trafo::strang_splitting_pic_vm_3d3v_cef_trafo ( class(sll_t_time_propagator_pic_vm_3d3v_cef_trafo), intent(inout)  self, real(kind=f64), intent(in)  dt, integer(kind=i32), intent(in)  number_steps  )

private

Finalization.

Strang splitting

Parameters

[in,out]	self	time propagator object
[in]	dt	time step
[in]	number_steps	number of time steps

Definition at line 122 of file sll_m_time_propagator_pic_vm_3d3v_cef_trafo.F90.