push_physics.hpp

Go to the documentation of this file.

#ifndef PUSH_PHYSICS_HPP
#define PUSH_PHYSICS_HPP
#include "push_controls.hpp"
#include "electric_field.hpp"
#include "perturbed_B_field.hpp"
#include "field_aligned_local_fields.hpp"
#include "push_diagnostic.hpp"
#include "local_equil_profiles.hpp"

// Overload derivs_single_with_e_elec with and without diagnostics
template<class Device, KinType KT, PhiInterpType PIT, MarkerType MT, MagneticFieldMode MFM, BFieldSymmetry BFS, typename GFPT, 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, GFPT> &gfpack,
        const PerturbedBField<Device> &perturbed_B_field, SimdParticles &part,SimdPhase &dy, 
        SimdGridWeights<Order::One, PIT>& grid_wts, const FieldAlignedLocalFields<KT, PIT, GFPT>& E_mag, double time, VFDiag<PDT>& vf_diag);

template<class Device, KinType KT, PhiInterpType PIT, MarkerType MT, MagneticFieldMode MFM, BFieldSymmetry BFS, typename GFPT>
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, GFPT> &gfpack,
        const PerturbedBField<Device> &perturbed_B_field, SimdParticles &part,SimdPhase &dy, 
        SimdGridWeights<Order::One, PIT>& grid_wts, const FieldAlignedLocalFields<KT, PIT, GFPT>& E_mag, double time);


template<class Device, MarkerType MT, MagneticFieldMode MFM, BFieldSymmetry BFS, typename GFPT, 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<GFPT>& fld, const LocalEquilProfiles<MT>& local_eq_profiles,
        SimdPhase &yprime, VFDiag<PDT>& vf_diag);

#include "push_physics.tpp"
#include "em_push_physics.tpp"

#endif