charge.hpp

Go to the documentation of this file.

#ifndef CHARGE_HPP
#define CHARGE_HPP
#include "space_settings.hpp"
#include "electric_field.hpp"
#include "perturbed_B_field.hpp"

// Charge class
template<class Device>
class Charge {
    public:

    // Constructor
    Charge(int nnode,int npsi00,int nrho,int nlarmor, int ntor, int nthreads);

    // Default constructor
    Charge(){}

    // Arrays for particle-to-grid operations
#ifdef XGC1
    Kokkos::View<double***,Kokkos::LayoutRight,Device> edensity;   
#ifdef EXPLICIT_EM
    Kokkos::View<double***,Kokkos::LayoutRight,Device> ejpar;   
    Kokkos::View<double****,Kokkos::LayoutRight,Device> ijpar;   
#endif
    Kokkos::View<double****,Kokkos::LayoutRight,Device> idensity;   
#else
    Kokkos::View<double**,Kokkos::LayoutRight,Device> edensity;   
    Kokkos::View<double***,Kokkos::LayoutRight,Device> ejpar_re;
    Kokkos::View<double***,Kokkos::LayoutRight,Device> ejpar_im;

    Kokkos::View<double***,Kokkos::LayoutRight,Device> idensity;   
    Kokkos::View<double***,Kokkos::LayoutRight,Device> ijpar_re;
    Kokkos::View<double***,Kokkos::LayoutRight,Device> ijpar_im;
#endif

    void reset_charge() const;
    void reset_ion_charge() const;
};

#include "charge.tpp"

#endif