ndg_moe_rossmanith.h

Go to the documentation of this file.

#ifndef WMNDG_MOE_ROSSMANITH_H
#define WMNDG_MOE_ROSSMANITH_H
 
// Wm includes
#include "dfem/tools/tools.h"
#include "dfem/variable_adjusters/variable_adjuster.h"
#include "dfem/variable_adjusters/limiter/ndg_moment_slope_limiter_base.h"
 
#include <memory>
 
namespace wxm
{
namespace dfem
{
namespace variable_adjuster
{
namespace limiter
{
 
class ndg_moe_rossmanith
    : public wxm::dfem::variable_adjuster::limiter::ndg_moment_slope_limiter_base
{
public:
    /* Constructor */
    ndg_moe_rossmanith() = default;
 
    virtual void setup(const WxCryptSet& wxc) override;
 
protected:
    real
    y_func(const real extremum, const real qi_bar, const real element_extremum) const;
    real phi_func(const real y) const;
    real alpha_func(const int elementIndex) const;
 
    virtual void limiting_kernel(const WxRange& element_scope_range,
                                 real time,
                                 std::vector<wxm::array::patch_array_t*>& input) override;
 
    virtual void allocate() override;
 
    std::vector<real>
        _neighbor_upper_bound; // result of equation 6.4 num_elems x num_vars
    std::vector<real>
        _neighbor_lower_bound;    // result of equation 6.5 num_elems x num_vars
    std::vector<real> _theta_M_i; // result of equation 6.6 num_elems
    std::vector<real> _theta_m_i; // result of equation 6.6 num_elems
    std::vector<real> _theta_i;   // result of equation 6.7 num_elems
 
    // relevent parameter
    real _alpha_coefficient;
 
private:
    ndg_moe_rossmanith& operator=(const ndg_moe_rossmanith& var) = delete;
 
    ndg_moe_rossmanith(const ndg_moe_rossmanith& var) = delete;
};
 
} // namespace limiter
} // namespace variable_adjuster
} // namespace dfem
} // namespace wxm
 
#endif