solver_poisson_axisym.hpp

Go to the documentation of this file.

#ifndef SOLVER_POISSON_AXISYM_HPP
#define SOLVER_POISSON_AXISYM_HPP
 
#include "solver.hpp"
 
class SolverPoissonAxisym : public Solver{
 
    public:
 
    bool get_is_iter_solver(NLReader::NamelistReader& nlr) const{
        nlr.use_namelist("sml_param");
        bool iter_solver = nlr.get<bool>("sml_iter_solver", true); 
                                                  // non-axisymmetric parts. The axisymmetric part has the
                                                  // flux-surface averaged potential on the RHS. Solve
                                                  // iteratively by lagging the RHS potential by one
                                                  // iteration. Stop after sml_iter_solver_niter iterations.
                                                  // If .false., split Poisson equation into
                                                  // flux-surface averaged part and rest (assuming that the
                                                  // polarization term commutes with the flux-surface
                                                  // average) and the rest; solve in two steps.
        return iter_solver;
    }
 
    SolverPoissonAxisym(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, P0, -1, spectral_in, get_is_iter_solver(nlr) ? 2 : 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.cbd0");
            write_solver_boundary(bd, "Efield_excluded_bnd.pbd0");
        }
    }
 
    void setup_boundaries(NLReader::NamelistReader& nlr, bool is_XGCa, bool explicit_electromagnetic, const MagneticField<DeviceType>& magnetic_field, const Grid<DeviceType>& grid);
};
 
#endif