neutrals.hpp

Go to the documentation of this file.

#ifndef NEUTRALS_HPP
#define NEUTRALS_HPP

#include "view_arithmetic.hpp"
#include "xgc_io.hpp"

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

    // Neutrals variables
    bool use_ion_loss;         
    View<double*,CLayout,Device> weight_sum_lost;   
    View<double**,CLayout,Device> weight_wall_lost;  

    View<double*,CLayout,HostType> weight_sum_lost_h;   
    View<double*,CLayout,HostType> weight_wall_lost_h;  
    View<double*,CLayout,HostType> node_dens_h;
    View<double*,CLayout,HostType> node_temp_h;

    // Constructor
    Neutrals( int nthreads, int nbins, int nnode);

    // "Null constructor"
    Neutrals(){}

    void write_checkpoint_files(const XGC_IO_Stream& stream) const;

    void read_checkpoint_files(const XGC_IO_Stream& stream);

    void copy_to_host() const;

    void clear_device_views() const;
};
#endif