solver_ampere_cv.hpp

Go to the documentation of this file.

#ifndef SOLVER_AMPERE_CV_HPP
#define SOLVER_AMPERE_CV_HPP
 
#include "solver.hpp"
 
class SolverAmpereCV : public Solver{
 
    double skin_norm;
 
    public:
 
    SolverAmpereCV(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_CV, A_CV_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) );
 
 
        // 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");
        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