1 #ifndef VELOCITY_GRID_HPP
2 #define VELOCITY_GRID_HPP
29 nvp = nlr.
get<
int>(
"f0_nvp", 15);
30 nmu = nlr.
get<
int>(
"f0_nmu", 31);
31 vp_max = nlr.
get<
double>(
"f0_vp_max", 3.0);
45 exit_XGC(
"\nError: Cannot run with f0_velocity_interp_use_pseudo_inv = true without using PETSc.\n");
48 exit_XGC(
"\nError: f0_velocity_interp_pseudo_inv_order must be 0 or larger.\n");
51 std::string msg =
"\nError: Cannot use f0_velocity_interp_pseudo_inv_order larger than ";
53 msg +=
", if you really want to use that large order then modify LAGRANGE_MAX_ORDER in the code.\n";
58 exit_XGC(
"\nError: f0_nmu must be divisible by f0_velocity_interp_pseudo_inv_order when f0_velocity_interp_use_pseudo_inv = true.\n");
61 exit_XGC(
"\nError: 2*f0_nvp must be divisible by f0_velocity_interp_pseudo_inv_order when f0_velocity_interp_use_pseudo_inv = true.\n");
70 KOKKOS_INLINE_FUNCTION
double vpar_n(
int ivz)
const{
75 KOKKOS_INLINE_FUNCTION
double vperp_n(
int ivr)
const{
81 return (ivr==0 || ivr==
nvr-1) ? 0.5 : 1.0;
94 if(vr_wt.
i>=
nmu || vz_wt.
i>=2*
nvp || vz_wt.
i<0){
107 vgrid_wts.
w_00 *= inv_smu0;
108 vgrid_wts.
w_01 *= inv_smu0;
109 vgrid_wts.
w_10 *= inv_smu1;
110 vgrid_wts.
w_11 *= inv_smu1;
115 double vol0 = vol*(vgrid_wts.
i_vr==0 ? 0.5 : 1.0);
116 double vol1 = vol*(vgrid_wts.
i_vr==
nmu-1 ? 0.5 : 1.0);
118 vgrid_wts.
w_00 *= vol0;
119 vgrid_wts.
w_01 *= vol0;
120 vgrid_wts.
w_10 *= vol1;
121 vgrid_wts.
w_11 *= vol1;
126 double inv_vol0 = inv_vol*(vgrid_wts.
i_vr==0 ? 2.0 : 1.0);
127 double inv_vol1 = inv_vol*(vgrid_wts.
i_vr==
nmu-1 ? 2.0 : 1.0);
129 vgrid_wts.
w_00 *= inv_vol0;
130 vgrid_wts.
w_01 *= inv_vol0;
131 vgrid_wts.
w_10 *= inv_vol1;
132 vgrid_wts.
w_11 *= inv_vol1;
int element_order
velocity grid finite element order (0 = nearest neighbor), (1 = linear), (2 = quadratic), (3 = cubic), ....
Definition: velocity_grid.hpp:22
double smu_max
max mu
Definition: velocity_grid.hpp:14
int nmu
n points in mu (not including zero)
Definition: velocity_grid.hpp:13
int nvp
n points in parallel velocity (not including zero)
Definition: velocity_grid.hpp:9
T get(const string ¶m, const T default_val, int val_ind=0)
Definition: NamelistReader.hpp:373
bool pseudo_inv_on
whether pseudo-inverse interpolation is used in velocity space
Definition: velocity_grid.hpp:24
double w_00
Definition: vgrid_weights.hpp:10
Definition: linear_weights.hpp:8
Definition: velocity_grid.hpp:8
KOKKOS_INLINE_FUNCTION void weight_by_inv_vol(double inv_vol, VGridWeights &vgrid_wts) const
Definition: velocity_grid.hpp:125
Definition: NamelistReader.hpp:193
int nvr
full grid size (including zero)
Definition: velocity_grid.hpp:19
double w_01
Definition: vgrid_weights.hpp:11
KOKKOS_INLINE_FUNCTION double vpar_n(int ivz) const
Definition: velocity_grid.hpp:70
double w_11
Definition: vgrid_weights.hpp:13
KOKKOS_INLINE_FUNCTION double mu_vol_fac(int ivr) const
Definition: velocity_grid.hpp:80
void use_namelist(const string &namelist)
Definition: NamelistReader.hpp:355
#define LAGRANGE_MAX_ORDER
Definition: lagrange_weights.hpp:18
int i
Definition: linear_weights.hpp:9
double vp_max
max parallel velocity
Definition: velocity_grid.hpp:10
int i_vr
Definition: vgrid_weights.hpp:8
KOKKOS_INLINE_FUNCTION void weight_by_vol(double vol, VGridWeights &vgrid_wts) const
Definition: velocity_grid.hpp:114
void exit_XGC(std::string msg)
Definition: globals.hpp:37
VelocityGrid(NLReader::NamelistReader &nlr)
Definition: velocity_grid.hpp:26
Definition: vgrid_weights.hpp:7
double dsmu
grid spacing in mu
Definition: velocity_grid.hpp:15
double inv_mu0_factor
Set value of lowest mu in grid –> 1/mu0_factor.
Definition: velocity_grid.hpp:17
int nvz
full grid size (including negative and zero)
Definition: velocity_grid.hpp:20
KOKKOS_INLINE_FUNCTION VGridWeights get_weights(double smu, double vp) const
Definition: velocity_grid.hpp:85
double dvp
grid spacing in parallel velocity
Definition: velocity_grid.hpp:11
double w_10
Definition: vgrid_weights.hpp:12
KOKKOS_INLINE_FUNCTION void weight_by_inv_dsmu(VGridWeights &vgrid_wts) const
Definition: velocity_grid.hpp:103
KOKKOS_INLINE_FUNCTION double vperp_n(int ivr) const
Definition: velocity_grid.hpp:75