volumes_and_areas.hpp

Go to the documentation of this file.

#ifndef VOLUMES_AND_AREAS_HPP
#define VOLUMES_AND_AREAS_HPP

#include "globals.hpp"
#include "domain_decomposition.hpp"

// Volume- and Area-related Host Views. Should find a different location for these so that these unused views aren't sent to GPU
struct VolumesAndAreas{
    View<double*,CLayout,HostType> tr_area_h;
    View<double*,CLayout,HostType> tr_vol_h;
    View<double**,CLayout,HostType> node_vol_ff_h; // Default to midplane, for now
    View<double**,CLayout,HostType> node_vol_ff_l_h; // Field following volume for left plane nodes
    View<double**,CLayout,HostType> node_vol_ff_r_h; // Field following volume for right plane nnodes
    View<double*,CLayout,HostType> node_vol_nearest_h;
    View<double*,CLayout,HostType> node_vol_h;
    View<double*,CLayout,HostType> node_area_h;
    View<double*,CLayout,HostType> vol00_h;

    VolumesAndAreas(){}

    VolumesAndAreas(int ntriangle_in, int nnode_in, int npsi00_in);

    void set_fortran_ptrs() const;

    bool saved_file_exists(const std::string& input_dir) const;

    void read(const DomainDecomposition<DeviceType>& pol_decomp, bool is_axisymmetric, const std::string& input_dir);

    void write(const DomainDecomposition<DeviceType>& pol_decomp, bool is_axisymmetric) const;
};

#endif