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, 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);


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);

#include "push_physics.tpp"

#if defined(EXPLICIT_EM) || defined(DELTAF_CONV)
#include "em_push_physics.tpp"
#endif

#endif