velocity_space_moment_taker.h

Go to the documentation of this file.

#ifndef variable_adjuster_velocity_space_moment_taker_h
#define variable_adjuster_velocity_space_moment_taker_h
 
// WARPXM includes
#include "dfem/variable_adjusters/variable_adjuster.h"
#include "lib/geometry/phase_space_geometry/phase_space_dg_geometry.h"
#include "dfem/tools/tools.h"
 
// STL includes
#include <memory>
 
namespace wxm
{
namespace dfem
{
namespace variable_adjuster
{
namespace moment_taker
{
 
class VelocitySpaceMomentTaker : public wxm::dfem::variable_adjuster::variable_adjuster_t
{
public:
    VelocitySpaceMomentTaker() = default;
 
    void setup(const WxCryptSet& wxc) override;
 
    void solve(real time, variables_type& input) override;
    void solve_old(real time, variables_type& input);
    void solve_new(real time, variables_type& input);       
 
    std::vector<size_t> convertToOutputIndexes(const std::vector<int>& indexes) const;
 
protected:
private:
    std::vector<std::vector<size_t>> _app_to_out_idcs_func_eval;
 
    // The variables
    std::vector<size_t> _moment_input_indexes;
    std::vector<size_t> _moment_output_indexes;
 
    // checking of indexes
    std::vector<size_t> _out_idcs;
    std::vector<size_t> _all_idcs;
 
    // scopes
    wxm::dfem::tools::scope_t _element_flux;
 
    // velocity space basis set
    std::unique_ptr<WmBasisArraySet> _velocity_space_basis_set;
 
    // number of element nodes in each velocity space dimension
    int _N_p_vx;
    int _N_p_vy;
    int _N_p_vz;
 
    // weights
    size_t _num_weights;
    std::vector<real> _weights_vx;
    std::vector<real> _weights_vy;
    std::vector<real> _weights_vz;
 
    // velocity space element determinant
    // ASSUMPTION:  All velocity space elements are same size and so this just needs to be
    // calculated once
    real _v_space_determinant;
 
    // velocity space range
    std::unique_ptr<WxRange> _velocity_space_element_range;
    std::unique_ptr<WxRange> _full_velocity_space_element_range;
 
    // Phase space geometry
    std::shared_ptr<geometry::phase_space::PhaseSpaceDGGeometry> _phase_space_dg_geometry;
 
    VelocitySpaceMomentTaker& operator=(const VelocitySpaceMomentTaker& var) = delete;
 
    real calcDeterminant();
 
    void generateVelocitySpaceElementRange1VPositive();
 
    void generateVelocitySpaceElementRange1VNegative();
};
} // namespace moment_taker
} // namespace variable_adjuster
} // namespace dfem
} // namespace wxm
 
#endif