ndg_tu_alibadi.h

Go to the documentation of this file.

#ifndef WMNDG_TU_ALIBADI_H
#define WMNDG_TU_ALIBADI_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_tu_alibadi
    : public wxm::dfem::variable_adjuster::limiter::ndg_moment_slope_limiter_base
{
public:
    /* Constructor */
    ndg_tu_alibadi() = default;
 
    virtual void setup(const WxCryptSet& wxc) override;
 
    // /**
    //  * @brief Solves the spatial system and puts result in _rhs
    //  */
    // void solve(real time, variables_type& input) override;
 
protected:
    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) override;
 
    virtual void limiting_kernel(const WxRange& element_scope_range,
                                 real time,
                                 std::vector<wxm::array::patch_array_t*>& input) override;
 
    void face_average(const WxRange& face_scope_range,
                      const real time,
                      const std::vector<wxm::array::patch_array_t*>& input);
 
    void face_gradient(const WxRange& face_scope_range,
                       const std::vector<wxm::array::patch_array_t*>& input);
 
    void unlimited_centroid_gradient(const WxRange& element_scope_range);
    void limited_centroid_gradient(const WxRange& element_scope_range);
    void reconstruct_primitives(const WxRange& element_scope_range);
 
    std::vector<std::vector<real>>
        _q_avg_face_nodes; // [num_vars] x [num_faces * num_face_nodes]
    std::vector<std::vector<real>>
        _w_avg_face_nodes; // [num_vars] x [num_faces * num_face_nodes]
 
    std::vector<std::vector<real>> _w_face_gradients; // [num_vars] x [num_faces*3]
    std::vector<std::vector<real>>
        _w_unlimited_centroid_gradients; // [num_vars] x [num_elements*3]
    std::vector<std::vector<real>>
        _w_limited_centroid_gradients; // [num_vars] x [num_elements*3]
 
    std::vector<std::vector<real>> _w_limited;  // [num_vars] x [num_elements * N_p]
    std::vector<std::vector<real>> _dw_limited; // [num_vars] x [num_elements * N_p]
 
    virtual void allocate() override;
    virtual void zero() override;
 
    // static std::mutex _mutex;
 
private:
    ndg_tu_alibadi& operator=(const ndg_tu_alibadi& var) = delete;
 
    ndg_tu_alibadi(const ndg_tu_alibadi& var) = delete;
};
// std::mutex ndg_tu_alibadi::_mutex;
 
} // namespace limiter
} // namespace variable_adjuster
} // namespace dfem
} // namespace wxm
 
#endif