ndg_moment_slope_limiter_base.h

Go to the documentation of this file.

#ifndef WMNDG_MOMENT_SLOPE_LIMITER_BASE_H
#define WMNDG_MOMENT_SLOPE_LIMITER_BASE_H
 
// Wm includes
#include "dfem/tools/tools.h"
#include "dfem/variable_adjusters/variable_adjuster.h"
 
#include <memory>
 
namespace wxm
{
namespace dfem
{
namespace variable_adjuster
{
namespace limiter
{
 
class ndg_moment_slope_limiter_base
    : public wxm::dfem::variable_adjuster::variable_adjuster_t
{
public:
    /* Constructor */
    ndg_moment_slope_limiter_base() = default;
 
    virtual void setup(const WxCryptSet& wxc) override;
 
    void solve(real time, variables_type& input) override;
 
protected:
    // void primitive_kernel(const WxRange& element_scope_range, real time,
    //                    const std::vector<wxm::array::patch_array_t*>& input);
 
    void primitive_kernel(const WxRange& element_scope_range,
                          real time,
                          std::string direction,
                          const std::vector<std::vector<real>>& input,
                          std::vector<std::vector<real>>& output);
 
    virtual void averaging_kernel(const WxRange& element_scope_range,
                                  const WxRange& face_scope_range,
                                  const real time,
                                  const std::vector<wxm::array::patch_array_t*>& input);
 
    virtual void limiting_kernel(const WxRange& element_scope_range,
                                 real time,
                                 std::vector<wxm::array::patch_array_t*>& input) = 0;
 
    virtual void allocate();
    virtual void zero();
 
    // variables associated with the solve
    std::vector<real>
        _global_averages; // num_elems x num_vars (averages for conserved vars q)
    std::vector<real>
        _global_centers; // num_elems x num_vars (centers for conserved vars q)
    std::vector<real>
        _global_maximums; // result of equation 6.3 num_elems x num_vars (for conserved q)
    std::vector<real>
        _global_minimums; // result of equation 6.3 num_elems x num_vars (for conserved q)
 
    std::vector<real>
        _global_w_averages; // num_elems x num_vars (average for primitive vars w)
    std::vector<real>
        _global_w_centers; // num_elems x num_vars (centers for primitive vars w)
    std::vector<real> _global_w_maximums; // result of equation 6.3 num_elems x num_vars
                                          // (for primitives w)
    std::vector<real> _global_w_minimums; // result of equation 6.3 num_elems x num_vars
                                          // (for primitives w)
 
    // Applications
    wxm::dfem::tools::scope_t _limiting_scope;
    wxm::dfem::tools::scope_t _averaging_scope;
    wxm::dfem::tools::scope_t _averaging_face_scope;
 
    // Spatial basis set name
    std::string _basisSetName;
 
    // Spatial basis set
    WmBasisArraySet _basis_set;
 
    // WmUDGGeometry
    std::unique_ptr<WmUDGGeometry> _dg_geometry;
 
    // primitive variables
    std::vector<std::vector<real>> _w;
    std::vector<std::vector<real>> _in_vec;
    std::vector<std::vector<real>> _out_vec;
 
private:
    ndg_moment_slope_limiter_base&
    operator=(const ndg_moment_slope_limiter_base& var) = delete;
 
    ndg_moment_slope_limiter_base(const ndg_moment_slope_limiter_base& var) = delete;
};
 
} // namespace limiter
} // namespace variable_adjuster
} // namespace dfem
} // namespace wxm
 
#endif