pseudo_inverse_diag.hpp

Go to the documentation of this file.

 
#ifndef PSEUDO_INVERSE_DIAG_HPP
#define PSEUDO_INVERSE_DIAG_HPP
 
#include "sml.hpp"
#include "grid.hpp"
#include "magnetic_field.hpp"
#include "species.hpp"
#include "sources.hpp"
#include "diagnostic.hpp"
 
 
class PseudoInvDiagnostics : public Diagnostic {
 
    public:
 
    PseudoInvDiagnostics() : 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_pseudo_inv_on", false); 
 
        if (is_on) Diagnostic::init(name, f_source_period);
    }
 
    void write(const Simulation<DeviceType> &sml, const Grid<DeviceType>& grid,
                         const MagneticField<DeviceType>& magnetic_field,
                         const Species<DeviceType>& species,
                         const DomainDecomposition<DeviceType>& pol_decomp,
                         uint64_t plane, int n_species,
                         int n_nodes_all, int n_nodes_pseudo_inv_converged_all,
                         int n_nodes_pseudo_inv_diverged_all, int n_nodes_pseudo_inv_skipped_all,
                         View<int*, CLayout, HostType> pseudo_inv_num_its,
                         View<bool*, CLayout, HostType> pseudo_inv_converged_node
                         );
};
 
#endif