loop_voltage_diag.hpp

Go to the documentation of this file.

#ifndef LOOP_VOLTAGE_DIAG_HPP
#define LOOP_VOLTAGE_DIAG_HPP

#include "sml.hpp"
#include "grid.hpp"
#include "electric_field.hpp"
#include "sources.hpp"
#include "diagnostic.hpp"


class LoopVolDiagnostics : public Diagnostic {

    public:

    LoopVolDiagnostics() : Diagnostic() {}

    void init(NLReader::NamelistReader& nlr, int f_source_period, const std::string& name){
        nlr.use_namelist("diag_param");
        bool is_on = nlr.get<bool>("diag_current_drive_on",false); 
        int loop_vol_diag_period = f_source_period*nlr.get<int>("diag_current_drive_period",1); 

        if (is_on) Diagnostic::init(name, loop_vol_diag_period);
    }

    void write(const Simulation<DeviceType> &sml, const Grid<DeviceType>& grid,
                     const ElectricField<DeviceType>& electric_field,
                     View<double*, CLayout, HostType>& jpara_diff_current
                     );
};

#endif