#include "globals.hpp"
#include "magnetic_field.hpp"
#include "grid.hpp"
#include "particles.hpp"
#include "species.hpp"
#include "perturbed_B_field.hpp"
#include "grid_field_pack.hpp"
#include "field_aligned_local_fields.hpp"
#include "loop_voltage.hpp"

Include dependency graph for push_physics.tpp:

This graph shows which files directly or indirectly include this file:

Functions
template<class Device , KinType KT, PhiInterpType PIT, PushDiagToggle PDT>
KOKKOS_INLINE_FUNCTION void	derivs_single_with_e_elec_w_diag (const Grid< Device > &grid, const PushControls &push_controls, const Species< Device > &species, const MagneticField< Device > &magnetic_field, const GridFieldPack< Device, PIT > &gfpack, const PerturbedBField< Device > &perturbed_B_field, SimdParticles &part, SimdPhase &dy, SimdGridWeights< Order::One, PIT > &grid_wts, const FieldAlignedLocalFields< KT, PIT > &E_mag, double time, VFDiag< PDT > &vf_diag)

template<class Device , KinType KT, PhiInterpType PIT>
KOKKOS_INLINE_FUNCTION void	derivs_single_with_e_elec (const Grid< Device > &grid, const PushControls &push_controls, const Species< Device > &species, const MagneticField< Device > &magnetic_field, const GridFieldPack< Device, PIT > &gfpack, const PerturbedBField< Device > &perturbed_B_field, SimdParticles &part, SimdPhase &dy, SimdGridWeights< Order::One, PIT > &grid_wts, const FieldAlignedLocalFields< KT, PIT > &E_mag, double time)

KOKKOS_INLINE_FUNCTION void	get_curl_B (const SimdVector &bfield, const SimdVector(&jacb)[3], double inv_r, double(&curl_B)[3], int i_simd)

KOKKOS_INLINE_FUNCTION double	get_nb_curl_nb (double over_B2, const SimdVector &bfield, const SimdVector(&jacb)[3], double inv_r, int i_simd)

KOKKOS_INLINE_FUNCTION double	get_vspace_jacb_factor (double rho_mod, double over_B2, const SimdVector &bfield, const SimdVector(&jacb)[3], double inv_r, int i_simd)

KOKKOS_INLINE_FUNCTION void	get_curl_nb (const SimdVector &bfield, const double(&grad_B)[3], double over_B, double over_B2, const double(&curl_B)[3], double inv_r, double(&curl_nb)[3], int i_simd)

KOKKOS_INLINE_FUNCTION void	get_grad_B (const SimdVector &bfield, const SimdVector(&jacb)[3], double over_B, double(&grad_B)[3], int i_simd)

KOKKOS_INLINE_FUNCTION void	calculate_Bstar (double drift_on, double rho, const SimdVector &bfield, double B, double over_B, const LocalFields &fld, const double(&curl_nb)[3], double(&Bs)[3], double(&Bsa)[3], int i_simd)

KOKKOS_INLINE_FUNCTION void	get_ExB_drift (double D, double drift_on, const SimdVector &bfield, const SimdVector &E, int i_simd, double over_B2, double inv_r, double(&yp_exb)[3])

template<class Device , PushDiagToggle PDT>
KOKKOS_INLINE_FUNCTION void	derivs_elec (const PushControls &push_controls, const Species< Device > &species, const SimdParticles &part, const Simd< double > &inv_r_vec, const SimdVector &bfield, const Simd< double > &B_mag, const SimdVector(&jacb)[3], const Simd< double > &psi, const SimdVector2D &gradpsi, const SimdVector &tdb, const LocalFields &fld, const LocalEquilProfiles< MT_GLOBAL > &local_eq_profiles, SimdPhase &yprime, VFDiag< PDT > &vf_diag)

Function Documentation

KOKKOS_INLINE_FUNCTION void calculate_Bstar	(	double	drift_on,
		double	rho,
		const SimdVector &	bfield,
		double	B,
		double	over_B,
		const LocalFields &	fld,
		const double(&)	curl_nb[3],
		double(&)	Bs[3],
		double(&)	Bsa[3],
		int	i_simd
	)

Here is the caller graph for this function:

template<class Device , PushDiagToggle PDT>

KOKKOS_INLINE_FUNCTION void derivs_elec	(	const PushControls &	push_controls,
		const Species< Device > &	species,
		const SimdParticles &	part,
		const Simd< double > &	inv_r_vec,
		const SimdVector &	bfield,
		const Simd< double > &	B_mag,
		const SimdVector(&)	jacb[3],
		const Simd< double > &	psi,
		const SimdVector2D &	gradpsi,
		const SimdVector &	tdb,
		const LocalFields &	fld,
		const LocalEquilProfiles< MT_GLOBAL > &	local_eq_profiles,
		SimdPhase &	yprime,
		VFDiag< PDT > &	vf_diag
	)

Get the phase derivatives of a vector of particles

Parameters

