solver_ampere.hpp

Go to the documentation of this file.

#ifndef SOLVER_AMPERE_HPP
#define SOLVER_AMPERE_HPP
 
#include "solver.hpp"
 
class SolverAmpere : public Solver{
 
    double skin_norm;
 
    public:
 
    bool em_mixed_variable;
 
    SolverAmpere(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, A, A_SPEC, spectral_in, 2, 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.jbdH");
            write_solver_boundary(bd, "Efield_excluded_bnd.AbdH");
        }
 
        // skin term correction (finite sphere)
        // note sml_flat_marker_cutoff1 (para) can in principle be different
        // from sml_flat_marker_cutoff2 (perp)
        nlr.use_namelist("sml_param");
        em_mixed_variable = nlr.get<bool>("sml_em_mixed_variable",true);
        double v_para_cutoff = nlr.get<double>("sml_flat_marker_cutoff1", 4.0);
        skin_norm = std::erf(v_para_cutoff) - (1.0 + v_para_cutoff*v_para_cutoff*2.0/3.0)*2.0/sqrt(PI)*v_para_cutoff*exp(-v_para_cutoff*v_para_cutoff);
    }
 
    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 rhs2_set_alpha_beta(const SolverInitData<HostType>& solver_data, double scale, const View<double*,CLayout,HostType>& alpha, const View<double*,CLayout,HostType>& beta) override;
 
    void setup_boundaries(NLReader::NamelistReader& nlr, bool is_XGCa, bool explicit_electromagnetic, const MagneticField<DeviceType>& magnetic_field, const Grid<DeviceType>& grid);
};
 
#endif