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"
#include "diag_f0_df.hpp"

extern "C" void initialize_diffusion(int* bd);

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

    public:

    Diffusion(){}

    Diffusion(NLReader::NamelistReader& nlr, const MagneticField<DeviceType>& magnetic_field, const Grid<DeviceType>& grid);

    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,
              DiagF0DF& diag_f0_df, const Moments& moments,
              const VGridDistribution<DeviceType>& f0_f, VGridDistribution<HostType>& f0_df0g,
              double dt);
};

#endif