gradient_interior_penalty_internal_approximation.h

Go to the documentation of this file.

#ifndef WMGRADIENT_INTERIOR_PENALTY_INTERNAL_APPROXIMATION_H
#define WMGRADIENT_INTERIOR_PENALTY_INTERNAL_APPROXIMATION_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 gradient
{
 
class gradient_interior_penalty_internal_approximation
    : public wxm::dfem::variable_adjuster::variable_adjuster_t
{
public:
    gradient_interior_penalty_internal_approximation() = default;
 
    void setup(const WxCryptSet& wxc) override;
 
    void solve(real time, variables_type& input) override;
 
protected:
    void alloc_fluxes();
    void zero_fluxes();
 
    void ex_kernel(const WxRange& face_scope_range,
                   std::vector<wxm::array::patch_array_t*>& input);
    void ldg_scheme(const int faceIndex,
                    const int local_face_node_index,
                    const int global_node_index_in,
                    const int global_node_index_out,
                    solverVariables_t* pSV,
                    std::vector<wxm::array::patch_array_t*>& input);
    void ip_scheme(const int faceIndex,
                   const int local_face_node_index,
                   const int global_node_index_in,
                   const int global_node_index_out,
                   solverVariables_t* pSV,
                   std::vector<wxm::array::patch_array_t*>& input);
 
    void in_kernel(const WxRange& element_scope_range,
                   std::vector<wxm::array::patch_array_t*>& input);
    void co_kernel(const WxRange& element_scope_range,
                   std::vector<wxm::array::patch_array_t*>& input);
 
    std::vector<size_t> _gradient_input_indexes;
    std::vector<size_t> _gradient_output_indexes;
 
    // Penalty for gradient solvers
    real _penalty_beta;
    real _penalty_eta;
 
    // Spatial basis set name
    std::string _basisSetName;
 
    // Spatial basis set
    WmBasisArraySet _basis_set;
 
    // WmUDGGeometry
    std::unique_ptr<WmUDGGeometry> _dg_geometry;
 
    // Fluxes
    std::vector<std::vector<std::vector<std::vector<real>>>> _numerical_fluxes;
    std::vector<std::vector<std::vector<real>>> _internal_fluxes;
 
    wxm::dfem::tools::scope_t _element_flux;
    wxm::dfem::tools::scope_t _face_flux;
 
    WxRange _chunk_size_element;
    WxRange _chunk_size_face;
 
private:
    gradient_interior_penalty_internal_approximation&
    operator=(const gradient_interior_penalty_internal_approximation& var) = delete;
 
    gradient_interior_penalty_internal_approximation(
        const gradient_interior_penalty_internal_approximation& var) = delete;
};
 
} // namespace gradient
} // namespace variable_adjuster
} // namespace dfem
} // namespace wxm
 
#endif // WMGRADIENT_INTERIOR_PENALTY_INTERNAL_APPROXIMATION_H