solver_poisson_turb.hpp

Go to the documentation of this file.

#ifndef SOLVER_POISSON_TURB_HPP
#define SOLVER_POISSON_TURB_HPP
 
#include "solver.hpp"
 
class SolverPoissonTurb : public Solver{
 
    public:
 
    SolverPoissonTurb(NLReader::NamelistReader& nlr, bool is_XGCa, bool is_em, const MagneticField<DeviceType>& magnetic_field, const Grid<DeviceType>& grid, std::string label_in, bool spectral_in, bool use_pade_in, bool update_solver, int update_solver_nstep)
      : Solver(label_in, H, H_SPEC, spectral_in, 1, grid.ntriangle, use_pade_in, update_solver, update_solver_nstep)
    {
        TIMER("EXTEND_BOUNDARY",
            setup_boundaries(nlr, is_XGCa, is_em, magnetic_field, grid) );
 
        if(is_rank_zero()){
            write_solver_boundary(bd_rhs, "Efield_excluded_bnd.cbdH");
            write_solver_boundary(bd, "Efield_excluded_bnd.pbdH");
        }
    }
 
    void lhs_set_alpha_beta(const SolverInitData<HostType>& solver_data, const View<double*,CLayout,HostType>& alpha, const View<double*,CLayout,HostType>& beta, double& scale) override;
 
    void setup_boundaries(NLReader::NamelistReader& nlr, bool is_XGCa, bool explicit_electromagnetic, const MagneticField<DeviceType>& magnetic_field, const Grid<DeviceType>& grid);
};
 
#endif