[in]	push_controls	Set of options/parameters used in the particle push
[in]	species	Contains species parameters
[in]	part	Vector of particles
[in]	inv_r_vec	Vector of 1/r (1/r_axis in cylindrical limit)
[in]	bfield	Vector of magnetic fields
[in]	B_mag	Vector of bfield magnitudes
[in]	jacb	Vector of jacobians
[in]	psi	Vector of psi coordinates
[in]	gradpsi	Vector of psi gradients (r,z)
[in]	tdb	Vector of perturbed magnetic field
[in]	fld	Contains electromagnetic fields at a set of points
[in]	local_equil_profiles	Contains f0 profile information is required for conventional delta-f method (in dw/dt)
[out]	yprime	Vector of derivatives of the particle phases
[out]	vf_diag	Optional vector of quantities for diagnostic

Here is the call graph for this function:

template<class Device , KinType KT, PhiInterpType PIT>

KOKKOS_INLINE_FUNCTION void derivs_single_with_e_elec	(	const Grid< Device > &	grid,
		const PushControls &	push_controls,
		const Species< Device > &	species,
		const MagneticField< Device > &	magnetic_field,
		const GridFieldPack< Device, PIT > &	gfpack,
		const PerturbedBField< Device > &	perturbed_B_field,
		SimdParticles &	part,
		SimdPhase &	dy,
		SimdGridWeights< Order::One, PIT > &	grid_wts,
		const FieldAlignedLocalFields< KT, PIT > &	E_mag,
		double	time
	)

Here is the call graph for this function:

Here is the caller graph for this function:

template<class Device , KinType KT, PhiInterpType PIT, PushDiagToggle PDT>

KOKKOS_INLINE_FUNCTION void derivs_single_with_e_elec_w_diag	(	const Grid< Device > &	grid,
		const PushControls &	push_controls,
		const Species< Device > &	species,
		const MagneticField< Device > &	magnetic_field,
		const GridFieldPack< Device, PIT > &	gfpack,
		const PerturbedBField< Device > &	perturbed_B_field,
		SimdParticles &	part,
		SimdPhase &	dy,
		SimdGridWeights< Order::One, PIT > &	grid_wts,
		const FieldAlignedLocalFields< KT, PIT > &	E_mag,
		double	time,
		VFDiag< PDT > &	vf_diag
	)

Compute field quantities needed for the push, then output the phase derivatives of a vector of particles

Parameters

[in]	grid	Spatial grid object
[in]	push_controls	Set of options/parameters used in the particle push
[in]	species	Contains species parameters
[in]	magnetic_field	Magnetic field object
[in]	gfpack	GridFieldPack object
[in]	perturbed_B_field
[in,out]	part	Vector of particles
[out]	dy	Vector of derivatives of the particle phases
[in]	grid_wts	Vector of previous triangle index (guess it hasnt moved, to reduce search time)
[in]	E_mag
[in]	time	Global simulation time, not used at the moment
[out]	vf_diag	Optional vector of quantities for diagnostic

Here is the call graph for this function:

Here is the caller graph for this function:

KOKKOS_INLINE_FUNCTION void get_curl_B	(	const SimdVector &	bfield,
		const SimdVector(&)	jacb[3],
		double	inv_r,
		double(&)	curl_B[3],
		int	i_simd
	)

Here is the caller graph for this function:

KOKKOS_INLINE_FUNCTION void get_curl_nb	(	const SimdVector &	bfield,
		const double(&)	grad_B[3],
		double	over_B,
		double	over_B2,
		const double(&)	curl_B[3],
		double	inv_r,
		double(&)	curl_nb[3],
		int	i_simd
	)

Here is the caller graph for this function:

KOKKOS_INLINE_FUNCTION void get_ExB_drift	(	double	D,
		double	drift_on,
		const SimdVector &	bfield,
		const SimdVector &	E,
		int	i_simd,
		double	over_B2,
		double	inv_r,
		double(&)	yp_exb[3]
	)

Get the phase derivatives of a vector of particles

Parameters

[in]	D
[in]	drift_on
[in]	bfield	Vector of magnetic fields
[in]	E
[in]	i_simd	vector index
[in]	over_B2	1/\|B\|^2
[in]	inv_r	1/r
[out]	yp_exb

Here is the caller graph for this function:

KOKKOS_INLINE_FUNCTION void get_grad_B	(	const SimdVector &	bfield,
		const SimdVector(&)	jacb[3],
		double	over_B,
		double(&)	grad_B[3],
		int	i_simd
	)

Here is the caller graph for this function:

KOKKOS_INLINE_FUNCTION double get_nb_curl_nb	(	double	over_B2,
		const SimdVector &	bfield,
		const SimdVector(&)	jacb[3],
		double	inv_r,
		int	i_simd
	)

Get normalized b dot curl of normalized b

Parameters

[in]	over_B2	1/\|B\|^2
[in]	bfield	Vector of magnetic fields
[in]	jacb	Vector of jacobians
[in]	inv_r	1/r
[in]	i_simd	vector index

Returns: double normalized b dot curl of normalized b

Here is the caller graph for this function:

KOKKOS_INLINE_FUNCTION double get_vspace_jacb_factor	(	double	rho_mod,
		double	over_B2,
		const SimdVector &	bfield,
		const SimdVector(&)	jacb[3],
		double	inv_r,
		int	i_simd
	)

Here is the call graph for this function:

Here is the caller graph for this function:

Functions

Function Documentation