divergence_integral_by_parts.h

Go to the documentation of this file.

#pragma once
 
// Wm includes
#include "dfem/tools/tools.h"
#include "dfem/variable_adjusters/variable_adjuster.h"
#include <memory>
 
namespace wxm
{
namespace dfem
{
namespace variable_adjuster
{
namespace auxiliary_variables
{
 
class DivergenceIntegralByParts : public wxm::dfem::variable_adjuster::variable_adjuster_t
{
public:
    DivergenceIntegralByParts() = default;
 
    void setup(const WxCryptSet& wxc) override;
 
    void solve(real time, variables_type& input) override;
 
protected:
    void rhs(real time,
             const WxRange& element_scope_range,
             std::vector<wxm::array::patch_array_t*> input);
 
    std::string _basisSetName;
 
    WmBasisArraySet _basis_set;
 
    std::unique_ptr<WmUDGGeometry> _dg_geometry;
 
    wxm::dfem::tools::scope_t _element_flux;
 
    real _penalty_beta;
    real _penalty_eta;
 
    std::string _source_quad_style;
 
    bool _has_fluxes;
    bool _has_sources;
    bool _sources_lgl;
    bool _sources_gq;
 
    std::vector<size_t> _rhs_idcs;
    std::vector<size_t> _all_idcs;
};
 
} // namespace auxiliary_variables
} // namespace variable_adjuster
} // namespace dfem
} // namespace wxm