phase_space_gradient.h

Go to the documentation of this file.

#ifndef variable_adjuster_phase_space_gradient_h
#define variable_adjuster_phase_space_gradient_h
 
// WARPXM includes
#include "dfem/tools/tools.h"
#include "dfem/variable_adjusters/variable_adjuster.h"
#include "lib/geometry/phase_space_element/phase_space_element.h"
#include "lib/geometry/phase_space_geometry/phase_space_dg_geometry.h"
 
// STL includes
#include <memory>
 
namespace wxm
{
namespace dfem
{
namespace variable_adjuster
{
namespace phase_space_gradient
{
 
class PhaseSpaceGradient : public wxm::dfem::variable_adjuster::variable_adjuster_t
{
public:
    PhaseSpaceGradient() = default;
 
    void setup(const WxCryptSet& wxc) override;
    
    void solve(real time, variables_type& input) override;
 
    /* /\** */
    /*  * @brief Get a range of elements associated with the entire */
    /*  * patch (periphery+interior) */
    /*  *\/ */
    /* const WxRange& element_full_range() const override */
    /* { */
    /*     return _element_flux._full; */
    /* } */
 
protected:
    // input and output indices
    std::vector<size_t> _gradient_input_indexes;
    std::vector<size_t> _gradient_output_indexes;
    
    // gradient method
    std::string _gradient_method;
 
    // gradient solver penalties
    real _penalty_beta;
    real _penalty_eta;
 
private:
    // One function needs it
    real _cfl;
    
    // Element range
    wxm::dfem::tools::scope_t _element_flux;
    
    // Velocity space range
    std::unique_ptr<WxRange> _velocity_space_element_range;
    std::unique_ptr<WxRange> _velocity_space_interior_element_range;
    
    // Phase space geometry
    std::shared_ptr<geometry::phase_space::PhaseSpaceDGGeometry> _phase_space_dg_geometry;
 
    // Phase space element type
    std::unique_ptr<geometry::phase_space_element::PhaseSpaceElement> _phase_space_element;
 
    // assignment and copy operators?
    PhaseSpaceGradient& operator=(const PhaseSpaceGradient& var) = delete;
 
    PhaseSpaceGradient(const PhaseSpaceGradient& var) = delete;
};
} // namespace phase_space_gradient
} // namespace variable_adjuster
} // namespace dfem
} // namespace wxm
 
#endif // variable_adjuster_phase_space_gradient_h