velocity_grid.hpp

Go to the documentation of this file.

#ifndef VELOCITY_GRID_HPP
#define VELOCITY_GRID_HPP

// Specifies the dimensions of the velocity grid
struct VelocityGrid{
    int nvp;        //< n points in parallel velocity (not including zero)
    double vp_max;  //< max parallel velocity
    double dvp;     //< grid spacing in parallel velocity

    int nmu;        //< n points in mu (not including zero)
    double smu_max; //< max mu
    double dsmu;    //< grid spacing in mu

    double inv_mu0_factor = 1.0/3.0; //< Set value of lowest mu in grid --> 1/mu0_factor

    int nvr;        //< full grid size (including zero)
    int nvz;        //< full grid size (including negative and zero)

    VelocityGrid(){}

    VelocityGrid(int f0_nvp, int f0_nmu, double f0_vp_max, double f0_smu_max) 
            : nvp(f0_nvp),
              nmu(f0_nmu),
              vp_max(f0_vp_max),
              smu_max(f0_smu_max),
              dvp(vp_max/nvp),
              dsmu(smu_max/nmu),
              nvr(nmu+1),
              nvz(nvp*2+1){}

};
#endif