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"

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, Simd<double> *vf_diag=nullptr );

#if defined(EXPLICIT_EM) || defined(DELTAF_CONV)
#include "em_push_physics.hpp"
#else
template<class Device>
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 SimdVector &efield,
#ifdef SONIC_GK
    const SimdVector &du2_E,
#endif
        SimdPhase &yprime, Simd<double> *vf_diag  );
#endif

#include "push_physics.tpp"
#endif