globals.hpp

Go to the documentation of this file.

#ifndef GLOBALS_HPP
#define GLOBALS_HPP
#include "space_settings.hpp"
#include "array_deep_copy.hpp"
#include "access_add.hpp"
#include "simd.hpp"

enum SpeciesType {
    ELECTRON = 0,
    MAIN_ION = 1
};

enum KinType{
    DriftKin=0,
    GyroKin=1
};

enum ParticlePhase {
    PIR = 0, 
    PIZ,     
    PIP,     
    PIRHO,   
    PIW1,    
    PIW2,    
    PTL_NPHASE
};

enum ParticleConsts {
    PIM = 0, 
    PIW0,    
    PIF0,    
    PTL_NCONST
};

const double UNIT_CHARGE = 1.6022e-19; 
const double EV_2_J=UNIT_CHARGE;       
const double J_2_EV=1.0/EV_2_J;        
const double PROTON_MASS = 1.6720e-27;            //< proton mass (MKS)
const double PI        = 3.1415926535897932;    //< pi
const double TWOPI     = 6.2831853071795862;    //< 2 pi

// Diagnostics
const int DIAG_HEAT_NVAR = 5;
const int N_VF_DIAG = 9;
#ifdef XGC1
  const int D1D_NPV1 = 16;
#else
  const int D1D_NPV1 = 18;
#endif

/*** SIMD structures (probably temporary) ****/


struct SimdVector2D{
    Simd<double> r;
    Simd<double> z;
};

struct SimdVector{
    Simd<double> r;
    Simd<double> z;
    Simd<double> phi;

    // Vector operation, sets SimdVector to 0.0
    KOKKOS_INLINE_FUNCTION void zero(){
        for (int i_simd = 0; i_simd<SIMD_SIZE; i_simd++){
            r[i_simd] = 0.0;
            z[i_simd] = 0.0;
            phi[i_simd] = 0.0;
        }
    }
};

/* Fortran pointers */
// This type is used to refer to fortran pointers that we never plan on using (i.e. dereferencing) inside the C++ code.
// The reason we don't use them in the C++ code is that the objects being referred to (spall, grid, psn), are very complicated.
// In order to use them, we would need to carefully create the equivalent object in C++, and since developers tend to add members
// to those Fortran objects on the fly, using them in C++ would be prone to bugs.
// Here we tell C++ that the pointer is to an integer, i.e. int *. This could be anything, but int works.
typedef int* FortranPtr;

#endif