general_maxwellian_two_stream.h

Go to the documentation of this file.

#ifndef wmapplication_function_kinetics_general_maxwellian_two_stream_h
#define wmapplication_function_kinetics_general_maxwellian_two_stream_h
 
// WARPXM includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace functions
{
namespace kinetics
{
class GeneralMaxwellianTwoStream : public WmApplication
{
public:
    GeneralMaxwellianTwoStream();
 
    ~GeneralMaxwellianTwoStream();
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _distribution;
    }
 
    void evaluate_function(const real* q,
                           const real* aux,
                           const solverVariables_t* pSV,
                           real* result) const override;
 
protected:
    std::vector<int> _distribution; // output variable
 
    // real _thermal_velocity;    // v_th
    // real _density0;            // n0
    // real _perturbation;        // A
    // real _wave_number;         // k
    // real _velocity_separation; // v'
 
 
    real _mass;
 
    real _n1;
    real _n2;   
    std::vector<real> _velocity1;
    std::vector<real> _velocity2;   
    real _pressure1;
    real _pressure2;    
    
private:
    GeneralMaxwellianTwoStream& operator=(const GeneralMaxwellianTwoStream& var);
    GeneralMaxwellianTwoStream(const GeneralMaxwellianTwoStream& var);
};
 
} // namespace kinetics
} // namespace functions
} // namespace apps
} // namespace wxm
 
#endif