calc_diffusion_coefficient.h

Go to the documentation of this file.

#ifndef variable_adjuster_calc_diffusion_coefficient_h
#define variable_adjuster_calc_diffusion_coefficient_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 diffusion_coefficient
{
 
class CalcDiffusionCoefficient : public wxm::dfem::variable_adjuster::variable_adjuster_t
{
public:
    CalcDiffusionCoefficient() = default;
    
    void setup(const WxCryptSet& wxc) override;
    
    void solve(real time, variables_type& input) override;
    
    real calcDeterminant();
    
    // *
    // *
    // */
    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;
    
    // 
    std::vector<size_t> _input_indexes;
    std::vector<size_t> _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;
    
    // basis set
    std::unique_ptr<WmBasisArraySet> _wavenumber_space_basis_set;
    
    // number of element nodes
    int _N_p_kx;
    int _N_p_ky;
    int _N_p_kz;
    
    // quadrature weights
    size_t _num_weights;
    std::vector<real> _weights_kx;
    std::vector<real> _weights_ky;
    std::vector<real> _weights_kz;
    
    // Jacobian determinant for k-space line subspaces
    real _k_space_determinant;
    
    // Phase space geometry
    std::shared_ptr<geometry::phase_space::PhaseSpaceDGGeometry> _phase_space_dg_geometry;
    
    // Wavenumber space geometry
    std::shared_ptr<geometry::phase_space::PhaseSpaceDGGeometry> _wavenumber_space_dg_geometry;
    
    CalcDiffusionCoefficient& operator=(const CalcDiffusionCoefficient& var) = delete;
};
}
}
}
}
 
#endif