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