sml.hpp

Go to the documentation of this file.

#ifndef SML_HPP
#define SML_HPP
#include "space_settings.hpp"
#include "NamelistReader.hpp"

// Simulation class
template<class Device>
class Simulation {
    public:

    // Constructor
    Simulation(NLReader::NamelistReader& nlr);

#ifdef XGCA
    static constexpr bool is_XGCa = true; 
#else
    static constexpr bool is_XGCa = false; 
#endif

#ifdef DELTAF_CONV
    static constexpr bool reduced_deltaf = true; 
#else
    static constexpr bool reduced_deltaf = false; 
#endif

#if defined(DELTAF_CONV_SEPARATE_N0) || !defined(DELTAF_CONV)
    static constexpr bool separate_n0 = true; 
#else
    static constexpr bool separate_n0 = false; 
#endif

#ifdef EXPLICIT_EM
    static constexpr bool explicit_electromagnetic = true; 
#else
    static constexpr bool explicit_electromagnetic = false; 
#endif

#ifdef F0_UPDATE_ANALYTIC_OLD
    static constexpr bool old_f0_update_analytic = true; 
#else
    static constexpr bool old_f0_update_analytic = false; 
#endif

    // Input parameters
    int bounce;                
    int bounce_zero_weight;    
    double bounce_buffer;      
    double dt;                 
    int sheath_mode;           
    bool sheath_adjust;        
    bool ptb_3db_on;           
    bool drift_on;             
    bool neutrals;             
    bool diag_heat_on;         
    bool deltaf;               
    bool xz_00_up;             
    int mstep;                 
    int nrk;                   
    bool electron_on;          
    bool f0_grid;              
    bool limit_marker_den;     
    int special;               
    int nthreads;              

    int restart_write_period;  

    bool field_solver_on=true;  

    bool split_weight_scheme;   
    bool dwdt_fix_bg;           // 1-w factor in dwdt
    bool dwdt_exb_only;         // excluding grad-b drift in dwdt

    bool em_mixed_variable;     
    bool em_control_variate;    
    int  em_pullback_mode;      
    bool em_control_variate_final_cv; // Should be local variable/argument
    bool em_control_variate_flag; // Should be local variable/argument
    int em_control_variate_niter; 
    bool em_b_para_eff;               
    bool em_es_step;            
    bool exclude_private;      

    bool grad_psitheta;        
    // Used in getf0
    bool update_poisson_solver; 
    double dpot_te_limit;      

    // Used in update_f0
    bool no_fp_in_f;            
    double f0_grid_alpha;       
    bool f0_update_analytic;     
    bool f0_update_analytic_local;        
    double f0_update_analytic_alpha;      
    double f0_update_analytic_damp_width; 
    int f0_grid_alpha_start;    
    bool ignore_f0g; 
    // Diffusion
    bool diff_on;

    // Resampling
    bool resamp_on; 
    bool resamp_for_final_restart_write; 
    bool resamp_restart_read; 

    // marker setup
    bool field_aligned_initial;
    bool flat_marker;
    double marker_temp_factor;
    double marker_temp_factor2;
    double marker_temp_factor3;
    double flat_marker_decay_start1;
    double flat_marker_decay_start2;
    double flat_marker_cutoff1;
    double flat_marker_cutoff2;
    double flat_marker_width1;
    double flat_marker_width2;
    double initial_deltaf_noise;
    double marker_min_temp;
    double low_mu_fill_population;

    // Simulation variables
    double time;               
    int gstep;                 

    bool multirate_timestepping; 
    int mr_factor[5]; 
    double mr_psi_max[5]; 

    // RMP penetration calculation
    int mstep_rmp_es;           
    int mstep_rmp_em;           
    int rmp_es_to_em_dt_ratio;  

    // Core-edge coupling
    bool cce_coupling_on;       

    // Push performance - could be elsewhere
    bool use_unfused_electron_push_kernel;

    // Loop-voltage (for driving the equilibrium current)
    double loop_voltage;        
    double loop_voltage_psimax; 
};

#include "sml.tpp"
#endif