diffusion.hpp

Go to the documentation of this file.

#ifndef DIFFUSION_HPP
#define DIFFUSION_HPP

#include "sml.hpp"
#include "grid.hpp"
#include "magnetic_field.hpp"
#include "electric_field.hpp"
#include "domain_decomposition.hpp"
#include "velocity_grid.hpp"
#include "plasma.hpp"
#include "vgrid_distribution.hpp"
#include "source.hpp"

class Diffusion : public Source{
    double inner_psi_bound;
    double outer_psi_bound;

    public:

    Diffusion(){}

    Diffusion(NLReader::NamelistReader& nlr, const MagneticField<DeviceType>& magnetic_field){
        nlr.use_namelist("diff_param");
        inner_psi_bound = nlr.get<double>("diff_bd_in",magnetic_field.inpsi);
        outer_psi_bound = nlr.get<double>("diff_bd_out",magnetic_field.outpsi);

        // Normalize by xpt
        // Default values are already normalized, so only normalize if not default
        if(nlr.present("diff_bd_in")) inner_psi_bound *= magnetic_field.equil.xpt_psi;
        if(nlr.present("diff_bd_out")) outer_psi_bound *= magnetic_field.equil.xpt_psi;
    }


    void apply(const Simulation<DeviceType>& sml, const Grid<DeviceType>& grid, const MagneticField<DeviceType>& magnetic_field,
              ElectricField<DeviceType>& electric_field,
              const DomainDecomposition<DeviceType>& pol_decomp, const VelocityGrid& vgrid, Plasma& plasma,
              const Moments& moments,
              VGridDistribution<HostType>& f0_f, VGridDistribution<HostType>& f0_df0g,
              VGridDistribution<HostType>& df0g_tmp,
              double dt);
};

#endif