increase_ptl_tor.hpp

Go to the documentation of this file.

#ifndef INCREASE_PTL_TOR_HPP
#define INCREASE_PTL_TOR_HPP

#include "grid.hpp"
#include "plasma.hpp"

struct IncreasePtlTorSettings{
    int increase_factor; 
    double rand_phi;
    double rand_w0;
    double rand_w1;
    double rand_vp;
    double rand_mu;

    IncreasePtlTorSettings(){}

    IncreasePtlTorSettings(NLReader::NamelistReader& nlr){
        nlr.use_namelist("sml_param");
        increase_factor = nlr.get<int>("sml_increase_ptl_factor_int", 1); 
        rand_w0 = nlr.get<double>("sml_increase_ptl_rand_w0", 1.0e-4); 
        rand_w1 = nlr.get<double>("sml_increase_ptl_rand_w1", 1.0e-4); 
        rand_mu = nlr.get<double>("sml_increase_ptl_rand_mu", 1.0e-4); 
        rand_vp = nlr.get<double>("sml_increase_ptl_rand_vp", 1.0e-4); 
        rand_phi = nlr.get<double>("sml_increase_ptl_rand_phi", 1.0e-1); 
    }
};

void increase_ptl_tor(const MagneticField<DeviceType>& magnetic_field, const Grid<DeviceType>& grid, Plasma& plasma, const IncreasePtlTorSettings& settings);

#endif